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Preface 


Beginning with Volume XX, the Deep Space Network Progress Report changed from 
the Technical Report 32- series to the Progress Report 42- series. The volume number 
continues the sequence of the preceding issues. Thus, Progress Report 42-20 is the 
twentieth volume of the Deep Space Network series, and is an uninterrupted follow-on to 
Technical Report 32-1526, Volume XIX. 

This report presents DSN progress in flight project support, tracking and data 
acquisition (TDA) research and technology, network engineering, hardware and software 
implementation, and operations. Each issue presents material in some, but not all, of the 
following categories in the order indicated. 

Description of the DSN 

Mission Support 

Ongoing Planetary/Interplanetary Flight Projects 
Advanced Flight Projects 

Radio Astronomy 

Special Projects 

Supporting Research and Technology 
Tracking and Ground-Based Navigation 
Communications— Spacecraft/Ground 

Station Control and Operations Technology • 

Network Control and Data Processing 

Network and Facility Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 
Quality Assurance 

Operations 

Network Operations 
Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Program Planning 
TDA Planning 

In each issue, the part entitled “Description of the DSN” describes the functions and 
facilities of the DSN and may report the current configuration of one of the seven DSN 
systems (Tracking, Telemetry, Command, Monitor & Control, Test & Training, Radio 
Science, and Very Long Baseline Interferometry). 

The work described in this report series is either performed or managed by the 
Tracking and Data Acquisition organization of JPL for NASA. 
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Network Functions and Facilities 

N. A. Renzetti 

Office of T racking and Data Acquisition 


The objectives, functions, and organization of the Deep Space Network are 
summarized; deep space station, ground communication, and network operations control 
capabilities are described. 


The Deep Space Network was established by the National 
Aeronautics and Space Administration (NASA) Office of 
Space Tracking and Data Systems and is under the system 
management and technical direction of the Jet Propulsion 
Laboratory (JPL). The network is designed for two-way 
communications with unmanned spacecraft traveling approxi- 
mately 16,000 km (10,000 miles) from Earth to the farthest 
planets and to the edge of our solar system. It has provided 
tracking and data acquisition support for the following NASA 
deep space exploration projects: Ranger, Surveyor, Mariner 
Venus 1962, Mariner Mars 1964, Mariner Venus 1967, Mariner 
Mars 1969, Mariner Mars 1971, and Mariner Venus-Mercury 
1973, for which JPL has been responsible for the project 
management, the development of the spacecraft, and the 
conduct of mission operations; Lunar Orbiter, for which the 
Langley Research Center carried out the project management, 
spacecraft development, and conduct of mission operations; 
Pioneer, for which Ames Research Center carried out the 
project management, spacecraft development, and conduct of 
mission operations; and Apollo, for which the Lyndon B. 
Johnson Space Center was the project center and the Deep 
Space Network supplemented the Manned Space Flight Net- 
work, which was managed by the Goddard Space Flight 
Center. The network is currently providing tracking and data 
acquisition support for Helios, a joint U.S./West German 
project; Viking, for which Langley Research Center provided 
the Lander spacecraft and project management until May, 


1978, at which time project management and mission opera- 
tions were transferred to JPL, and for which JPL provided the 
Orbiter spacecraft; Voyager, for which JPL provides project 
management, spacecraft development, and is conducting 
mission operations; and Pioneers, for which the Ames 
Research Center provides project management, spacecraft 
development, and conduct of mission operations. The network 
is adding new capability to meet the requirements of the 
Gahleo mission to Jupiter, for which JPL is providing the 
Orbiter spacecraft, and the Ames Research Center the 
probe. In addition, JPL will carry out the project management 
and the conduct of mission operations. 

The Deep Space Network (DSN) is one of two NASA 
networks. The other, the Spaceflight Tracking and Data 
Network (STDN), is under the system management and 
technical direction of the Goddard Space Flight Center 
(GSFC). Its function is to support manned and unmanned 
Earth-orbiting satellites. The Deep Space Network supports 
lunar, planetary, and interplanetary flight projects. 

From its inception, NASA has had the objective of 
conducting scientific investigations throughout the solar sys- 
tem. It was recognized that in order to meet this objective, 
significant supporting research and advanced technology devel- 
opment must be conducted in order to provide deep space 
telecommunications for science data return in a cost effective 
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manner. Therefore, the Network is continually evolved to keep 
pace with the state of the art of telecommunications and data 
handhng. It was also recognized early that close coordination 
would be needed between the requirements of the flight 
projects for data return and the capabilities needed in the 
Network. This close collaboration was effected by the appoint- 
ment of a Tracking and Data Systems Manager as part of the 
flight project team from the initiation of the project to the 
end of the mission. By this process, requirements were 
identified early enough to provide funding and implementa- 
tion in time for use by the flight project in its flight phase. 

As of July 1972, NASA undertook a change in the interface 
between the Network and the flight projects. Prior to that 
time, since 1 January 1964, in addition to consisting of the 
Deep Space Stations and the Ground Communications 
Facility, the Network had also included the mission control 
and computing facilities and provided the equipment in the 
mission support areas for the conduct of mission operations. 
The latter facilities were housed in a building at JPL known as 
the Space Flight Operations Facility (SFOF). The interface 
change was to accommodate a hardware interface between the 
support of the network operations control functions and those 
of the mission control and computing functions. This resulted 
in the flight projects assuming the cognizance of the large 
general-purpose digital computers which were used for both 
network processing and mission data processing. They also 
assumed cognizance of all of the equipment in the flight 
operations facility for display and communications necessary 
for the conduct of mission operations. The Network then 
undertook the development of hardware and computer soft- 
ware necessary to do its network operations control and 
monitor functions in separate computers. A characteristic of 
the new interface is that the Network provides direct data flow 
to and from the stations; namely, metric data, science and 
engineering telemetry, and such network monitor data as are 
useful to the flight project. This is done via appropriate ground 
communication equipment to mission operations centers, 
wherever they may be. 

The principal deliverables to the users of the Network are 
carried out by data system configurations as follows: 

(1) The DSN Tracking System generates radio metric data, 
i.e., angles, one- and two-way doppler and range, and 
transmits raw data to Mission Control. 

(2) The DSN Telemetry System receives, decodes, records, 
and retransmits engineering and scientific data gener- 
ated in the spacecraft to Mission Control. 

(3) The DSN Command System accepts spacecraft com- 
mands from Mission Control and transmits the com- 
mands via the Ground Communication Facility to a 


Deep Space Station. The commands are then radiated 
to the spacecraft in order to initiate spacecraft func- 
tions in flight. 

(4) The DSN Radio Science System generates radio science 
data, i.e., the frequency and amplitude of spacecraft 
transmitted signals affected by passage through media 
such as the solar corona, planetary atmospheres, and 
plantetary rings, and transmits this data to Mission 
Control. 

(5) The DSN Very Long Baseline Interferometry System 
generates time and frequency data to synchronize the 
clocks among the three Deep Space Communications 
complexes. It will generate universal time and polar 
motion and relative Deep Space Station locations as 
by-products of the primary data delivery function. 

The data system configurations supporting testing, training, 
and network operations control functions are as follows: 

(1) The DSN Monitor and Control System instruments, 
transmits, records, and displays those parameters of the 
DSN necessary to verify configuration and validate the 
Network. It provides the tools necessary for Network 
Operations personnel to control and monitor the Net- 
work and interface with flight project mission control 
personnel. 

(2) The DSN Test and Training System generates and 
controls simulated data to support development, test, 
training and fault isolation within the DSN. It partici- 
pates in mission simulation with flight projects. 

The capabilities needed to carry out the above functions 
have evolved in the following technical areas: 

(1) The Deep Space Stations, which are distributed around 
Earth and which, prior to 1964, formed part of the 
Deep Space Instrumentation Facility. The technology 
involved in equipping these stations is strongly related 
to the state of the art of telecommunications and 
flight-ground design considerations, and is almost com- 
pletely multimission in character. 

(2) The Ground Communications Facility provides the 
capability required for the transmission, reception, and 
monitoring of Earth-based, point-to-point communica- 
tions between the stations and the Network Operations 
Control Center at JPL, Pasadena, and to the JPL Mis- 
sion Operations Centers. Four communications dis- 
ciplines are provided: teletype, voice, high-speed, and 
wideband. The Ground Communications Facility uses 
the capabilities provided by common carriers through- 
out the world, engineered into an integrated system by 
Goddard Space Flight Center, and controlled from the 
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communications Center located in the Space Flight 
Operations Facility (Building 230) at JPL. 

The Network Operations Control Center is the functional 
entity for centralized operational control of the Network and 
interfaces with the users. It has two separable functional 
elements; namely, Network Operations Control and Network 
Data Processing. The functions of the Network Operations 
Control are: 

(1) Control and coordination of Network support to meet 
commitments to Network users. 

(2) Utilization of the Network data processing computing 
capability to generate all standards and limits required 
for Network operations. 

(3) Utilization of Network data processing computing 
capability to analyze and validate the performance of 
all Network systems. 

The personnel who carry out the above functions are located 
in the Space Flight Operations Facility, where mission opera- 


tions functions are carried out by certain flight projects. Net- 
work personnel are directed by an Operations Control Chief, 
The functions of the Network Data Processing are: 

(1) Processing of data used by Network Operations Control 
for control and analysis of the Network. 

(2) Display in the Network Operations Control Area of 
data processed in the Network Data Processing Area. 

(3) Interface with communications circuits for input to 
and output from the Network Data Processing Area. 

(4) Data logging and production of the intermediate data 
records. 

The personnel who carry out these functions are located 
approximately 200 meters from the Space Flight Operations 
Facility. The equipment consists of minicomputers for real- 
time data system monitoring, two XDS Sigma 5s, display, 
magnetic tape recorders, and appropriate interface equipment 
with the ground data communications. 
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Pioneer 1 1 Saturn Encounter Mission Support 

R. B. Miller 
Mission Support Office 


The Pioneer 11 spacecraft will fly by the Planet Saturn in August and September 
1979, The Pioneer 11 Saturn encounter activities are described, followed by the DSN 
plans foPsupporting the event. 


I. Pioneer 11 Saturn Encounter Description 

The Pioneer 11 spacecraft was launched on April 6, 1973, 
and flew by the planet Jupiter, with the closest approach on 
Decembers, 1974. After nearly 6-1/2 years of flight, the 
Pioneer 11 spacecraft will have a closest approach to the 
planet Saturn on September 1, 1979. At that time, the 
spacecraft will be over 10 astronomical units from Earth, and 
the radio transmission will take 2 hours and 53 minutes from 
Earth out to the spacecraft and back to Earth again. 

The near-encounter trajectory is pictured in Figs. 1 through 
3 (which were provided by the JPL Pioneer Navigation Team). 
Figure 1 is a view of Saturn from Earth during the closest 
approach with the spacecraft entering from the upper left and 
penetrating the ring plane at the point marked “inbound 
descending node,” which is at a distance of 2.87 Saturn radii 
(R^) from the center of Saturn. The ring plane penetration is 
symmetrical and takes place at closest approach ±2 hours at 
the same distance from the center of Saturn. The 2.S7-R^ 
distance was selected because this is the distance that Voyager 
2 will have to penetrate the ring plane if it is to go on to the 
planet Uranus. Therefore, Pioneer 11 will provide two data 
points to determine whether that is a safe distance to 
penetrate the ring plane. Periapsis occurs at 1.36 R^, just 
one-third of a radius of Saturn above the cloud tops, which is 


81 ,426 km from the center of the planet. Periapsis passage will 
occur at approximately 1800 GMT, ground observed time. 
Figure 2 shows the trajectory in a polar view and includes the 
orbit of the inner moon Mimas, and the visible ring. The line 
marked “occultation zones” shows the location of the 
spacecraft as it passes into the shadow of Saturn; this location 
is essentially the same as the point where Earth is occulted by 
the planet. Notice that periapsis occurs just before enter 
occultation, right under the inner visible ring (referred to as 
the C-ring). Figure 3 is a scale drawing (except for ring 
thickness), of Saturn, with rings edge on, which rotates with 
the trajectory to show the spacecraft path with respect to the 
ring plane and the disk of Saturn. At the inbound and 
outbound crossing points, labeled “a” and “c”, the spacecraft 
is at the distance that duplicates Voyager 2’s Uranus trajec- 
tory. This will also be approximately the periapsis distance for 
the Voyager 2 mission. When the spacecraft passes under the 
outer edge of the A-ring, it will be only 2,000 km from the 
ring plane, while at the closest approach point, labeled “b”, 
the spacecraft will be approximately 9,000 km below the ring 
plane. The ring plane penetration times of closest approach 
±2 hours also corresponds to the closest distance for which it 
is expected that Pioneer 1 1 will be able to successfully image 
the planet. This is because the relative motion of the target 
body will be too great for the spin scan imaging system. 
However, while under the ring plane, several imaging single- 
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spin scans of the ring plane will be performed to gain, 
hopefully, information about other smaller-scale features of 
the rings, such as smaller gaps than can be detected by 
Earth-based observation. 

The encounter time period is defined as August 3 to 
October 2, 1979. In this time period, approximately 15,000 
commands will be transmitted to the spacecraft, all of which 
will be acted upon at the time of receipt. The majority of 
these commands are required for the operation of the Imaging 
Photopolarimeter, To understand the operation of this instru- 
ment and why it requires such an extensive amount of 
commanding, see Ref. 1 . The Earth occultation will last from 
closest approach to E + 79 minutes. 

The imaging will exceed Earth-based resolution (which is 
approximately 1 ,200 km) at ±6 days around closest approach. 
The best resolution of Saturn is expected to be on the order of 
80 to 100 km achieved at E ± 2 hours. The rings will fill the 
field of view at 39 hours before closest approach, and the disk 
of the planet will fill the field of view 33 hours before closest 
approach. The trajectory allows for good imaging after closest 
approach. Approximately 50 images better than Earth-based 
are expected, with 20 of them being before and 30 after 
closest approach. Several of the satellites will be observed, 
both in imaging and polarimetry modes, with the best target 
being Titan. Closest approach to Titan will be 356,000 km at 
25-1/2 hours after closest approach. At this distance, the 
resolution is expected to be about 180 km. 

II. Voyager’s Expectations for the 
Pioneer 1 1 Saturn Encounter 

There are several areas in which the scientists involved in 
the sequence planning activities of Voyager hope to benefit 
from the data returned by Pioneer 1 1 . First and foremost is to 
determine if 2.87 is a safe distance to penetrate the ring 
plane. If Pioneer 1 1 is not able to safely pass through the ring 
plane at this distance, it is quite likely that the Uranus option 
of Voyager 2 will be deleted. 

There are several areas in which Voyager hopes to benefit in 
the general category of characterizing the fields and particle 
environment. Earth-based observation indicates that the fields 
and particle environment around Saturn is one-tenth or less as 
intense as it is for Jupiter; however, there is great uncertainty 
in the Earth-based observations. It is hoped that the Pioneer 
1 1 data will enable the determination of the magnitude and 
orientation of the magnetic field around Saturn. In particular. 
Voyager is interested in the orientation of the magnetic 
equator with respect to the planet’s equator in order to 
properly time Voyager’s fields and particle roll maneuvers. 
These are maneuvers where the entire Voyager spacecraft is 


rolled in order to act as a “scan platform” for the fields and 
particle instruments at key times in passing through Saturn’s 
magnetosphere. It is desired to have these roll maneuvers 
symetrically on either side of the magnetic equator and one 
right on the magnetic equator. Another item is to know 
whether the corotating magnetosphere of Saturn extends as far 
as the moon Titan. If the corotating magnetosphere extends as 
far as Titan, there will be a requirement to observe the 
magnetospheric wake behind Titan that will prevent the 
spacecraft from executing a roll maneuver in order to look 
back at Titan after closest approach to that moon. Character- 
izing the intensity and energy levels of the fields and particle 
environment will help Voyager determine what special pre- 
cautions may be necessary to keep the spacecraft and its 
instruments safe, as well as help in determining optimum gain 
settings for the fields and particle instruments. 

If the Pioneer 11 imaging is able to identify any unique 
visual features of Saturn, Voyager would use this information 
to refine the imaging plans. In particular, the Voyager team 
hopes that the Pioneer 11 imaging of Titan will enable a 
decrease in the uncertainties of the radius of that body. 
Current uncertainties in the radius of Titan are ±300 km; it is 
hoped that the Pioneer 1 1 data will reduce that uncertainty to 
on the order of 100 km. This would enable Voyager to 
produce a more efficient mosaicing plan of the moon Titan, as 
well as to assure the return of a limb picture with one-third as 
many pictures dedicated to that purpose. 

Lastly, Pioneer 1 1’s polarimetry of Saturn and its satellites 
may be the only polarimetry into the distant future, since this 
function has failed on the Voyager 1 spacecraft and may not 
be used on the Voyager 2 spacecraft because of the fear that 
the rotating filter wheels might stick in an undesirable 
position. 

III. DSN Support Plans 

Even though the Pioneer 1 1 Saturn encounter will involve 
the same level of activity as the Pioneer 10 and 11 Jupiter 
encounters, very limited special training on the part of the 
DSN is necessary to prepare for the event because it will look 
to the DSN very similar to the configurations and level of 
activity required for supporting the on-going Pioneer Venus 
orbital operations. The training will, therefore, concentrate on 
the special functions that are added to support the Saturn 
event. 

Only a single significant implementation in addition to the 
configuration supplied for the Pioneer 10 and 11 Jupiter 
encounters was necessary for the Pioneer 1 1 Saturn encounter. 
This implementation was needed to push the telecommunica- 
tions performance to the absolute limit available with current 
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technology: necessitated by attempting a Saturn encounter 
with a spacecraft that was designed for the Jupiter distance. 
The specific implementation was the installation of low-noise 
S-band masers at Australia and Spain; these installations were 
completed on May 6, 1979, The masers give an improvement 
in performance of approximately 0.8 dB when operated in a 
diplex mode (transmit and receive simultaneously). A similar 
low-noise maser had already existed at the Goldstone 64-meter 
station. An additional implementation done on a short-term 
basis is the installation of a programmable receiver at Station 
12 to accommodate the high doppler rates that will be 
experienced at the periapsis passage. 

In addition, to gain additional performance for ±7 days 
around closest approach (which covers the time period when 
Earth-based resolution is exceeded for the imaging instru- 
ment), the listen-only mode will be utilized at the 64-meter 
sites; this will give an additional 0.7 dB performance improve- 
ment. Because of the extensive command activity required, a 
second smaller station will have to be scheduled with each 
64-meter site for this 14-day period. The stations that will 
support during these 14 days are: DSSs 12 and 14 (34- and 
64-meter stations at Goldstone), DSSs 44 and 43 (Honey- 
suckle and Tidbinbilla in Australia), and DSSs 62 and 63 
(Cebreros and Robledo near Madrid, Spain). 

An R&D real-time arraying system, which has been installed 
at Goldstone on a demonstration basis for the Voyager 2 
Jupiter encounter, will be utilized for the 14 days around 
closest approach. This system will receive the spacecraft signal 
at the 64-meter site and 34-meter site simultaneously, with a 
microwave transmission of the 34-meter data over to the 
64-meter site where the signals are automatically phase- 
matched and added to achieve a performance improvement 
that is expected to be greater than one-half a dB. 

Against a Project requirement for continuous 64-meter 
coverage for closest approach ±30 days, the actual coverage 


negotiated between Voyager, Pioneer Venus, and Pioneer 1 1 
will be 10 to 19 hours a day from 30 to 15 days before closest 
approach, and then continuous coverage from that point in 
time until 30 days after closest approach. 

There will be some special training of each shift at all 
stations to support extensive, contingency manual command- 
ing required for the Saturn encounter. Since the imaging 
sequence is executed by the receipt of real-time commands, 
the possibility of a loss of communications circuits or a 
computer outage at the Ames Research Center, where the 
commands originate, is covered by plans that allow uninter- 
rupted command capability for the stations by voice instruc- 
tion from Ames, The level of commanding required is much 
greater than is usually used for manual commanding, so 
practice exercises have been scheduled to cover all affected 
shifts. 

At least two arraying tests will be scheduled to measure the 
actual performance improvement for the Pioneer 11 bit rate 
and coding type. In addition, it is expected that one or two 
periapsis passage training exercises will be, scheduled for 
Goldstone and Madrid because of the complexity of the 
operations required during that time period. Complex receiver 
tuning will be necessary to maintain receivers in lock because 
of the high doppler rates that will be experienced. This is the 
reason that the programmable receiver had to be installed at 
Station 1 2 so that that station could retain its receivers in lock 
to support the arraying. Programmable receivers already exist 
at the 64-meter sites. There will also be complex requirements 
to support radio science experiments, and to operate the 
occultation data assembly equipment and associated special 
recording functions. Finally, to guard the supercritical data, 
particularly at ring plane crossing and exit occultation, an 
activity is underway to investigate the feasibility of utilizing 
the precarrier detection telemetry recovery technique that was 
implemented for the Pioneer Venus Multiprobe entry event; 
this would provide a back-up analog recording in the event 
receiver or telemetry system lock was lost. 
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Fig. 1. Pioneer 11 Saturn flyby - balanced ring plane crossings at 
2.87 view from Earth 
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Fig. 2. Pioneer 11 Saturn flyby ~ balanced ring plane crossings at 
2.87 Rg, north polar view 
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The Gyrotron: A High-Frequency Microwave Amplifier 

A. Kupiszewski 

Radio Frequency and Microwave Subsystems Section 
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for future generations of millimeter high-power uplinks to spacecraft and planetary radar 
transmitters. Basic electron-electromagnetic field interaction theory for RF power gain is 
explained^ and the starting point for general analytical methods leading to detailed design 
results is presented. 


\. Introduction 

The next generation of spacecraft uplinks and planetary 
radar transmitters will probably operate at higher frequencies 
in the 20- to 35-GHz range. Present crowding at S- and 
X-bands as well as increased resolution and ranging available at 
higher frequencies make such a shift desirable. However, a 
cursory examination of readily available microwave amplifiers 
shows that there are no conventional high-power (>10kW) 
tubes useful for such purposes. Recently, much effort has been 
directed towards development of these higher frequency 
amplifiers. One such device of high potential is the gyrotron, a 
cyclotron resonance maser utilizing magnetic coupling be- 
tween cavity or waveguide fields close to cutoff and relativistic 
electrons gyrating in large dc magnetic fields. Its basic 
advantage over klystrons and TWTs is that the electronic 
interaction volume does not scale down with wave length, 
hence high power can be obtained without the conventional 
heat transfer problems inherent in small interactive volumes. 

The history of the gyrotron is in itself quite interesting. 
Theory for the device was developed as far back as 1958—59 
by R. Q. Twiss and J. Schneider (Ref. 1), but power output 
levels were then measured in milliwatts. Interest in the U.S. 


lagged, and funding was shifted elsewhere. In the U.S.S.R., 
however, research continued unabated, and some impressive 
results were achieved, for example 22 kW CW at 2 mm 
(Ref. 2). Due to such success, attention in the U.S. has again 
focused on this class of amplifier with contemplated applica- 
tions including plasma heating and mm-wave radar. 

II. Description 

Gyrotrons may be operated over a wide range of frequen- 
cies from 2 to 235 GHz, depending only on selection of cavity 
and magnetic field strength. Figure 1 shows a typical arrange- 
ment for a gyrotron. Electrons from an annular cathode biased 
at several thousand volts form a hollow beam. (Hollow, 
circular beams are used in the interest of efficiency, since there 
is no RF component on the axis of symmetry. Breakup is 
inhibited by the short drift space and relatively high voltages 
involved, Ref. 3). The beam is then compressed by an 
axisymmetric magnetic field according to the adiabatic invari- 
ant v\jB = constant, where is the Larmor orbit velocity and 
B is the magnetic field strength. This is a magnetic mirror 
effect similar to that used to confine plasmas, where perpen- 
dicular, azimuthal electron energy is increased at the expense 
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of logitudinal drift energy. In the interaction space, which has 
a gently varying cross section to maximize efficiency, the 
electrons are guided by highly uniform fields. Magnetic 
coupling to the RF cavity fields is responsible for the 
microwave amplification, and, after completing interaction, 
the electrons settle on an extended collector surface in a 
region of weaker magnetic field. 

In these respects, the gyrotron is just like klystrons and 
TWTs. An electron beam provides energy for microwave 
amplification. Spent electrons impact on an extended, water- 
cooled collector surface, which is the bulkiest, heaviest part of 
the tube. There the similarity ends. Instead of longitudinal 
bunching for microwave amplification, the gyrotron uses 
cyclotron resonance and azimuthal bunching. 

One of the most important characteristics of gyrotrons is 
the high efficiencies theoretically obtainable. With proper gun 
design and good electron optics, the axial symmetry favors 
circumstances where all electrons interact with the RF fields 
under identical conditions. High efficiencies are contingent 
upon such relative uniformity. Losses in the magnetic mirrors 
and the interaction space are minimized due to solenoid 
design, which ideally reduces magnetic field variations to 
within a few tenths of one percent of maximum field strength. 
An example of a high-efficiency tube design with all of the 
aforementioned characteristics is the gyrotron TWT designed 
by Chu, Drobot, Granatstein, and Seftor (Ref. 4). The 
important parameters given are an output power of 342.5 kW 
CW with 51.0 percent calculated efficiency for a beam voltage 
and current of 70.82 kV and 9.48 A, respectively, and a 
magnetic field strength of 17.87 kG. The output frequency is 
35 GHz with a maximum gain of 20 dB and a *3 dB 
bandwidth of 910 MHz. 


ill. Interaction Mechanism 

To gain an understanding of the interaction mechanism, 
including theoretical results for bandwidth and efficiency, it is 
necessary to employ classical calculations based on Vlasov 
theory. The scope of such an endeavor, however, exceeds the 
aims of this introductory report, so only an overview of the 
process for arriving at a condition for amplification in a 
cylindrical cavity gyrotron will be given. More details will be 
addressed in future articles. 

Starting with the Vlasov Equation (a continuity equation 
relating spacial changes in the electron distribution function to 
variations in time), one has the following: 

|7+vV/+^(E + vXB) • V/= 0 (1) 


where 

/ is the electron distribution function 
V is the velocity vector 
q is the charge on the electron 
E is the electric field 
B is the magnetic field 
m is the electron mass 

The distribution function and the cavity field can be split into 

/ = /o(x,v)+/j (x,v,/) (2) 

E=Ej (the electrons have already been 
accelerated before entering the cav- 
ity) 

B = Bq + (Bq is the dc magnetic field) 

(3) 

Using perturbation theory and linearizing to first order, the 
expression for the perturbed distribution is obtained. 

= ^(Ej + vXBj)- (4) 

And, for the TEqj ^ mode in cylindrical geometry, 

- ^11 

Bj = rE^ ^11^ 

- zE^—J^ {kj) sin z sin (5) 

Ej = 0 E^J co% (6) 

The time averaged power absorbed by the electrons is given by 
{P ^ ={Nq / av v/, • E, > . (7) 

abs time ^ *^1 1 time ^ ^ 
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Further integration yields the following: 




(14) 


M >,• 

, abs time 


q^El-np {k^a) v\ 


2mk^^ V 


3 

II 


2k, V, 


(8) 


Although the process outlined above gives detailed, correct 
results for the lineshape function, it is easier to understand the 
amplification mechanism if it is viewed quantum mechanically. 
Power absorbed by the electrons traversing the cavity can be 
expressed (Refs. 5 and 6) as 


P = N a((jj p ,,“CO -p ,) (15) 

n,n + l ^n,n+l ^ 


where is the RF field amplitude 

Q is the cyclotron frequency in the rest frame 

CO is the frequency in radians of the cavity field 

N(r, 0) is the electron density 

c is the speed of light in vacuum 

V|| is the drift velocity, is the parallel wave 
number 


where N is the number of electrons 

1i is Plank’s constant divided by 2ti 

(j^n is the frequency in radians of an upward transition 
from level n to n + I 

co^ „_i frequency in radians of a downward transi- 

tion from level n to n - 1 

M + 1 respective transition probabilities 


is the Larmor orbit velocity, k^ is the perpendicu- 
lar wave number 

a is the cavity radius 


p{k^a) = 



de f rdrN(r,e)j\{k^r) 

Jq 


(9) 


For free-free, higher-order transitions (i.e., n very large), it 
follows that p„ „ + i ^ Pn n-\' Amplification, i.e., P<0, 
occurs only if ^ mildly 

relativistic electron beam ('^lOkV), such a condition can be 
easily achieved due to 'Lorentz corrections to the cyclotron 
frequency. Looking at the energy levels, one finds that 
differences can be expressed as Lorentz corrected percentages 
of the cyclotron frequency. 


Q (X) = X- 
^ 2GdX 


ky 


( 10 ) 


cos 


ttX 


G(Ji2C) = lineshape function 


(1 - 


( 11 ) 


X = 


CO - 12 

k,v. 


( 12 ) 




k,cvj 


(13) 


By inspection, it can be seen that the condition for amplifica- 
tion is 


Energy] 




— “ ilco 




/ 


2|\l/2 


The lower energy levels have a smaller Vy than the higher 
levels, which produce an unequal energy spacing between 
levels. This unequal spacing accounts for the emission exceed- 
ing absorption. It should be noted, however, that pursuit of 
this line of reasoning leads to Lorentzian lineshapes and not 
G^(J0 j which is confirmed by experiments (Ref. 7). Hence, 
present work in this field is based almost exclusively on 
classical calculations utilizing relativistic corrections. 
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IV. Summary 

Gyrotrons are cyclotron resonance masers that show great 
promise in becoming the next generation of deep space 
transmitters due to high efficiencies at K-band frequencies 
with high-power outputs. Similarities to klystrons enhance the 
basic appeal of gyrotrons. Handling and maintenance technolo- 


gies should not differ much except in the area concerning 
cyrogenically-cooled superconducting magnets. Further re- 
search continues to produce impressive results with respect to 
power and efficiency at ever higher frequencies. Future reports 
on this subject will address new advances, efficiency calcula- 
tions, and potential integration problems with existing sys- 
tems. 
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Two-Station narrow-band /IVLBI requires phase connections between consecutive 
scans. This article presents an efficient computer-aided scheme for this purpose. This 
scheme is an iteration process alternating between a grand fit on many scans and integer 
quantization of the phase-shift cycles to be assigned to the scans. Only linear simultaneous 
equations of a few unknowns need to be solved. A simulation analysis indicates that 
faultless phase connection can be expected when there is no localized systematic noise. 
When systematic noise of moderate level exists, the possible incorrect phase connection 
can be detected and corrected for by comparing the connected phases from the two 
alternating observations, after removing the residual diurnal effects. 


I. Introduction 

Narrow-band AVLBI (Differential Very Long Baseline 
Interferometry) has been considered as one of the possible 
accurate spacecraft tracking systems. A demonstration plan 
has been made for Voyager spacecraft during Jupiter encoun- 
ters in 1979 (Ref. 1). As is well known, narrow-band AVLBI 
requires a long pass (3-4 hr) of tracking so that its diurnal 
signature can resolve, in two orthogonal directions, the angu- 
lar separation between spacecraft (S/C) and an angularly 
nearby extragalactic radio source (EGRS). On the other hand, 
with 2-station mode, which is most likely the case, S/C and 
EGRS observations are taken alternately and gaps in each 
data stream are inevitable. The resulting interferometric 
phase of each continuous data segment (scan) has a 2nn 
ambiguity with different n for different segments. To recover 
the correct diurnal signature of each data stream, the inter- 
ferometric phases of the corresponding scans have to be con- 
nected with correct 2mr phase shift assigned to each scan. 


Such phase connection, if done manually by “eyeballing” 
the variations of the phases, is not only time consuming and 
tedious but may result in incorrect integer n since the phase 
change across the gap may be as large as tens of cycles. 

This article presents an efficient computer-aided phase- 
connection scheme. A simulation analysis shows that the 
scheme works faultlessly when there is no systematic noise 
such as clock drift or localized media effects. When such 
localized fluctuation exists, the scheme may assign incorrect 
integer cycle adjustments to some scans. This can easily be 
detected and corrected for by a phase validation scheme. 

It should be pointed out that for a successful phase connec- 
tion, either manually or by a computer, the effect of the 
earth’s spin should be modeled out and only the “residual 
phases” need be connected. 
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II. The Grand Fit 

When manually connecting the phases, one would fit each 
scan by a straight line and then try to connect two consecu- 
tive scans at one time by shifting one of them up or down by 
2 a27t such that they look continuous. This appears no problem 
only when the two fitted straight lines have the same slopes; 
ambiguity arises when the slopes are different, as may happen 
in practice. Higher-order fit may seem to resolve this problem. 
However, an independent higher-order fit on each scan is 
prone to error due to a limited number of noise-affected data 
points in a scan. Consequently, a grand fit on as many scans 
as possible simultaneously is desirable for a computer-aided 
phase-connection scheme. 


To minimize (2) its partial derivatives with respect to 

flj, , , . y and Xj are set to zero. With some 

manipulations these become 


M j 

E E 


m = 0 /■= 1 



a 


m 




EE 



X. = 
! 



hi 1=1 


(3) 


It is well known that most, if not all, fitting (optimization) 
algorithms require the solved-for parameters to be continuous 
over the region of interest. However, the integers n of concern 
are discrete numbers. A grand fit with such constraint 
imposed^ becomes nonlinear and annoying. The present 
scheme does the grand fit free from such constraint and allows 
the parameters n to assume any value. An iterative process is 
then performed to successively quantize these values of n into 
integers. The iterative process will be discussed in Section III. 


and 




(2nl.)x. = = 

i=l 


(4) 


Let the measured residual phase of the ith data point in the 
yth scan be denoted as let there be J scans to be included 
in a single grand fit and Ij data points in the /th scan. Then for 
a grand fit of order M one needs to minimize the following 
function with respect to and Xj: 


Equations (3) and (4) are M + 7 simultaneous linear equations 
for the Af + / unknowns and x-. A further reduction of this 
system of linear equations can be done by substituting each 
Xj from (4) into (3). This yields 


J h 


M 


hi i=l \ m = 0 


E 


m = 0 



where t-j is the time. The symbol Xy has been used in place of 
n to denote its role as a continuous variable. 

Since only relative phases are needed we can arbitrarily 
set Xj = 0. Then (1) is replaced by 



Xi 1 <t>i, + E 


/=! 




/=2 L/=l 


" i.J IJ 


/ = E 


/-I 




J /. 

* E E 

h 2 /= 1 



/J 




r = 0,1,.,.,M 


( 2 ) 


Hence, a system of Af + 1 linear equations for is resulted. 
This can be written in the following matrix form: 


^For instance, a penalty function (1 - cos^'” tta:) can be added onto 

the function to be minimized (cf, Eq. 2 ). [ 5 ]{a}={^} ( 6 ) 
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where [ ] denotes a square matrix and { } a column matrix. 

The elements of [s] and {b } are 

S =S = 

r,m m.r 


PA 


+ 


/=! 


7=2 L/=l 



(7) 


and 


Z <,1 ^,M ^ E 


i=l 


7 = 2 L/=l 


^7 


Z-r 1.7 ^*,7 



( 8 ) 


In practice, M = 2 or 3 will suffice and the solutions of 
are straightforward. Once the coefficients have been deter- 
mined, the phase shift parameters Xj are calculated according 
to (4): 




M 

/ \ 

■ 



E 

_m =0 

fe'«) 

m ^i,j 

i= 1 

,/ = 2,3,. 



(9) 


III. Iteration Process of Integer Quantization 
for xj 

The phase-shift parameters Xj determined by a grand fit in 
the preceding section will in general differ from integers. If 
the differences are all small, one may conceivably set them 
equal to the integers they approximate. However, if the 
differences are not all small, one may hesitate to equate x^ to 
an integer. However, the following iteration bears out the 
method. 


From (1) it is obvious that any variation in Xj from correct 
integers induces corresponding changes in the fitted coeffi- 
cients . Let the change in each be m = 0, 1 , . . . , M. 
Then the variation in x^ can be expressed, according to (9), as 



( 10 ) 


This implies that the variation Xj from the correct integers 
increases with j. In other words, the first few Xy are much 
closer to the correct integers than the remaining Xy. Hence, 
one can comfortably set the first few Xj to the nearest integers. 
With these integers fixed, another grand fit is performed to 
improve and the remaining Xj. The first few improved Xj 
will now be closer to the correct integers and can be quantized 
with greater confidence. This process is repeated until all x- 
are quantized into integers. 

After each grand fit, the criterion of setting the first x^ to 
the nearby integer can be more relaxed, as this Xj is more 
likely to be very nearly an integer. On the other hand it is 
more probable for the remaining Xy’s to miss the correct 
integer values by an amount greater than 1/2 and the criterion 
should be made more stringent. In the proposed scheme the 
allowable variation is chosen to be 0.4 for the first x- after 
each grand fit; for each of the remaining Xy, both Xj and Xy_ ^ 
must have a variation < 0.1 . 


IV. Simulation Analysis 

The flow chart of the phase-connection scheme is shown in 
Fig. I . The order of the polymonial to be fitted is assigned. In 
most cases a second- or third-order polynomial will be suffi- 
cient. The simulated phase data are generated with the follow- 
ing parameters: 

data density = 1 per minute 

gap width"= W= SLT\ integer multiple of a minute 

number of data points per scan = /, a variable parameter 

data noise = a, a variable parameter 

pass length = 3 hours 

number of scans = 7 = (3 X 60)I(W + f) rounded to the next 
lowest integer 

= 0.5 rad/s 

< 7 ^ =0,-5 X 10"*, -1 X lO""^ rad/s^ 

=0,5 X 10”^ rad/s^ 

3 

(f) . = a o - 2n.n with /i. to be such that 

^/,7 m i,j 7 7 

m = l 

.<27T. 

Both M = 2 and M = 3 are tried for the fit. The following 
inferences are drawn from the simulation: 

(1) While higher order M for the polynomial to be fitted 
may be chosen even when “ ^3 “ 0 in the simulated 
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phase data, a lower-order polynomial is less susceptible 
to data noise. Hence one may start with a higher-order 
polynomial and decrease the order upon failure in 
phase connection. 

(2) As few as 2 data points per scan are allowable for a 
successful phase connection provided the gap width 
and the data noise are reasonably srnall. 

(3) Given a data density and a gap-width-to-data-span ratio 
the success in phase connection degrades as data noise 
a increases but seems independent of the gap width, at 
least for gap width <10 minutes. This relaxes the gap 
width limit for minimum data loss when time offset 
between S/C and EGRS observations is called for 
(Ref. 2). 

(4) The maximum data noise above which phase connec- 
tion may fail decreases with increasing gap-width-to- 
data-span ratio. Figure 2 is an example with 5-minute 
data span (per scan) over three hours. When there is no 
localized fluctuation due to systematic noise, a typical 
AVLBI pass^ will fall well within the convergence 
range, as shown in the figure. 


V. Grouping of Scans in a Single Grand Fit 

As pointed out in Section II, it is desirable to include as 
many scans as possible in a single grand fit. However, a prob- 
lem arises when sizable localized phase fluctuations exist 
among scans included in a single grand fit: A low-order poly- 
nomial may lose track of these fluctuations while a polynomial 
of higher order may degrade the convergence of the iteration 
process. Hence an appropriate grouping of scans in a single 
fit is essential. 

A criterion to determine the grouping of scans is the change 
in slope of the straight lines fitted to the consecutive scans. 
Two consecutive scans are to be grouped together if the slopes 
of their fitted straight lines satisfy 


\Sf - S ._ J I - (of + of ^ I mHz (1 1 ) 


where and Oj are the slope and its uncertainty of scan /. 
The number of scans / in a group increases until (11) fails. 
When this number of scans J is determined, the order of 
polynomial is selected to be 


M = J, if y<5 

M = 5, if J>S 


( 12 ) 


^S/N = 5 corresponding to cr = 0.2 rad, and gap width /data span = 1.2. 


The order of fit M is to be decreased by 1 upon failure of con- 
vergence in the iteration and a new fit is tried. If the iteration 
does not converge for ail M > 2 the last scan in the group is 
removed from the fit, the largest M according to (12) is 
selected and the process repeats. 

It is obvious that phase connection is redundant for a group 
with fewer than two scans. When the number of scans reduces 
to 2, due to the failure of either (11) or in the iteration pro- 
cess, the slope test of (1 1) is to be ignored. Also, the allowable 
deviation of the fitted Xj away from an integer is to be relaxed 
to 0.5 to exclude any failure in the iteration process. The 
actual deviation of each fitted Xj from its quantized value 
can serve as an indication of the degree of confidence in the 
resulting phase connection. 


VI. Validation of Connected Phases 

The phase-connection scheme discussed in the preceding 
sections works faultlessly when there is no localized phase 
fluctuations. Because of irregularities in transmission media 
(ionosphere in particular) and instability of electronic signal 
path, localized fluctuations do exist. When such fluctuations 
are sizable, not only do we need to divide the scans into 
several groups in connecting phases, but also some of the scans 
may be assigned incorrect cycle adjustments. To guard against 
such danger, the following phase validation process will be 
necessary. 

Since the residual phases from S/C and EGRS are each 
independently connected, an approach of phase validation is 
to bring the two connected phases together and examine 
whether their variation is continuous over the pass. However, 
such examination is indicative only if the two residual phases 
have identical residual model errors. The residual model 
errors in frequency offset are identical since common fre- 
quency standards are used for both signals. Baseline error is 
not identical even though the same baseline is used, owing to 
different baseline projections onto the planes of sky in the two 
different directions of the two sources. Angular position (right 
ascension and declination) error is also different. Since both 
baseline error and angular position error appear in the form of 
residual diurnal variation, they can be easily removed by 
fitting A cos r + 5 sin t to each of S/C and EGRS 
residual phases, where co^ is the earth’s spinning rate. In 
practice, two more terms, Ct + D, are needed to remove the 
residual frequency offset and phase offset. This is a simple 
linear least-squares fit with four degrees of freedom.^ 


^Alternatively, a simultaneous fit for S/C and EGRS phases with 
common frequency offset term can be adopted. The degrees of free- 
dom become 7. 
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After removing the residual model errors, the S/C and 
EGRS phases are plotted on the same graph. This can be done 
handily by print-plotting on standard computer print. Now 
comes the time for inspection. Unfortunately, this has to be 
done by eyeballing at the present time. This is best illustrated 
by an example of the actual AVLBI pass. Figure 3 shows the 
connected phases of Voyager 1 and OJ 287 with their residual 
model errors removed. Two consecutive scans from the same 
source (S/C or quasar) are considered “continuous” if the 
phase change across the gap is consistent with the phase 
change in the intervening scan from the other source. A quick 
glance over the graph reveals that there are discontinuities 
between scans 2 and 3 of S/C phase and between the last two 
scans of quasar phase. When these two 1 -cycle adjustments 
are made, the phases after removing the updated residual 
model errors are shown in Fig. 4. No obvious discontinuity 
remains and the adjustments are validated. Normally only one 
or two iterations of adjustment are needed to eliminate all 
offensive discontinuities. 

There has been concern that when a clock without ex- 
tremely high stability is used, phase connection may fail. By 
the above phase validation scheme, incorrectly connected 
phases can still be detected and corrected as long as the phase 


drift due to clock instability is continuous (i.e., no phase 
jumps occur). Figure 5 is an example where rubidium clock 
was used at one end of the baseline. A discontinuity is detected 
between the last second and third scans of the S/C phase. 
After a 1 -cycle adjustment the phases are shown in Fig. 6. 
Again, no obvious discontinuity remains. It should be noted 
that the gross effect of the clock drift has been absorbed by 
the fitting functions and only the residual “random walk” 
effect remains. 


VII. Conclusions 

A computer-aided, phase -connection scheme has been 
developed. A simulation analysis shows that, when there is 
no sizable localized phase fluctuation, this scheme connects 
VLBI residual phases faultlessly. When there are localized 
phase fluctuations the scheme may assign incorrect cycle 
adjustments to some of the scans. This can easily be detected 
and corrected for by a phase validation scheme. Currently 
this validation process requires an eyeballing inspection. A 
means to validate the connected phases automatically is cur- 
rently being sought. 
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This article compares the presently used JPL hydrogen maser suspended atomic 
storage bulb to a new, rigid, single-plane mounted bulb. The new bulb incorporates three 
major design changes: 

(1) Mounting design. 

(2) Alterations to the collimator. 

(3) Decrease in mass. 

These design changes are expected to increase the long-term stability of the frequency 
standard by reducing its sensitivity to vibration and thermal effects. 


I. Introduction 

In June of 1976, an experiment to test the gravitational 
redshift using a hydrogen maser in a suborbital rocket probe 
was successfully performed jointly by the Smithsonian Astro- 
physical Observatory (SAO) and NASA’s Marshall Space FHght 
center (Ref. 1). The design of the space maser was the result of 
a program begun under NASA sponsorship in 1962 (Ref. 2), 
and one of the results of this design effort was the single-plane 
mounted atomic hydrogen storage bulb. This mount consisted 
of a solid quartz collar with one rim fused to the quartz 
storage bulb and the other fixed to the RF cavity base plate. 
By lightening and refining this basic collar design, SAO 
designed the tetrahedral spoke support system that is used in 
the SAO VLG -1 1 today (Ref. 3). 

The success of the solid-mounted SAO storage bulb design 
prompted JPL to test this concept in one of two prototype 
masers that was disassembled for repair. The rebuilt prototype 
masers (designated and P 3 ) will be used as reference 


standards in the Frequency Timing System (FTS) test facility 
at JPL. The rigid mount bulb to be tested is shown in Fig. 1 . 


II. Two-Point Suspension Storage Bulb 

In the present JPL Hydrogen Maser design, the storage 
bulbs are suspended in the microwave cavity by a two-point 
mounting system (Fig. 2). The upper support is a quartz rod 
fixed to the center line of the storage bulb. This rod extends 
through the cavity frequency adjustment plate into a spring- 
loaded mounting fixture. The spring pressure on the rod holds 
the storage bulb securely against the cavity base plate. The 
lower support is a Teflon^ plug that is fixed to the cavity base 
plate and inserted into the neck of the storage bulb. The press 
fit of the Teflon plug into the neck of the storage bulb pro- 


* Teflon is a Dupont name for tetrafluorethylene fluorinated ethylene 
propolene copolymer. 
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vides support on a radial and transverse plane. This plug is 
center bored so that it also serves as a collimator for the stor- 
age bulb. One difficulty with this mounting scheme is that the 
bulb moves if the maser physics package receives a slight jolt. 
The magnitude of the displacement is determined by the clear- 
ance required in the slip joint of the spring mount, and it is 
this clearance that allows the storage bulb to rock slightly on 
the Teflon plug. Although the displacement of the storage 
bulb is small (0.0025 cm, 0.001 inches), it is sufficient to 
change the output frequency by approximately a part in 
10”^^. Since it is practically impossible to avoid some form of 
physical shock at different times during maser operation, this 
shift poses a real problem. 

The collimator consists of a center-bored Teflon plug with a 
0.953-cm (0.375-inch) diameter orifice. The dimensions of the 
collimator are primary factors in determining the atomic 
storage time in the microwave cavity. The equation for storage 
time for the special case of a spherical storage bulb with a 
cylindrical collimator is: 


where: 

is the storage time 

C is the length of the collimator 

b is the radius of the collimator 

a is the radius of the storage bulb 

V is the mean velocity of the hydrogen atoms entering 
the storage bulb 

Since the collimator is part of the storage area, it has an effect 
on the wall shift. In fact, since the collimator is a solid plug of 
Teflon, it has a different wall shift than the storage bulb, 
which is coated with Teflon, The equation for wall shift in a 
storage area coated with a single material (Ref. 4) is 

Aw = ^ (2) 

0 

where : 

Aoj is the shift in line frequency 
0 is the phase shift in atomic wave function per collision 
is the mean time between collisions 


III. Rigid Single Plane Mount Bulb 

Comparing Figs. 2 and 3, it can be seen that with the 
removal of the support rod, there is no direct connection 
between the frequency adjustment plate and the rigid single- 
plane mount bulb. All support is through the base plate. Also, 
it can be seen in Fig. 3 that the collimator is no longer used for 
support. Freed from this duty, the collimator can be made a 
permanent part of the storage bulb. This improves the 
characteristics of the wall shift since the entire storage bulb 
substrate is quartz, coated with a continuous film of Teflon. 
Due to the rigid mounting design, resettability of the maser is 
improved since there is no chance of changing 1) the alignment 
of the collimator with respect to the storage bulb, and 2) the 
quartz (dielectric) position within the cavity electromagnetic 
field. 

Following the SAG design, the collimator has been length- 
ened and its inside diameter decreased to 0.559 cm (0,220 
inch). From Eq. (1), it is seen that this will increase the storage 
time by a factor of five, approximately. This increase in 
storage time will manifest itself as an increase in atomic line Q, 
From the equation for cavity pulling (Ref. 4): 

Aw = Aw^-^ (3) 

where: 

Aco is the shift in the output frequency 

Aco^ is the shift in the cavity resonant frequency 
is the atomic line g 
is the cavity Q 

it is seen that the higher the line Q, the smaller the resulting 
shift in the master output frequency (i.e., the output 
frequency is more resistant to physical changes in the cavity). 

Preliminary test results have shown a factor of 2 increase in 
line Q using this new bulb configuration. 

One additional modification has been incorporated into the 
rigid mount bulb. The two-point suspension storage bulbs 
weigh between 280 and 300 grams. Small thermal fluctuations 
(0.00 1°C) present in the maser can change the dielectric 
constant of this mass. Since the storage bulb resides in a region 
of high electromagnetic fields in the resonant cavity, these 
changes can affect the output frequency of the maser. By 
reducing the amount of quartz to 230 grams, it is hoped to 
reduce this frequency shift. 
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IV. Conclusion 

The P 3 prototype maser is now under construction with the 
rigid storage bulb mount. The following modifications are 
necessary to fit the new bulb to the JPL maser: 

( 1 ) The material of the cavity bottom plate has been 
changed from aluminum to quartz (this modification 


also improves the thermal properties of the microwave 
cavity by decreasing the coefficient of thermal expan- 
sion of the base plate). 

(2) The thickness of the base plate has been increased. 

(3) A base support groove has been added to the quartz 
plate (Section A- 1 , Fig. 3). 
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Fig. 1. JPL hydrogen maser 
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A generalized modulation scheme, which includes minimum shift keying (MSK) and 
staggered quadriphase shift keying (SQPSK) as special cases, is analyzed. The general 
modulator can be realized as a one-input, two-output sequence transducer whose outputs 
select the carrier signal for each baud. This form of the modulator has the practical 
advantage of not requiring any RF filtering since there is no actual mixing of the carriers 
with the modulating signals. 

It is shown that the optimum demodulator (whether hard-decision or soft-decision) 
always can make its decisions from the received waveform over two bauds when the 
interference is additive white Gaussian noise, thus generalizing a well-known result for 
hard-decision demodulation of MSK and SQPSK signals. The power spectra of MSK and 
SQPSK signals are derived to isolate the role played by coherency between the 
modulating signals and the carriers. 


I. Introduction 

In this report, we study a generalized modulation scheme 
that includes minimum shift keying (MSK) and staggered 
quadriphase shift keying (SQPSK) modulation as special cases. 
The generalized modulator, whose form has certain practically 
advantageous features, is introduced in Section II as a specific 
one-input, two-output sequence transducer whose outputs 
select the carrier signal for each baud. The resulting modulated 
signal is described by means of a waveform “trellis” associated 
with the “trellis” of the sequence transducer. The particular 
choice of carrier signals to obtain SQPSK and MSK modula- 
tion are specified. 


In Section III, we use the waveform trellis of the modulator 
in a “magic genie” argument to prove that the optimum 
demodulator, whether hard-decision or soft-decision, can 
restrict its observations to two consecutive bauds when the 
interference is additive white Gaussian noise. This generalizes a 
well-known fact about MSK and QPSK modulation. 

In Section IV, we derive the power spectra for both MSK 
and SQPSK signals, paying special attention to the required 
relationships between the carriers and the modulating signals. 
Finally, in Section V, we make some concluding remarks and 
give some suggestions for future research. 
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il. The Modulator 

A. Sequence Transducer and Trellis 

Let d = • ' • d 2 , ' ' ' hQ \hQ data sequence that 

constitutes the modulator input, where each d- has value +1 
(corresponding to a binary 0) or - 1 (corresponding to a binary 
1). The heart of our modulator structure will be the sequence 
transducer, shown in Fig. 1 , whose output sequences x and y 
are given by 

x.=^{d. + d._^) (la) 

and 

yi=\‘^i^d.-d._^) (lb) 

for < I < + oo, where a is the alternating sequence defined 
by 

a. = (-!)'"» (2) 

From Eqs. (1) and (2), we see that the components of the 
output sequences x and y take values in the set {-*1, 0, + 1}, 
i.e., these sequences are ternary-valued. Moreover, for each /, 
either x- = 0 (corresponding to d-=^ ^/-i) or “ 0 (corre- 
sponding to d- = d-_^), but not both. Thus x- and_y. can (and 
soon will) be used to amplitude-moderate the “quadrature 
components” of a carrier in such a way that one and only one 
component will be present for each i. 


labels on this path specify the resulting output sequences x 
and y. 

The trellis of Fig. 1 is, of course, closely akin to the 
“trellis” introduced by Forney to represent convolutional 
codes (Ref. 1). The only essential difference is that the 
“sections” between the nodes at each depth are identical in 
convolutional code trellises because the encoder is a time- 
invariant transducer. The sequence transducer of Fig, 1, 
however, is time-varying because of the effect of the alternat- 
ing sequence a, which has period 2. In fact, the trellis of Fig. 2 
is periodic with period 2 in the sense that the section 
beginning at time instant i is the same as that beginning at time 
instant / + 2 for all i. 

We remark that the sequence transducer shown in Fig. 3 is 
equivalent to that in Fig. 1, as is easily seen from the facts that 
d\ = aid^ and that = -a^d^_^, since a^ = 

Note that if d is a random data sequence, i.e., if the random 
variables d^ are statistically independent and identically dis- 
tributed (i.i.d.) with P{d^ = +l) = P(J^-= -1)= 1/2, then d' is 
also a random data sequence. From Fig. 3, we can see that the 
sequences x and y will have identical statistics in this case. 

B. Modulator and Trellis 

We are now ready to introduce our proposed modulator, 
whose structure is shown in Fig, 4. Here, c^{t) and Cy{t) are 
carrier waveforms that will be specified later for specific 
modulation schemes. The baud length will be denoted as T. 
The function of the “RF selector switch” in Fig. 4 is to select 
one of the carrier inputs for transmission in each baud 
according to the following rule: 


A convenient way to display the input /output relationship 
of the sequence transducer of Fig. 1 is by means of the trellis 
shown in Fig. 2. The nodes of this trellis at depth i correspond 
to the possible states of the sequence transducer at time 
instant i. From Eq. (1), we see that this state can be chosen as 
the value of d ^_^ ; the nodes in Fig. 2 have been labeled with 
this choice of state. The branches leaving each state at depth / 
correspond to the possible state transitions, the upper branch 
being the transition when the input d^ equals “ 1 and the lower 
branch being the transition when the input d^ equals +1; the 
branches are labeled with the value of (x corresponding to 
the output pair for that transition. We see, for instance, that if 
the transducer is in state +1 at time instant 0, then the input 
Jq = - 1 will cause a transition to state ~ 1 at time instant 1 and 
the accompanying outputs will be (Xq, = (0, +1); similarly, 
the input = - 1 will also cause a transition from state +1 at 
time 1 to state - 1 at time 2, but the accompanying outputs 
will be (Xj, = (0, ~1). Thus, each path through the entire 
trellis corresponds to a particular input sequence d, and the 


5(0 = 



,iT<t<iT+T (3) 


We emphasize that no addition of RF signals and no RF 
filtering is required in the modulator, since there is no actual 
mixing of the carriers with modulating signals. 

Because of Eq. (3), we see that the trellis of Fig. 2 can be 
modified, as shown in Fig. 5, to show the input /output 
structure of the modulator of Fig. 4. The trellis in Fig. 5 
differs from that in Fig. 4 only in that the node depth is 
labeled by the time iT instead of by the time instant /, and 
that the transitions are now labeled with the value of s{t) in 
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the baud iT<t< (i + \)T rather than with (x,-, y>). Thus, each 
path through the entire trellis of Fig. 5 still corresponds to a 
particular input data sequence d, but the labels on the path 
now specify the resulting modulated signal s{t) from the 
modulator of Fig. 4. 


But we see from Eqs. (4) and (5) that the phase difference at 
time t = iT between and Cy{t) is (Aco)/T= iir, which is a 
multiple of 27 t when i is even; thus the phase of s(t) is 
continuous for all four of the above transitions. 


C. MSK Operation 

We show now that the modulator of Fig. 4 realizes binary 
minimum shift keying (MSK) modulation (Ref. 2) when the 


Suppose on the other hand that i is odd. We see now from 
Fig. 5 that the only carrier-switching transitions that are 
possible at r = iT for s(r) are: 


carriers are selected in the manner 


(i) 

from “C^(0 

to 





(ii) 

from +^^^(0 

to 

-Cyif), 

c^(t) = A sin 

(4a) 






l\ 2 / j 

(iii) 

from ’C (/) 

to 

+C;,(0, or 

and 



(iv) 

from +c (r) 

to 



Cy{t) = A sin 



(4b) 


where A and 0 are an arbitrary amplitude and arbitrary phase, 
respectively, where coq is the carrier center frequency, and 
where 


But, from Eqs. (4) and (5), we see that the phase difference at 
time t = iT between c^{t) and ~c^(0 is (Aco)rT+ tt = (/ + 1)7t, 
which is a multiple of 27 t when i is odd; thus again the phase of 
s(r) is continuous for all four possible transitions. This 
completes the demonstration that the modulation is indeed 
binary MSK. 


(A6j)r=7r (5) 

Since, in each baud, the modulated signal s{t) will be one of 
+c^(f), -c^{t\ ’^Cy{t) and ~Cy{t), it suffices to show that the 
phase of s(/^) is continuous at the transitions between bauds 
(i.e., at times t- iT for all i and for all choices of the data 
sequence d). 

From the trellis of Fig. 5, we see that the phase of s{t) is 
certainly continuous at r = iT if s(t) does not switch between 
the two carriers at time iT, since then s{t) will be the same one 
of +c^(Oj or ^Cyif) in the two bauds adjacent 

at time iT. It remains to show that the phase is continuous 
when s{t) switches from either +c^(0 or ~c^{t) to either +c^(0 
or -Cy{t\ or from either +c^(0 or ~Cy{t) to either +c^(0 or 
-c^{t), at r = iT. 

Suppose first that i is even. We see then from Fig. 5 that the 
only such carrier-switching transitions possible at r = iT for s(^) 
are: 


(i) 

from “C^(0 

to 


(ii) 

from +c^(0 

to 


(iii) 

from -Cy{t) 

to 

-cjf ) . or 

(iv) 

from +c^(f) 

to 



More precisely, the above argument shows that the modu- 
lated signal s(/^) is the same as for MSK modulation. However, 
in what is commonly considered to be MSK modulation, the 
data sequence controls the modulated waveform in the manner 
that when the data digit is +1 or -1, then the modulated 
signal in the baud is ±A sin [(coq + Acol2)t + 0] or ±/l sin 
[(coq - A(jjj2)t + 0 ] , respectively, when the appropriate sign is 
chosen to maintain phase continuity. In other words, the data 
digit directly controls the baud frequency. We see from 
Eq. (la), however, that x\ (which controls the sinusoid of 
frequency coq + Ao)/2) is nonzero if and only if d^ = and 
that y^ (which controls the sinusoid of frequency coq ~ Acj/2) 
is nonzero if and only if d^ ^ d^_ j . But the bit in the first 
difference of the binary data stream is 0 (corresponding to +1) 
when di = d^_ ^ and is 1 (corresponding to - 1 ) when d^ d^_ ^ . 
Thus, the modulator of Fig. 4 actually realizes differential 
MSK since the differences in the data sequence form the 
sequence that actually controls the baud frequency. However, 
we shall continue to say simply “MSK modulation,” rather 
than the more precise “differential MSK modulation,” to 
describe the modulation performed by the modulator of Fig. 4 
when the carriers are specified as in Eq. (6). 

D. SQPSK Operation 

We show now that the modulator of Fig. 4 realizes 
staggered-quadriphase shifpkeying (SQPSK) modulation 
(Ref. 3), also called offset-keyed quadriphase shift-keying 
(OKQPSK), when the carriers are selected in the manner 
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c^(t) = A cos (Wq/ + 6) (6a) 

and 

Cy(t) = A sin + d) (6b) 

where again >1 and d are an arbitrary amplitude and arbitrary 
phase, respectively. 

By definition, SQPSK is four-phase modulation in which 
the phase can change by either 0 or ±7t/2 between bauds, but 
never by ±n. (This is usually accomplished by modifying an 
ordinary QPSK modulator for a baud of length 27" in which 
two independent data sequences binary-antipodally modulate 
the quadrature carrier components so that one data sequence 
makes its transitions mid-way between the other’s transi- 
tions — the advantage over ordinary QPSK is a reduction of 
cochannel interference due to nonideal effects in demodula- 
tion.) But this restricted phase-changing is obviously achieved 
by the modulator of Fig. 4 since, as we see from Fig. 5, s(t) 
can never make a transition between +c^(r) and -c^(r) or 
between +c^(Q and “C^(0, which are the only transitions 
according to Eq. (6) for which the phase of s(t) would change 
by ±7T. 

Again, we remark that we have really shown here only that 
the modulated signal s(t) is the same as for SQPSK modula- 
tion. Just as was the case for MSK modulation, however, when 
the correspondence between the data sequence and the 
modulated signal is considered, one finds that the modulator 
of Fig. 4 actually realizes differential SQPSK. 


III. The Demodulator 

We now consider demodulation for the modulator of Fig. 4 
when the received signal is 

KO = s(t) + «(0 (7) 

where n(t) is additive white Gaussian noise (AWGN). We will 
show the rather remarkable fact that the hard-decision 
demodulator shown in Fig. 6 is optimum (in the sense of 
minimizing the probability of error in the decision for d^, 
for all 0 regardless of the choice of c^(t) and Cy(t) (i.e., these 
two carriers need have no special orthogonality properties in 
each baud) provided only that, for every i, c^{t) and Cy(t) have 
the same energy in the baud iT<t< iT + T (but the energy 
could depend on /) and that d is a random data sequence. 
Henceforth, we assume that this energy condition and data 
condition are satisfied. 


To demonstrate the optimality of the hard-decision demod- 
ulator of Fig. 6, we exploit the “magic genie” approach of 
Wozencraft and Jacobs (Ref. 4, p. 419). Suppose we wish to 
estimate d^ where i is even. Suppose further that the genie is 
kind enou^ to tell us both the state, a,-, of the modulator at 
time instant i and also the state, at time instant i + 2. If 
the genie says = +1 and a ^*+2 “ instance, we see from 

the trellis of Fig. 5 that d^ = +1, would imply s(/) = for 
iT < t <iT -^2T, whereas d.-- \ would imply s(t) = +c^(0 for 
r < iT + IT. Moreover, any permissible choice of s(t)y for 
t< iT and for t > iT 2T, when d^ = +1 is also permissible 
when d-= -1, since the only requirement is that the state be 
+1 at time instant i and again +1 at time instant i + 2. Thus, 
the decision problem for d^ (with genie’s help) reduces to 
deciding only whether c^(t) or Cy(t) was transmitted in the 
interval iT<t< iT + 2T. But this is the classical problem of 
deciding between two equally-likely equal-energy signals in the 
presence of AWGN, and the well-known (Ref. 4, pp, 238-239) 
optimum decision rule is: Choose d- = -\-\ if and only if 

J iT+2T fiT’^2T 

r{t)cjj)dt > I r{t)Cy{t)dt (8) 

iT J iT 

Defining 

^/T+r 

Xi = r(f)cjt)dt (9a) 

J iT 

and 

Yi = I r(j)Cy(t)dt (9b) 

J iT 

we see that Eq. (8) can be written as 

(>0) 

But now suppose instead that the genie had told us that 
Uy = “ 1 and Recalling that i is even, we see from the 

trellis of Fig. 5 that d^ = +1 corresponds to s(r) = -Cy(jf) for iT 
< t < iT+T and s{t)= c^(t) for iT-^T < t < /r+2r. 
Similarly, d> = - 1 corresponds to s(r) = -^^^(0 for iT<t < 
iT + T and ^(r) = Cy{t) for /T + T < r < /7" + 2T. Thus, the 
optimum (genie-aided) decision rule is: Choose c?,- = +1 if and 
only if 


r iT-«- T 
J iT 


r{t)Cy{t) dt + 


^iT+2r 

JiT+r 
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>- 


riT-\^T 

I ' 

JiT 


r(t)c(t) dt + 


J, 


iT+2T 


r+r 


r(t)Cy(t)dt (11) 


whether and Cy(t) are orthogonal over each baud as they 
are in both MSK and SQPSK. 


Using Eq. (9), we can write Eq. (1 1) as 


I /+! i 1+1 


( 12 ) 


which we see is precisely the same condition as Eq. (10)! 


Similar analyses for the case o^ = +1 and a ^+2 = ”1 and for 
the case “ 1 j 2 “ " 1 show that the optimum genie*aided 
decision rules are again: Choose = +1 if and only if Eq. (10) 
is satisfied. But the four cases considered exhaust the possible 
values for a,, and o^^ 2 ' conclude that we can exorcise the 
genie \ we have no use for his information since the optimum 
decision rule for d^, when i is even, is independent of his 
information and is just: Choose d- = if and only if Eq. (10) 
is satisfied. Note that this is precisely the rule for the 
demodulator in Fig. 6 since, for i even, ~a^ = +1 so that 

and j , where and are defined as 


A. 


X. + X 


i+i 


(13a) 


We have argued elsewhere (Refs. 5 and 6) that demodula- 
tors should be designed to maximize the cut-off rate, /?q, of 
the discrete channel created by the modulator, waveform 
channel and demodulator, rather than to minimize bit error 
probability for a hard-decision demodulator. Thus, we find it 
much more satisfactory than showing that the demodulator of 
Fig. 6 is the optimum hard-decision demodulator to show 
that: 


When A^ - B^ is taken as the output, the demodulator 
of Fig. 6 is optimum in the sense of maximizing the 
cut-off rate, of the discrete channel between the 
modulator input and demodulator output (and also in 
the sense of maximizing the capacity, C, of this 
channel). 

To prove this claim, we must show that the demodulator with 
output A^ - B^ preserves the likelihood ratio for the decision 
on d^, since any operation on r{t) reduces (and also C) 
unless and only unless this likelihood ratio is preserved 
(Ref. 5). 


= 03b) 

Next, suppose i is odd and the genie tells us that o-- +1 
and a, + 2 “ From Fig. 5, we see that d^ - +1 corresponds to 
5(f) = Cj^(f) for iT <: t < iT + 2T, whereas d^ = -l corresponds 
to s(t)= ’-Cy(t) for iT<t< iT + 2T, The same argument as 
before now shows the optimum genie-aided decision rule to 
be: Chose d^“ if and only if 


We begin by letting r, s and n be the vector representations 
of r(t), s{t) and «(f), respectively, in some appropriate 
Euclidean space (“signal space”)- Given that r is received and 
that a magic genie informs the receiver that d^ = 5y for all / i, 
the likelihood ratio for the decision on d^ is 

p(r\ d =H,d. = S foTj^i) 

A, =— ! (14) 

' p(i\d. = -\,d. = d. fox j¥=i) 


X. + X.^, 

I I+l 


>-Y - K , 

I 1+ 1 


where p is a conditional probability density function. Equa- 
tion (14) can be rewritten as 


Analysis of the three remaining choices of and a ^-^_2 shows 
that this same decision rule is also optimum. Thus, the genie 
can again be exorcised. We note further that this optimum 
decision rule is the rule for the demodulator in Fig. 6 since, for 
i odd, -a^ = - 1 so that A^ = X- + X^^ j and B^ = - F- - Y-^ ^ . 

This completes our proof that the hard-decision demodu- 
lator in Fig. 6 is optimum in the sense of minimizing 
Pr (d^ ¥=d^ when the noise is AWGN, when d is a random data 
sequence and, when for each /, the carriers c^(f) and Cy{f) have 
the same energy in the baud iT<t< iT + T. Note that the 
optimum decision d^ for d^ depends only on r(t) over the 
two-baud interval rT<f < iT + 2T; this is a well-known fact 
for both binary MSK and SQPSK, but we have now demon- 
strated that this “two-baud optimality” is independent of 


) 

A. = — 

P(r|s_,) 


(15) 


where j and s_ ^ are the signal space representations of the 
waveforms 5^.j(f) and 5_j(f) assumed by s{t) for the data 
sequence dj= 8- for /^/ with d^ equal to +1 and -1, 
respectively. Since r = s + n and the noise n{t) is AWGN, Eq. 
(15) becomes 




A. 

' 1 




( 16 ) 
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(21) 


where is the one-sided noise power spectral density. The 
assumption that, in each baud, Cjj(f) and c^(f) have the same 
energy implies that s+,(0 and s_j(r) have the same energy or, 
equivalently, |s+j = |s_jp. Thus,Eq. (16) becomes 


-s_,)] 

A. = e “ 

I 


(17) 


By the correspondence of “dot product” in signal space to 
correlation of waveforms, 


r-(s^, - s_j) = 



KO [s+i(0- (it 


(18) 


2 


TT (A.-B.) 


A. 


which shows that the output does indeed preserve the 

value of the likelihood ratio, as was to be shown. 

We note that, since the demodulator of Fig. 6 mthA^ - B- 
taken as the output preserves the likelihood ratio (and hence 
all the statistical information about d^), the optimum (in the 
sense of maximizing or C) soft-decision modulator with 
any specified number of demodulator output levels can be 
obtained simply by appropriately quantizing the output Z- (cf. 
Refs, 5 and 7 for details of this procedure). 


Because the data sequences yielding s^j(f) and s_^{t) differ 
only in the value of d^, we see from the trellis of Fig. 5 that 
1 (t) = j(r) for t < iT and t> iT + 2T so that Eq. (1 8) can 
be written as 

r • (s+, - s_,) = I* '■(0 Ki(^)- s_i(0] dt (19) 

JiT 


IV. Modulation Spectra 

A. General Considerations 

We now examine the spectral properties of the signal 
produced by the modulator of Fig. 4. 

With the sequence x of Eq. (la), we associate the pulsed 
waveform jc(r) defined by 


It remains to evaluate Eq. (19). From Eqs. (2) and (9) and 
the trellis of Fig. 5, we find that, for all four possible values of 
and d^+^, 


r iT+ T 
JiT 


r (t) , (0 - s_ j (0] dt = X.+ a. Y. 


and 


riT+2 

JiT+T 


r(t) [s+,(0- s_j(0] dt = 


1+1 1+1 /+1 


= X. 


/+1 


^/^/+1 


x{t) = X. iT<t<{i^\)T (19) 

In the same* way, we associate pulsed waveforms >^(^), d{f)^ and 
a{t) with the sequences y, d, and a, respectively. Note that a{t) 
is a “square wave” of period 2T. 

For purposes of spectral analysis, we need now to randomize 
the time origin so that all time signals under consideration 
become stationary random processes. We accomplish this by 
replacing all time functions [e.g., x(f) and s(0] with their 
translations by t^ seconds [e.g., x{t + Iq) and s{t + Tq)] where 
tQ is a random variable uniformly distributed over the interval 
0 < ^0 ^ statistically independent of the data pulse 

values. We will write B^{j) and 7 ^^(t), for instance, to denote 
the resulting autocorrelation functions 


This, together with Eq. (13), implies that Eq. (19) can be 
written as 


s_j) A. B. (20) 


r^(t) = E[x (r + r^ + r)x(r+ r^)] 
= E’[x(rQ + r)x(r^)] 


(which shows again that we can exorcise the genie, since >1,- 
and B^ can be formed without the genie’s help.) Finally, 
substitution of Eq. (20) into Eq. (17) gives 


R^ir) = £’[s(r + r^j +r)s(r + rQ)] 

= E’[5(rQ + r)5(/Q)] 
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respectively. It should be noted that, since all time signals are 
translated by seconds, this randomization of the time origin 
neither introduces nor removes any “coherence” between the 
modulating signals and the carriers. 

Note that the data signal d{t + is the familiar telegraph 
wave (with baud length T and amplitude ±1). Thus, is 

the triangular autocorrelation function shown in Fig. 7; and 
the power spectral density (the Fourier transform of ^^(t)) is 
the familiar 

s^if) = rsinc"* ifT) (20) 

where 

. V sin(7T2) 
sine (z) = ^ 

7TZ 

The spectrum is also shown in Fig. 7. 

Recalling Eqs. (la) and (lb), we see that 

x{t) =^[d{t) + d{t- T)] (21a) 

and 

y{t) = y a{t) [d (0 - d{t-T)] (21b) 


From Eqs. (21a) and (21b), we find by a simple calculation 
that 

= I R,(r ♦ n + i (tH fi (7 - r) 

= i (i) ca 


which further implies 

= Syif) = S^{2f) = rsinc^ (2/T) (23) 

These autocorrelation functions and spectra are shown in 
Fig. 8. Inasmuch as x{t) and y(t) are ternary- valued (0, +1 or 
-1), it is a bit surprising that their autocorrelation functions 
are of the same form as the binary-valued (±1) data signal d{t). 


Recall from our discussion in Section II about the 
sequences x and y that, for each /, one of and y- will be 0 
and the other will be +1 or -1. Thus, the random processes 
x(r+ Tq) and y{t'\- t^) are highly dependent. It is quite 
surprising then to find that the random processes x{t + r^) and 
y{t + ^q) are uncorrelated, i.e., that 

= E [jc(r + + T)y + r)] = o (24) 

as can be verified by a direct (but somewhat tedious) 
calculation from Eqs. (21a) and (21b), 

We now turn to the main task at hand, finding the spectrum 
of the modulated signal, i.e., of the random process s(r + t^). 
From Eq. (3), we see that 

s(r) = at(0 c^(t) + y{t) Cyif) (25) 

which is our starting point. To proceed further, we must 
consider specific carrier signals. One general note of caution is 
in order, however, at this point. Because of Eq. (24), it might 
seem that s(^ in Eq. (25) can be treated as the sum of two 
independently modulated signals. This is not true in general, 
however, because some types of coherence between c^{t) and 
Cy{t) can cause x{t + t^) c^{t + t^) to be correlated with j^(r + 
^o) though x(r+ /q) and.y(r+ t^) are 

uncorrelated. We shall soon see that this type of coherence 
between c^{t) and Cy{t) does not occur with SQPSK modula- 
tion, but does occur with MSK modulation. 

B. SQPSK Spectrum 

As we saw in Subsection II.D., the signal 5 (r) is an SQPSK 
modulated signal when the carriers are specified by Eqs. (6a) 
and (6b). From these equations and Eq. (25), we find in this 
case that 

s{t ^-t^) = A x{t + 1^) cos [wq {t + ?q) + 0] 

+ .4j(^(f + /jj)sin [Wp (? + ?q) + 0] (26a) 

It follows from Eq. (26a) that 

Ejj) = + 

= -^ R^{j) cos (cOjjT) + R^(t) cos (cOj^t) 

(27) 
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(Where the vanishing of the double carrier frequency terms 
under expectation strictly speaking requires T to be an integer 
multiple of the carrier period but is virtually the case 

whenever » I IT.) By virtue of Eq. (24), we see that the 
term resulting from the cross-correlation of the two carriers 
vanishes (i.e., Eq. (26) represents independent modulation of 
orthogonal carriers) so that Eq. (27) becomes 

R^(t) cos (oJjjT) (28) 

where we have also made use of Eq. (22). From Eqs. (28) and 
(23), it follows that 

5 (/) {sinc^ [2(C- /„) T] + sine' [Hf+Q T \ } 

(29) 

where is the carrier frequency in hertz. Equation 

(29) gives, of course, the well-known spectrum of SQPSK, 
which coincides with that of ordinary QPSK. The purposes of 
our deriving this result here are to emphasize the role played 
by the vanishing of R^y(r) and to point out that, because the 
phase angle 6 in Eq. (26) is arbitrary, no special phase 
relationship is required between the carriers and the modulat- 
ing signals. The spectrum S^(f) of Eq. (29) is sketched in 
Fig. 9. 

C. MSK Spectrum 

When the carriers in Eq. (25) are specified as in Eqs. (4a) 
and (4b), s(t + t^) becomes the MSK signal 



Applymg familiar trigonometric identities, we can rewrite 
Eq. (26b) as 

s(t = A [;ic(r + tg)+y(t + r^)] 

• cos (f + /„)] sin [Wq (t + ?(,) + e] 

+ A [x(r + /o)-j'(/ + ro)l 

• sin |^^(f + ?o)] ( 26 c) 


We are thus led to define the signals 

d^(t) = x(r)+j^(0 (30a) 

and 

d^(t) = x(r) - y(t) (30b) 

With the aid of Eqs. (21a) and (21b), we see from Eqs. (30a) 
and (30b) that 

! d{t), iT<t< (/ +1)7’ and i odd 

d(t - T), r < (j + 1) r and / even 

(30c) 

and 

( d{t - T), iT<t< (i +1)7’ and i odd 

d,(t) = 

( d{t), iT<t<(i+ l)T and i even 

(30d) 

In other words, d^(t) and d.^(t) are independent amplitude 
random data signals of baud length 2T, staggered so that the 
transitions in d^{f) occur at times iT with i odd whereas the 
transitions in d^it) occur at times iT with i even. The data 
values (+1 or -1) in d^(jt) correspond to those of d(t) in the 
bauds /T < t< (z + \) T with i odd, while the data values in 
d^{t) correspond to those of d{t) in the bauds /T < ^ < (/ + 1) 
T with i even. In particular, it follows that 

= /?^^(r) = 7?/r/2) (31) 

Substituting Eqs. (30a) and (30b) into Eq. (26c), we find 

+ e] + A d^it + sin (t + fj,)J 

cos[<o^(t + tQ) + e] (32) 

It follows from the statistical independence (and zero means) 
of the amplitudes of d^(t + r^) and d^(t + t^) that the 

cross-correlation between the two terms on the right in 
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Eq. (32) vanishes. Moreover, each of the two terms clearly has 
the same autocorrelation function so that 

= E [s(r + t^+r) s(t + fg)] 

= >1^ cos(co^t)^^(t) (33) 

where is the autocorrelation function of the random 

process 


m(t + = dj(r + t^) cos (^ + fp)j (34) 

We now proceed to find 

Em ^ ^ + /„ + t) w (/ + t^)] (35) 


where we have also used the fact that (Aco)r = tt. Taking 
Fourier transforms (recalling that R^(r) = R^(~t)) in 
Eq. (38) yields ^ 


T 

S (/) = ^ 


sine +sinc {ifT + jj 


sSinc 


2(/-^)r]t.inc[2(/tf)r 


(39) 


where 


^ = _L 

•' 27t 2T 


(40) 


= E[m{t^ + T)m(t^y\ 


Equation (39) is our desired expression for the spectrum of 
(36) m{t). 


It is convenient to take first the conditional expectation in Eq, (36) for a fixed t , which gives 


’ j w(r„ + r) m(t^) I (tg + r) {t^) cos (?„ + r)J cos J | 


COS 


(fo + cos ^ fo] ’ 0<fg<r-Tand7’<tQ<2r 


0, T-T<t<T 


(37) 


Recalling that is uniformly distributed over 0 < < 2T, we find that with the aid of Eq. (37) that 





t>2T 





0<t<2T 


(38) 
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( 42 ) 


From Eqs. (33) and (39), we have finally 


5(/)=i^'7’<:sinc 




+ sine 


+^>l2 7’<sincl2l/+/, 


[■( 


0 2 
2 


)r 




+ sine 




0 2 
2 


)r] 


(41) 


whieh is the well-known speetrum of MSK, but in a form 
perhaps more transparent than the usual one. This speetrum is 
sketehed in Fig. 9 for eomparison with that of SQPSK. 


We remark here that the phase angle 0 in Eq. (26e) is 
arbitrary. Thus, the MSK speetrum of Eq. (41) does not 
require any speeial phase relationship between the modulating 
signal and the earriers cj^i) and Cy{i). 

The MSK speetrum does, of eourse, depend on the 
eoherency between the modulating signals, x(r) and y{t), and 
the earriers, sin [(co^ + Aco/2) t d] and sin [(cj^ - Aco/2) t + 
0] , in Eq, (25). If there were no eross-correlation between the 
two terms on the right in Eq. (26b), then the autoeorrelation 
funetion of s(r + r^) would be simply 


1 


(r)eos CO + 


= \a^R^{t) |cos r + cos | 


and henee the power speetrum of S{t + would be 

+ sWc’ 


r |smc " [2 rj 

Comparing Eqs. (41) and (42), we see that the eross-eorrela- 
tion between the two terms on the right in Eq. (27) has the 
effeet of changing the spectrum from a sum of the squares of 
two sine functions (as in Eq. (42)) to the square of the sum of 
the same two sine funetions. Moreover, outside the eentral 
lobe of S^if) (i.e., outside the range - 3/4T < /< f^ + 
3/4 T), one of the sine funetions in Eq. (41) is positive and the 
other is negative for all /, as follows from the faet that their 
arguments differ by 2(A/) T=l. Thus, one of these sine 
functions always tends to cancel the other outside the central 
lobe, which is the reason that the MSK spectrum drops off 
much more rapidly than that of SQPSK outside the central 
lobe. 

We have already mentioned that the phase angle 6 between 
the carriers and the modulating signals has no effect on the 
MSK spectrum. In fact, more is true. The phase angle 6 in 
Eq. (26b) can even be time-varying, i.e., a random process 0{t 
+ r^), and the above analysis goes through unchanged provided 
only that 6(t) not change appreciably over one baud (say, by 
0.1 radians or less.) This latter condition is required to insure 
the vanishing of the expectation of the double carrier 
frequency terms when computing 


V. Concluding Remarks 

It is well-known (Ref. 8) that MSK modulation can be 
viewed as “offset QPSK” modulation in which the modulating 
signals are sinusoidal pulses in each band rather than rectangu- 
lar pulses — this follows immediately from Eq. (32) above. 
Simon (Ref. 9) has offered a generalization of MSK, based on 
this fact, by considering more general pulse shaping. 

By expanding on the left and using the usual trigonometric 
identities, one easily proves the identity 

sine + sine 

r 4 / 

(43) 

With the aid of Eq. (43), Eq. (39) can be rewritten as 


_ 1 1 + cos (27 tZ) 

2ir^ / 72 _J_V 


S (f) = — ^ ^ ( 4^/0 

TT^ (I6f 


(44) 
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which is the usual baseband form of the spectrum of MSK and 
shows that the spectrum drops off as 1//^ for large / (in 
contrast to that of QPSK which drops off only as 1//^). 

Our aim in this report has been to study MSK in particular 
from a new standpoint with the hope of isolating the 
fundamental structure that accounts for its spectral improve- 
ment over QPSK. The “sequence transducer” approach that 
we have taken appears amenable to new generalizations of 
MSK signalling by different choices of the sequence transducer 
and carriers. This seems to us to warrant further investigation. 

Finally, we remark that several recent papers, e.g., (Refs. 10 
and 1 1), have considered a “phase trellis” approach to MSK 


and other forms of continuous phase modulation. (Our 
approach is closely related to these but somewhat more 
general since we can admit arbitrary “carrier signals” when 
converting the trellis of the sequential transducer to a 
“waveform trellis.”) Anderson and Taylor (Ref. 11) in parti- 
cular show substantial power improvements (2 to 4 dB) over 
QPSK together with a smaller bandwidth. However, their 
modulation scheme introduces “convolutional code type” 
constraints into the modulated signal so that demodulation 
over a fixed number of bauds is not optimal (as it is for MSK). 
At this writing, it is not clear whether these schemes actually 
exploit the potential “coding gain” of simpler modulation 
formats and thus cannot be further encoded efficiently, or 
whether in fact they actually utilize bandwidth more effi- 
ciently than simpler modulation schemes. It seems to us that 
this is an important question to settle. 
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S.(f) = Ts5nc^ (fT) 



Fig. 7. The autocorrelation function (r) and the power spectral 
density S(§ (f) of the random data signal d (f + to) 


R^(r)=Ry(T) = l/2R^(T/2) 
1 1/2 


s (f) = SJf) = T line^ (2fT) 
* y 

It 


-2T 


2T 




-3/2T -V2T -1/2T -/2T 2/2T 3/2T 

Fig. 8. The autocorrelation functions and power spectra of the 
outputs X (f) and y (f) of the sequence transducer of Fig. 1 



Fig. 9. Power spectra of MSK and SQPSK showing attenuation relative to spectral peak at sidelobe center frequencies 
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Internal Noise of a Phase-Locked Receiver With a 
Loop-Control ied Synthesizer 

C. A. Greenhall 

Communications Systems Research Section 


We propose a local oscillator design that uses a digitally programmed frequency syn- 
thesizer instead of an analog VCO. The integral of the synthesizer input, the '^digital 
phase, '' is a convenient measure of integrated doppler. We examine the internal noise of 
such a receiver. At high carrier margin, the local oscillator phase noise equals that of the 
Block IV receiver, about 2 deg rms at S-band, whereas the digital phase noise is about 
0.5 deg rms. 


I. Introduction 

A simplified diagram of the carrier tracking function of the 
Block IV receiver is shown in Fig. 1 . To obtain the local oscil- 
lator output, the doubled output of the voltage-controlled 
crystal oscillator (VCXO) is mixed with the output of the 
Dana Synthesizer and multiplied up to the first heterodyne 
frequency at S-band or X-band. The synthesizer is pro- 
grammed to follow the previously estimated carrier frequency, 
and the difference between the estimated and actual carrier 
phases is made up by the loop-controlled VCXO. 

It has been proposed that the VCXO be replaced by a loop- 
controlled digital frequency synthesizer such as the Dana 
Synthesizer already used in the programmed oscillator. This 
idea leads to the local oscillator design shown in Fig. 2. The 
50-MHz reference clocks the synthesizer and is also mixed 
with the synthesizer output; the purpose of mixing here is to 
reduce the factor by which synthesizer phase noise is multi- 
plied. The gain applied to the digital loop error signal is chosen 
to make the loop transfer function essentially the same as that 
of the present Block IV receiver, for a chosen bandwidth 
setting. 


An important feature is the digital phase. At present, 
doppler phase is extracted by mixing, the local oscillator out- 
put with a signal from the exciter. In the proposed design, the 
frequency programmed into the synthesizer is integrated 
digitally to produce a phase that is in step with the contents of 
the synthesizer’s internal phase register. The synthesizer itself 
is a phase-locked loop that tries to keep its output phase in 
step with this internal register. Thus the digital phase can be 
expected to stay close to the local oscillator phase. 

The main purpose of this article is to compare the digital 
phase noise of the proposed receiver to the local oscillator 
phase noise of the present Block IV receiver. The models 
include only the “internal” phase noises due to the Block IV 
VCXO and the Dana Synthesizers; in particular, thermal noise 
and the jitter of the station reference timing signal are ignored. 
Effects of sampling and quantization in the proposed receiver 
are also ignored. 

We shall see that, for high carrier margin (at least 40 dB 
above threshold), most of the noise in the Block IV local 
oscillator comes from the Dana Synthesizer. Thus, the local 
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oscillator phase noises of the present and proposed receivers 
are almost equal, about 2 deg rms at S-band. The rms noise of 
the digital phase of the proposed receiver is about 0.5 deg, but 
grows very slowly as integration time increases. 


II. Properties of the Dana Synthesizer 

Since the Dana Synthesizer appears in both the present and 
the proposed local oscillators, we shall first set down our best 
estimates of its transfer function and phase noise spectrum. 
Assume that the reference timing signal has no jitter. Let v{f) 
be the frequency programmed into the synthesizer. (Actually, 
v{t) is sampled every 10 /is but we shall treat time as contin- 
uous.) Then the synthesizer output phase is 

27t — — v{t) + n^(t) 


where L^{s) is the phase -to-phase transfer function of the 
device, and n^{t) is the additive phase noise. 

A gross approximation of the transfer function is given by 






( 1 ) 


where = 165 jus. The one-sided fiducial bandwidth of the 
filter L^{s) is 0.75/r^ = 4500 Hz, which is large enough so that 
we may replace L^{s) by 1 when using the synthesizer as a 
loop-controlled oscillator in a narrow-band tracking loop. 


where 


\.SX 10~®. 

a 

This noise is nonstationary because of the 1// term at low 
frequencies. The noise variance for Fourier frequencies / 
above 1 Hz is 


= 2.3 X 10“* rad^ (3) 

Jitter in the reference is not included. Both and o^ have a 
±3 dB uncertainty. 

The following explains how these values were arrived at; 
the uninterested reader can skip to Section III. 

We first carried out a theoretical computation of the Dana 
Synthesizer phase noise spectrum. The model consisted of a 
VCO with 1//^ phase noise spectrum (Ref. 1), controlled by a 
second-order loop with perfect integrator. This yielded a 
synthesizer phase noise spectral density that behaves like / at 
low frequencies. Because two references (Refs. 1 and 2) report 
spectral densities approximately like l//at low frequencies, we 
were forced to abandon our model and rely wholly on an 
understanding of the synthesizer as a black box. 

The shape of the spectrum in (2) comes from a paper by 
G. Gillette of Dana Laboratories (Ref. 1); he plots the spectral 
density down to / = 1 Hz. Meyer and Sward of JPL (Ref. 2) 
measured a spectral density like 1//^’^^ for 4 Hz </< 80 Hz. 
We shall assume a 1// behavior down to zero frequency. This 
assumption may be dangerous. 


Our estimate of the one-sided spectral density of the phase 
noise in rad^/Hz is given by 


f ’ 

0</< 100 Hz 



Tto’ 

100 Hz </<10Hz 



u 

i 

10kHz</< 50 kHz 

f 

0, 

/> 50 kHz 


( 2 ) 


Because a frequency synthesizer derives its time base from 
a reference oscillator, whether internal or external, one must 
enquire whether published reports of synthesizer phase noise 
include the jitter of the reference. Meyer and Sward used a 
technique that cancels reference jitter. Dana’s manual for its 
synthesizer (Ref. 3) gives a method for measuring total phase 
noise variance (discussed below); this method also cancels 
reference jitter. We therefore assume that Gillette’s graph does 
not include reference phase noise. 

Since Refs. 1 and 2 deal with an old model of the Dana 
Synthesizer, we used Gillette’s graph to get the shape of the 
phase noise spectrum, but used phase noise specifications of a 
current model to obtain the normalizing constant Refer- 
ence 3 specifies a total phase noise of -54 dB = 4 X 10“^ rad^ 
in the band ±15 kHz, excluding/^ ±1 Hz, where is the 
programmed output frequency. Figure 3 shows in simplified 
form the setup used to verify this specification. The reference 
signal of synthesizer 2 is slaved to that of synthesizer 1 . 
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Suppose that the synthesizers have outputs 
x,(/) = cos (2n/|f + 0,(0 + 


x^{t) = cos + e^it) . 

Nvhere and are the programmed frequencies, B ^ (f) and 
B^{t) are independent, zero-mean phase noises, and is a 
constant phase offset. The rms voltmeter sees the signal 

= /: n/ 2 cos (2jr (/, - /j)/ + 0(0 + , 

where B{t) -B ^2(0> ^ is an unknown constant. The 

first thing we do in the test is to set /j “/2 “i* i Then the 
voltmeter reading is K. This is the reference level. Next, we 
set /i = /2 /> where 0.05 Hz </< 0.1 Hz, and observe the 

peak reading of the rms voltmeter, which now sees the signal 

y{t) = K \J2 [cos (27r// + 0^) cos 0(0 


- sin (27r/r + B^) sin B (0] . 


The phase noise specification is defined to be the square of the 
ratio of peak reading to reference level. Over the short run, say 
up to one minute, we can assume 0(0 1 rad; hence 

y{t)^K y/2 [cos(27T/r + 0Q)- 0(Osin(27r/r + 0^)] . 

The voltmeter, insensitive below 1 Hz (Fig. 3), does not 
respond to the first term. The second term causes the volt- 
meter reading to rise and fall every 1/(2/) seconds; the peak 
reading, K \fT a(0), occurs when the synthesizers are in 
quadrature. Dividing by the reference level K, we conclude 
from the published phase noise specification that 

2a^(0) = 2(a^(0,)+a^(0j)) = 4X 10"Sad^ 

Therefore, the larger of a^(0,) and 0^(0^) is between 10"® 
rad^ and 2 X 10"* rad^. We shall use 


Integrating (2) gives - 131.4 jV^, whence - 1.5 X 
10 "®. 

A more direct estimate of comes from the single side 
band (SSB) phase noise spectral density specification in Ref. 3. 
The phase noise in the band /q ±15 kHz, excluding /q ±10 Hz, 
is given as -90 dB/Hz maximum. The SSB spectral density, 
denoted by (0 in the literature, is asymptotically related 
to the baseband phase noise spectral density by 

provided /is large enough so that 



5^(2tt/V/' « 1 rad^ 


(Ref. 4). From (2) we set 5’^(207 t) = 10~^. Then NJ20 = 
10“^; = 2 X 10“®. This agrees well with the 1.5 X 10“® 

computed earlier. We settle on the latter figure with a ±3-dB 
uncertainty. 

Finally, the value of in (3) is obtained from by 
including the extra power from 15 kHz to 50 kHz. 


III. Internal Noise of Present Block IV 
Receiver 

The rms internal local oscillator phase noise at S-band is 
given in Table 1 as a function of design point bandwidth bj^Q 
and carrier margin. The total noise is the RSS of the noises 
from the VCXO and the Dana Synthesizer. At high carrier 
margin, most of the noise, about 1.7 deg, comes from the 
Dana Synthesizer and is almost independent of bandwidth. 
Most of the local oscillator noise is outside the loop passband. 
Derivations of these results follow. 

From Fig. 1 we see that the noise in the local oscillator 
phase 0, with the VCXO input shorted, is 


a 


2 

spec 


2X 10"* rad^ 


M + 2ny)) , 


To compute N we set 


r, 


SX27Tf)df = . 

^ ^ ^ spec 


where and are the phase noises of the Dana Synthesizer 
and the VCXO, and Af = 20 for S-band. The noise is 

d 

described in Section II. We shall assume that n has a one- 
sided spectral density A^^//^ ; the constant A^^ will be calibrated 
from receiver specifications. 
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When the phase-locked loop is allowed to control the local 
oscillator, its phase noise becomes 




where 


L(s) = 


1 + T^S 

1 + TjS + ’ 


( 4 ) 


r, = 0.75/»^„, 

and r depends on the one-sided threshold loop bandwidth 
bho the carrier margin. At threshold, r = 2. At high 

margin. 


r = 32 for 2^?^^ = 1, 10, 100 Hz 

= 18.5 for 2Z?^^ =3, 30, 300 Hz 

We are using the perfect integrator approximation of the loop 
filter, F{s) = (1 + t^s)I{t^s), 


The specifications of the Block IV Receiver-Exciter Subsys- 
tem give a rms phase jitter of 9 deg when T 2 = 1 .5 s, r = 2, 
M = 20. (This results from operating the VCXO in a phase- 
locked loop with these parameters and a strong signal, so that 
thermal noise is absent.) This calibrates and allows us to 
fill in the VCXO noise column of Table 1, which lists 2Mo^ at 
threshold and at high margin. When the receiver is switched to 
wide mode, is multiplied by 100, is divided by 10, and 
r is unchanged. Therefore, is multiplied by 10 and a is 
unchanged. 

Table 1 also lists the total rms internal noise a(^, where 


IV. Internal Noise of Proposed Receiver 

We now have to consider two phases, the local oscillator 
output phase B and the digital phase The rms noise a(B) is 
about 1.7 deg at S-band, the same as the Dana Synthesizer 
component of the noise in the present Block IV receiver. 
Since the VCXO is absent, o(6) does not rise as carrier margin 
decreases. Again, most of the 6 noise is outside the loop 
passband. 


The high-pass filter I - L (s), acting on cancels the 1// 
spectrum at low frequencies and passes most of the rest. 
Therefore, the rms contribution of the Dana Synthesizer to 
local oscillator phase noise is essentially 

Mo^ = 3.0 X 10“^ rad = 1.7 deg . 

We shall ignore the decrease at the higher loop bandwidths 
since it is at most about 0.1 deg. 

The variance of (1 - L (s)) is given on p. 32 of Ref. 5, 
It is convenient to recast those formulas as follows: 


On the other hand, most of the noise in the digital phase H 
is inside the loop passband. This noise turns out to be nonsta- 
tionary, with a spectral density like 1// at low frequencies. 
Thus it makes no sense to quote a single variance o^(6). 
Instead, we compute 

oHlT) var [?(f + T) - ?(f)] 




where D = I - 4jr. 


which is independent of t because S/dt is stationary. The 
factor 1/2 is needed to make a fair comparison with a(0), for 
if we apply the operator (•, r)to the stationary process 
B , we get 

{6, r) = ^ var (fi{t + 7} - e{T)) 

= R(0)-R{T) 

-»■ R(0) = a^(d) as , 
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where R is the autocovariance function of d. On the other For computation, we shall approximate by the Block 

hand, we shall see that IV transfer function L, given by (4). Then 


(?, as . 


Table 2 lists the rms noise a(0, T) at S-band as a function 
of r/r^ and the loop threshold bandwidth where = 
0.75/Z>^^, for high carrier margin. The values decrease sli^tly 
as margin decreases. Since d is “almost” stationary, a^(0, T) 
grows slowly with T, in fact, like log T. We must warn the 
reader, however, that this result holds only if the spectral 
density of the Dana Synthesizer phase noise behaves 

like 1// at least down to / = 1/(47). Derivations of these 
results follow. 


1 

yvarj 


{e. T) 


-f 


(I- cos uT)S^i<o)^ 


/: 


\ (1 - cos cor) \L (ico) P ^ 


Although the loop error signal of the proposed receiver is 
digital, we shall proceed with an analog analysis. A phase- 
locked receiver using the local oscillator of Fig. 2 has the 
closed-loop transfer function. 

1 + TS 

LAs) = 

1 + Tji + T\s'^l(rLjs)) 


where SJicS) is given by (2). This leads to the further approxi- 
mation 

a^(e, T) = M'^N^ (7, +/j) 

where 


where Lj^s), given by (1), is the transfer function of the Dana 
Synthesizer. The one-sided bandwidth of 7/^(5) is^bout 4500 
Hz. The internal noise in the local oscillator phase 6 is 



(1 “ cos u)T) \L(ico) 


du) 

Cl) 


( 6 ) 


(1 


(5) is the main portion, and 


where njjt) is the Dana Synthesize^ phase noise. As with the 
Block IV receiver, the rms value o{6) is approximately = 
1.7 deg for S-band. 

To obtain the internal noise in the digital phase 0, suppose 
that n^{t) is the only active influence on the loop. The fre- 
quency (Hz) programmed into the synthesizer is (1/27 tM) 
d'd/dt. The local oscillator phase (rad) is 


e = M 




•^2 = 200^ f (7) 

•'2007T 

is a correction that is significant only for 2Z)^^ = 100 Hz or 
300 Hz. In the Appendix, these integrals are evaluated in 
closed form in terms of the elementary functions and the 
exponential integral / (e^/z) dz. Numerical results are given in 
Table 2. 


= L^{s)e^Mn^ . 


Replacing d by (5), we get 


L^is) 


L,{s) 


Mn^ . 

a 


V. Conclusions 

We have proposed a design for a receiver local oscillator 
containing a loop-controlled synthesizer. The transfer function 
of the proposed receiver and its performance in the presence 
of thermal noise are essentially the same as those of the 
present Block IV receiver. The proposed design provides a 
convenient digital estimate of carrier phase, obtained by 
integrating the synthesizer input frequency. 
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For S-band and high carrier margin, the rms internal phase 1 .7 deg 

noises of the present (Block IV) and proposed receivers are 
approximately as follows: 

0.3 - 0.7 deg 

1.7 — 2.0 deg Present local oscillator phase noise. 

Increases as carrier margin decreases. 


Proposed local oscillator phase noise. 
Insensitive to carrier margin. 

Proposed digital phase noise. 

Decreases as margin decreases. 
Increases as integration time increases. 
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Appendix 

Evaluation of Two Integrals 


The integrals in question are/^ and given by (6) and (7). To evaluate we evaluate the integral 
We rewrite them as follows: 


•'o 


1 - cosax 


f(x)dx 


r“ 1 - cos tx 

Kt) = I —dx 

Jft a: (1 + x^) 


In terms of the exponential integrals 


■-f 

j. 


— I fix) dx 


^ t 


where 


m = 


1 + rx^ 


1 + (r - 2) + x'* 


Ei(0 




a = s/rT/T^, Xp = 200?; tJs/t. 


which are tabulated in Ref. 6. We use Raabe’s integrals 
(Ref. 7): 


From now on, we assume that r > 4. 


The partial fraction expansion of f{x) is 


" J " i p£j(0 + e-'Ei(o] . 


fix) 


^ ^2 


X^+Xj X^+Xj 


r ” X cos tx 1 r . .1 

I, it) = J dx = ^ [e'F, (0 - e-' Ei(oJ 


where 


(The integral form of / it) will not play a part here.) Since 


1 - rXj 1 - rx;: 

- "*2 = 

^2-^1 . 


^1 X ^ 


f 

Jo 


^( 0=1 lSM)du, 


( 8 ) 


x\ =y - 1 + ^ Vl - 4/r, x] = 1/x^ . 


Integrating /(x), we get 


4-1 it) = lit)- — , 

dtc ^ t 


we have, for some constant C, 


•'2 = 


x^ix\ - 1) 


^Xj- rjlan + 


^0^2 




I{t) = /^(0 + lnr + C, />0. 


As / 0, I{i) ^ 0 from (8). Limiting expressions for E^{t) 

and Ei(r) give 
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I^{t) + lnr -»-7 zst-*-0 , 


/j = (1 + k)I^(ax^) - kl^(a/x^) + 7 + (1 + 2A:) In + In a 


where y is Euler’s constant, 0.577216. Therefore, C= 7 . where 

For large t. 


I(t) = lnr + 7“ 1 /f^ + 0 (l/r^). 


k = 



+ 1 


- 1 


Finally, we have 


These formulas were spot-checked by numerical integration. 
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Table 1. RMS internal noise of Block iV receiver at S-band (total 
noise is RSS of VCXO noise with Dana synthesizer noise = 
1.7 deg) 


Hz 

Carrier 

margin 

VCXO noise, 
2A/a^, 
deg 

Total noise, 
a(0), 
deg 

Narrow 

Wide 

1 

10 

1 

9 

9 




0.9 

2.0 

3 

30 

1 

3.0 

3.5 



“ 

0.5 

1.8 

10 

100 

1 

0.9 

2.0 



eo 

0.09 

1.7 

30 

300 

1 

0.3 

1.8 




0.05 

1.7 


Table 2. RMS internal noise of digital phase d of proposed receiver 
at S-band as a function of integration time T (high carrier margin 
assumed) 


r/Tj,T 2 = 0.75/62,0 



1 

102 

lO'* 

10 ® 

1 

0.29 

0.42 

0.52 

0.60 

3 

0.28 

0.41 

0.51 

0.59 

10 

0-29 

0.42 

0.52 

0.60 

30 

0.29 

0.42 

0.52 

0.60 

100 

0.42 

0.52 

0.60 

0.67 

300 

0.50 

0.58 

0.66 

0.72 


o&, T), deg 



S-BAND 
INPUT ■ 


X-BAND 
INPUT ■ 


50-MHz 

REFERENCE 



Fig. 1. Block IV carrier tracking receiver 


50-MHz 



M = 20 S-BAND 
85 X-BAND 


Fig. 2. Proposed local oscillator with Dana Synthesizer as loop-controlled oscillator 



INPUT 


Fig. 3. Measurenoent of synthesizer phase noise variance 
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A Prototype DSN X-S Band Feed: DSS 13 Application Status 

(Third Report) 

W. Williams, D. Nixon, H. Reilly, J. Withington, and D. Bathker 
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This article, the third in this series discussing a new prototype X-S band horn feed for 
future use at various DSN sites, deals with the testing of the final fabricated feed at DSS 
13. Measured feedhorn patterns are presented and efficiencies calculated. Preliminary 
results of system noise temperature and 26-m antenna system gain measurements are 
presented. Also discussed are some measurements leading to an improved second 
generation feed. 

The preliminary results of the field measurements indicate that this horn will perform 
as originally specified and required. Also the tests for the second generation feed have 
indicated the potential cause of minor X-band moding. 


I. Introduction 

The first report (Ref. 1) of this series discusses the design of 
a dual-band (X-S) corrugated horn for DSN applications. The 
radiation patterns and efficiencies of a half-scale model were 
presented in Ref. 1 . The second report (Ref. 2) discusses the 
design of a diplexer or X-S combiner for injecting signals into 
the horn. A combiner for the X-S horn was fabricated. Also, 
efficiency calculations were made for the half-scale model 
horn patterns scattered from the Venus (DSS 13) site 
subreflector. 

The final full-scale horn fabrication has been completed, 
and full-scale measurements were made. These measurements 
indicate some additional problems in hybrid moding, which 
must be resolved. The full-scale horn was assembled in a 
feedcone with other electronic equipment and shipped to the 


DSS 13 site, where preliminary temperature measurements 
were made at X-band frequencies. The cone was mounted on 
the 26-m antenna for final evaluation. 

The final full-scale horn pattern measurements were used in 
the JPL theoretical scattering program (Ref. 3) to determine 
final efficiencies after reflection from the DSS 13 site 
hyperboloid with vertex plate. The program (spherical wave 
expansion) takes into account the fact that the subreflector is 
at a range somewhat less than D^/X. 

II. The Final X-S Horn 

The final Phase I X-S horn is shown in Fig. 1 , and measured 
radiation patterns of this horn are shown in Figs. 2, 3, and 4. 
Although these patterns are acceptable, there are some 
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differences from the measured half-scale model at frequencies 
in the high end of the X-band region (Fig. 2). There is a 
marked difference in the E and H plane patterns, the E plane 
being broader with a 1 -dB null (or dip) on axis. This results in 
more forward spillover from this broader E-plane pattern, and 
a somewhat reduced efficiency. 

The E and H plane patterns (Fig. 3) at 7.175 GHz overlay 
much better, although there is evidence of some difference 
near boresight where slight dips occur. Figure 4 indicates a 
shght difference in the E and H plane patterns at 2.295 GHz, 
but so minor as to be ignored. 

It is next necessary to predict the final RF efficiency of 
DSS 13 with this feedhorn in place. This is done by 
decomposing the patterns into Fourier-type components using 
the spherical wave expansion technique of Ref. 3. These 
components are then used to determine currents upon the DSS 
13 subreflector at the proper range from the horn and hence 
determine the final scattered pattern that illuminates the 26-m 
paraboloid. Fifty-eight of these expansion coefficients were 
determined for each X-band horn pattern, and 45 coefficients 
determined for the S-band horn pattern. The resulting scatter 
patterns are shown in Figs. 5, 6 , and 7. The result of an 
unequal E and H horn pattern is noted in Fig. 5 since the 
scattered E and H plane patterns are also unequal. 

Also it is apparent that the DSS 13 vertex plate was not 
designed for X-band. Geometric shadov^ng by the subreflector 
extends to about 8 deg. From Figs. 5 and 6 it is noted that this 
vertex plate achieves its purpose to about 13 to 15 deg, 
wasting a significant part of the inner paraboloidal region. 

Next, the efficiency of these patterns is calculated as 
illuminators of the 26-m paraboloid. Values are tabulated in 
Table 1 . Observe that these totals do not include spar blockage 
or surface tolerance. These are not a function of the X-S feed 
unit. 

It is interesting to understand qualitatively the reasons for 
the slight differences in efficiencies at the various frequencies. 
The horn patterns (Figs. 2, 3, and 4) show the angle of 32.7 
degrees, which is subtended by the total subreflector, thus 
representing the total portion of each pattern that is inter- 
cepted. The first item noted is the greater spillover at S-band 
than at X-band. The reason is evident when looking at these 
radiation patterns: the S-band pattern is triangular in shape, 
reaching the -40 dB level at nearly 36 deg, while X-band 
patterns have much steeper sides (slopes) and reach -40 dB at 
a smaller angle. Thus a greater fractional energy is intercepted 
at X-band. Next is noted the differences in illumination 
efficiency. The superiority of S-band can be noted in Figs. 5 , 
6 , and 7. In this case, the large X-band hole due to the vertex 


plate has drastically reduced the effectiveness of the reflector 
area compared to S-band. The superiority of 8.425 GHz 
illumination efficiency over 7.175 GHz can be seen in Figs. 5 
and 6 . The taper at 7.175 GHz (Fig. 6 ) can be seen as greater 
than 8.425 GHz (Fig. 5), the result of the greater taper on the 
hyperbola intercept (-16 dB, Fig. 3) compared to -14 dB at 
8.425 GHz (Fig. 2 ). Also, phase efficiency is superior at 
S-band, again the effect of the vertex plate having the more 
serious effect at the shorter wavelengths. It becomes apparent 
that the subreflector could be redesigned for this antenna, if 
one desired further X-band improvement, to increase illumina- 
tion and phase efficiency by, perhaps, a total of 5 percent. 


III. The Feedcone 

The new X-S feedhorn and other necessary equipment 
(power lines, brackets to hold masers and other electronic 
gear, etc.) were installed in a feedcone at JPL. The completed 
unit is shown in Fig. 8 . This assembly was sent to DSS 13 for 
further measurement and installation on the 26-m antenna. 

Figure 9 is a photograph of this new installation at DSS 1 3 
and Fig. 10 depicts the older reflex feed as applied to the 34-m 
modifications. 

The plan to evaluate the common aperture X-S horn, as 
installed in a feedcone system at DSS 13, was to measure total 
operating system noise temperature in a normal 

feedcone-on-ground test mode followed by (feedcone 
antenna mounted) and 26-m antenna area efficiency ( 17 ) 
measurements. Two fundamental problems have prevented the 
gathering of reliable data. 

First, the cone-on-ground S-band measurements are clouded 
by a problem discovered after the feedcone was antenna 
mounted. Although JPL Mesa Antenna Range measurements 
at S-band were done with a proper assembly of the feedhorn 
yoke assembly (power division for phasing the four S-band 
horn inputs), an inadvertent reassembly within the feedcone 
introduced a phasing problem. Resultant 26-m antenna main 
beam initial measurements showed a highly squinted and 
elliptical beam prior to proper assembly of the yoke. Thus full 
confidence was not obtained in the present S-band (23 K) 
since the expected 3 K feedcone OFF/ON antenna value (due 
to spillover and scatter) has not been demonstrated. The 
approximate S-band area efficiency under these uncertain 
conditions is 62 percent. 

Secondly, initial feedcone-on-ground X-band measure- 
ments appeared high (21 K) and not stable. Later, with the 
feedcone antenna mounted, the measurements were very high 
(about 30 to 40 K) and very unstable. Additionally, different 
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T^p measurement schemes produced different results. The 
problem was traced to the X-band TWM, an old unit that 
suffered from years of thermal cycling. The travelling wave 
comb structure had become distorted and short circuited in 
places. Now repaired and reinstalled in the antenna-mounted 
feedcone, different measurement schemes all indicate T^p in 
the low 20 K class. During the period with a suspect maser, 
X-band antenna efficiency appeared to be 48 percent. This 
value now requires another evaluation. 

In summary, despite the above problems, the X-S common 
aperture horn (initial receive-only model) appears to have no 
major problems; a full evaluation is planned just prior to (and 
after) the next opportunity for feedcone removal from the 
26-m antenna. At that time, feedcone-on-ground T^p measure- 
ments will be carefully compared with feedcone-on-antenna 
T^p measurements to discern the expected spillover and 
scatter contributions. 

IV. Further Horn Study, Second Generation 

The X-band horn has exhibited evidence of an operational 
problem as noted in the unequal E and H plane patterns of 
Fig. 2. There appears to be some small amount of moding, i.e., 
the generation of other corrugated hybrid modes than the 
desired HEj j mode. For instance, the EHj j mode may be 
generated at low X-band when input grooves are less than X/4, 
or higher modes, the HEj 2 and EHj 2 , when the horn is 
sufficiently large and horn interruptions (discontinuities) may 
occur. It was not obvious what caused these problems. 

Toward the end of solving this problem, several tests were 
made using segments of the full-scale horn that were left over 
from the first model development. An X-band input section 
was available, to about 12.7-cm (5-inch) aperture, and an 
additional section was built for enlarging the sample and to 
allow various tests to be made. Also, an additional input piece 
was made for the X-band, which consisted of a gradual radius 
change from the X-band cylindrical guide to the 17.1 -deg horn 


flare instead of the abrupt 17.1-deg horn attached to cylindri- 
cal waveguide. 

Basically, the following tests were made: 

(1) Metallic taping off two grooves to simulate the com- 
biner input that results in the pitch change for the 
corrugations. 

(2) Building the additional section at a different flare angle 
to simulate the change of flare angle for this horn. 

(3) Using the gradual flare change at the X-band input. 

The following results were noted: 

(1) A significant (very noticeable) change occurs when the 
gradual taper is used. The abrupt change from cylindri- 
cal to the 17.1-deg flare is a major cause of moding. 

(2) Neither the small flare angle change nor the corrugation 
interruption caused a significant moding, although 
some was noticed. 

The conclusion was that a new horn will be constructed using 
the gradual change to the 17.1-deg flare, over about a 10.2-cm 
(4-inch) transition length. Further investigation into the 
optimum groove depth and corrugation pitch is being carried 
on. 

Also in this second-generation study, a new diplexer/ 
combiner must be developed so that the S-band system will 
have sufficient bandwidth for transmission as well as low noise 
reception. Although the first generation diplexer (now at DSS 
13) offered sufficient X-band isolation (Ref. 2), it had an 
operating bandwidth of less than 100 MHz in S-band. The 
second-generation diplexer being designed will use the same 
approach as the first-generation design, except that the radial 
transmission line for S-band injection will be about 70 percent 
wider to obtain the bandwidth, and the X-band radial line 
chokes will be increased from 2 to 6 in number to maintain 
the high X-band isolation. 
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Table 1. Total efficiency as a function of frequency 


Efficiency Type 


Frequency 


8425 GHz 

7175 GHz 

2295 GHz 

Total spillover 

0.96660 

0.97139 

0.88693 

Illumination 

0.86173 

0.82389 

0.88610 

Cross-polarization 

0.98861 

0.99430 

0.99590 

Phase 

0.89627 

0.89301 

0.92438 

Blockage (subreflector 

1.0000 

1.0000 

0.96900 

Totals 

0.73804 

0.70687 

0.70107 



Fig. 1. Final first-generation X-S horn 












Fig. 7. The final 2.295 measured patterns, 45 nruxSes, scattered from DSS 13 subreflector 
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Fig. 9. Th« first-gen«r«tion X-S hom at DSS 13 
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Coding for Optical Channels With Photon-Counting 

R. J. McEliece and L. R. Welch 
Communications Systems Research Section 


We study the problem of coding for Piercers recent model for optical communication. 
We conclude that for any positive rate p (measured in nats per photon), the best code of 
length n has an error probability bounded by an exponentially decaying function of n; 
we exhibit explicit practical schemes for p and give evidence that p ^ 1 may be 

the "^practical limit** for optical communication. 


I. Introduction 

In a recent paper Pierce (Ref. 1) has shown that if one uses 
photon-counting techniques for communication at optical fre- 
quencies, and low noise temperatures, the rate at which one can 
operate is at most hf/kt = 4.80 X 10" ^ ^ ifl'T) nats per photon, 
where / is the photon frequency in hertz and T is the noise 
temperature in kelvin. This number is typically fairly large; 
for example / = 6 X 10^^ Hz, T = 400 K gives a limit of 
72 nats per photon. Pierce observed, however, that if one 
attempts to design a practical communication system for such 
a channel, one encounters severe problems of coding com- 
plexity long before the limiting rate of hf/kT is reached. 
Moreover, these problems are not caused by thermal noise, 
but rather by the nature of the photon-counting process 
itself. We shall now describe Pierce’s “noiseless” photon chan- 
nel, and devote the rest of our article to a discussion of the 
coding problem for this channel. 

We assume that the data to be transmitted consists of a 
long stream of O’s and I’s: x^,x^,x ^, . . . , and that the time 
interval during which these bits are to be transmitted is divided 
into small intervals (“slots”) of duration t^ seconds each. 
The transmitter is a semiconductor laser which is pulsed 
during the ith time interval if and only if Xj = 1 . The expected 
number of photons emitted during such a pulse we denote by 


X. The receiver is a photon counter; it emits a 1 if it is struck 
by one or more photons in a given time slot, and a 0 if it is 
not. If the laser is not pulsed, of course it emits no photons; 
hence a transmitted 0 is always received correctly. On the 
other hand, if the laser is pulsed, because of the Poisson 
statistics which govern photon emissions, there is a nonzero 
probability that no photons will be emitted. Hence a 
transmitted 1 is received incorrectly as a 0 with probability 
e~^. In effect, then, the photon channel is just the “Z-channel” 
familiar to information theorists (see Fig. 1), where the cross- 
over probability e is equal to e"^. 

In Ref. 1 Pierce considered the following coding scheme 
for this channel. There are n distinct binary codewords, each 
of length n. Each codeword has a 1 in one coordinate only. 
Since there are n codewords, each carries log n nats^ of 
information. Furthermore, the transmission of each codeword 
requires exactly one pulse and hence an average of X photons. 
Hence the rate of transmission, measured in nats per photon, 
is given by p = (log n)/\. The probability that the codeword 
will be received incorrectly is just the probability e~^ that 
no photons will be emitted by the laser during one pulse 
period. Since X = p“^ log n, it follows that if we use this 


^Throughout the paper all logarithms will be natural. 
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( 4 ) 


scheme at a rate of p nats per photon^ the error probability 
will be 

( 1 ) 

Hence for any desired rate of transmission p, however large, 
by choosing n large enough, one can make the receiver error 
probability arbitrarily small. 

In practice, however, Pierce’s codes are not of much use, 
unless p is quite small. For example, at p = 1 nat per photon, 
and = 10'^, Eq. ( 1 ) implies that n must be at least 10^. 
The problem is of course that the decay of in n and p 
is very slow. In Ref. 2 it was shown that there exist codes 
of length n for which P^ decays exponentially as a function of 
p"^, p -> 0. In this paper we shall show that for any fixed 
value of p, however large, there exists a sequence of codes of 
increasing length all having rate at least p nats per photon, 
whose error probabilities decay exponentially as a function of 
the code length. In Section II we will prove this fact (the main 
idea of the proof is to use Pierce’s coding scheme as a modula- 
tion scheme for more complex codes); in Section III we will 
exhibit some codes which are practical and give good perfor- 
mance for p < “^1 nat per photon. (These simple codes will 
outperform any possible coding strategy if one uses coherent 
signal amplification rather than photon-counting techniques, 
because p = 1 nat/photon is the capacity of the optical chan- 
nel when linear amplification is used (Ref. 1 ). Finally, in 
Section IV we will give evidence (but no proof) that it is 
probably very difficult to design practical systems that operate 
at rates much larger than 1 nat per photon. 

II. An Exponential Error Probability Bound 

In this section we will show that there exists a function 
£’(p), which is positive for all p > 0 , such that for any p > 0 
there exists a sequence of codes Cp C 2 , C 3 , . . of increasing 
block length, each having rate at least p nats per photon, such 
that the decoding error probability for the ki\v code satisfies 

Pe ^ exp [-n^Eip)] , (2) 

where rij^ is the length of the A:th code. 

Thus let p > 0 be fixed. Choose a positive integer q 
satisfying 

q>eP. (3) 

Having selected q, choose a real number R, 0 < R < \, 
satisfying 


- P)<P ^ log 

(This is possible since by Eq. (3) p"i log p > 1 ; and the left 
side of Eq. (4) is an increasing function of R and approaches 1 
as/^ 0). 

Now consider a p-ary erasure channel with erasure proba- 
bility e. This channel has input alphabet = {1,2, . . ., q} 
and output alphabet U {e}, where “e” is an erasure 

symbol. The channel’s transition probabilities are 

{ 1 ~ c if X 
€ if y = e, 

0 otherwise 

Thus each symbol is transmitted correctly with probability 
1 - e, and “erased” with probability c. This channel has 
capacity (1 ~ e) log p, and so according to the noisy-channel 
coding theorem (see Ref. 3, Chapter 5), for any R <(l - e), 
and any positive integer there exists a code of length with 
at least q^^ codewords, whose decoding error probability 
satisfies 

P^ < exp [-NE^(R )] , (5) 

where the random coding exponent for the channel, 

is a function which is positive for all 7? < 1 - e. 

Our object is now to use these codes, on the “photon 
Z-channel.” To do this we must have a “modulation” scheme 
which will allow us to transmit symbols from^^ = (1,2, ... , 
p } over the Z-channel. The modulation scheme we use is as 
follows: we assign each symbol from a unique binary 
vector of length p containing one “1,” For example, with 
p = 4 the assignment could be 

1 : 1000 
2 : 0100 
3 : 0010 
4 : 0001 

(This is exactly Pierce’s quantized pulse position modulation 
described in Ref. 1 and Section 1 .) 

When one of these length p binary patterns is transmitted 
over the Z-channel, the “demodulator” looks for a “1” among 
the p received bits; if a “1” is detected in the A:th position, the 
demodulator knows for certain that the transmitted symbol 
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was k. If, however, no I’s are detected, the demodulator has 
no idea what was sent and outputs the erasure symbol “e.” 
It follows that when this particular modulation scheme is 
adopted, the Z-channel becomes a ^-ary erasure channel 
whose erasure probability e is equal to the crossover proba- 
bility of the Z-channel. 

The code of length N has at least codewords and so 
each codeword conveys at least RN \ogq nats of information. 
Since transmission of each codeword requires exactly N pulses, 
we are getting at least R \ogq nats per pulse. If we are trans- 
mitting at rate p nats per photon, this means that we need at 
most, 

X = R\ogq photons per pulse. (6) 

Hence if we operate the laser at X photons per pulse, where 
X is given by Eq. (6), our coding scheme will operate at p nats 
per photon. Now according to the Poisson statistics governing 
photon emission, this means that the erasure probability e 
for the ^-ary erasure channel is 

e = q-^l? (7) 

Hence if the code’s (f 7 -ary) rate R is less that the channel’s 
(^-ary) capacity 1 - e, we can get exponential decay of the 
error probability. This condition is 

R < 1 (8) 

which is just a restatement of Eq. (4). Hence it follows that for 
every N there exists a code of rate R, whose error probability 
satisfies Eq. (5). This code, viewed as a code for the Z-channel, 
has length n = qN, and so its error probability satisfies 

< ^xp[~nEJiR)lq], (9) 

Since for fixed p > 0, and q,R satisfying Eqs. (3) and (4) we 
have E^{R)!q > 0, it follows that there exists a function 
E{p) such that 

<exp[-nE(p)], (10) 

which is what we set out to prove. 

We have made no attempt to give numerical bounds on 
E(p); it would be interesting to do this, however. 

In the next section, we will show how the technique just 
described can be used to design practical coding schemes 
for the photon channel in the range p < nats/photon; 


III. Practical Schemes for p < 1 

In this section we will give explicit examples of the coding 
schemes described in Section II; these will turn out to be 
easily implemented and have low error probabilities, provided 
p is not much larger than 1 nat/photon. 

As before, for a given rate p, choose a positive integer 
q > exp(p). Using the modulation scheme described in Sec- 
tion II, the Z-channel becomes a ^ 7 -ary erasure channel, and 
we must design a code over the alphabet capable of cor- 
recting as many erasures as possible. An obvious choice is the 
class of Reed-Solomon codes, which require only that q be 
a power of a prime. 

Over the alphabet (which we now view as a finite field), 
Reed-Solomon codes are linear cyclic codes of length 1, 

arbitrary dimension \ < k < n - \, with minimum distance 
d = n - k 1. They will correct any pattern of up to n - 
erasures, and very efficient encoding and decoding procedures 
are known, when q is sl power of 2 (see Ref. 5, Chapter 8 for 
details). 

If we use an (n,k) Reed-Solomon code for the present 
application, each of the q^ codewords carries k log q nats of 
information, and each codeword requires n pulses. Thus if we 
are transmitting at p nats/photon, the average number of 
photons per pulse is at most (cf. Eq. ( 6 )). 

X = — • p“Mogp photons per pulse. (11) 

It follows that the erasure probability for the corresponding 
^ 7 -ary erasure channel is (at worst) 

e = = q-^!f>^ ( 11 a) 

vvhere R = kjn is the code rate. Since the RS code can correct 
all patterns of up to « - erasures, it follows that the decod- 
ing error probability P^ satisfies 

E (iy (>-<>"-'■ (12) 

/=n-k+l ' 

In Fig. 2 we have plotted P^ vs p for four typical RS 
codes. The curve labelled <7 = 16 is a (15,8) RS code with 
q = 16; the others are (31,16) q = 32; (63,32), q = 64; and 
(127,64), ^7 = 128. It is seen that the performance of all of 

these codes degrades rapidly as p increases beyond 1. (It is 

to be recalled that as codes for the photon channel, the length 
is actually 16 * 15 = 240 for the ^ = 16 codes; « = 31 • 32 = 
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992 for the q = 32 code; n = 4032 for q = 64; and n = 16256 
for q = 128.) Each of these codes is the best of its length, at 
least in the limit as p 0, and so no significant improvement 
could be expected merely by altering the code’s rate. We 
conclude that for p < '^l practical schemes exist. In the 
next section we give evidence that p > ^1 may be very hard 
to achieve. 


Thus the erasure probability for the Z-channel is 

6 = exp(-/?/p5) (16) 

and the inequality Eq. (13) becomes 


IV. /?o for the Photon-Detection Channel 

Associated with any discrete memoryless channel there is 
a number C, called channel capacity, which represents the 
theoretical maximum rate at which one can communicate 
reliably over the channel. However, it is usually extremely 
difficult to operate at rates very close to C, and it has been 
conjectured (see e.g,, Ref. 4) that another number, which 
is associated with every discrete memoryless channel, is a more 
realistic measure of the maximum rate of practical reliable 
communication. In this section we shall show that for our 
photon-counting Z-channel, the value of p corresponding to 
Rq is exactly 1 nat/photon; this number, incidentally, is the 
limiting value (capacity) for coherent optical communication 
systems (Ref. 1). 

Given a discrete memoryless channel with input alphabet 
output alphabet Ay, and transition probabilities p(ybc), 
we can state the “Rq - conjecture” as follows. “Any practical 
communication system, in which the relative frequencies with 
which the input symbols x E A^ are used are described by 
the probability distribution p{x), must operate at a rate R 
satisfying 


R < -log 1 1 - 26(1 - 6) [1 - exp(-i?/2p5)][ (17) 


Our result is the following. 


Theorem: There exist real numbers (R,6) 0 < R < log 2, 
0 < 6 < 1 satisfying Eq. (17) if and only if p < 1 . Hence any 
system for the photon-counting Z-channel for which p > \ 
must violate Eq. (13). 


Proof: First assume p < 1, and define a- Rj2b. Then a > 0 
and Eq. (17) becomes 


2a6 < -log 1 1 - 26(1 - 6) [1 - exp(-o;/p)] [ . (18) 


Now choose a satisfying 


a < 1 - (19) 


R < ~-log ^ 

ye A y 




XeA , 


(13) 


In the specific case of the Z-channel with erasure proba- 
bility e, with p(x=0) =1-6, p(x=l) = 6, the quantity on the 
right side of Eq, (13) is given by 

= -log [1- 26(1 -6)(l-e‘/")]. (14) 


(This is possible since p < 1.) For this fixed value of a the 
right side of Eq. (18) is equal to 25(1 - + 0(6^). It 

follows that Eq. (18) is satisfied for all sufficiently small 
5 > 0. Thus with/? =2a6 and 5 small enough, Eq. (17) is 
satisfied. 


Now assume p > 1. Then exp(-a/p) > exp(-a), and the 
right side of Eq. (18) is bounded as follows: 


If we are using a code of length n and rate R for the photon- 
counting Z-channel, in which the average number of I’s per 
codeword is 5«, a calculation similar to that which led to 
Eq. (16) and Eq. (11) shows that the average number of 
photons per pulse must be 

5 = /?/p5 photons per pulse. (15) 


-logjl - 26(1 - 6)[1 - exp(-a/p)] [ 

<-log I 1 - 26(1 - 6)[1 - exp(-a)][ (20) 

If we replace the right side of Eq. (18) by its upper bound 
(Eq. 20)), a little rearrangement yields 
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<28(1-6). 

1 - exp(-a) 

Lemma: Let^>0,/(v) = (1 - j^^)/(l - Then 

(k for A: < 1 
inf /O) = < 

0<>'<i (l for k>\. 


(The proof of the lemma is elementary and is omitted.) 

Now let y = exp(”a), k - 28. Then the left side of Eq. (21) 
is (1 - y^)l(\ - y). By the lemma, this is >28 if 26 < 1 ; thus 
Eq. (21) implies 26 < 25(1 - 26), a contradiction. If 26 > 1, 
then Eq. (21) gives 1 < 26(1 - 5), another contradiction, 
since 25(1 - 5) < 1/2 for 0 < 6 < 1. This completes the proof 
of the theorem. 
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In this article, a combination of the Chinese Remainder Theorem and Winograd's 
algorithm is used to compute transforms of odd length over GF(2^ ). Such transforms are 
used to compute the syndromes needed for decoding BCH and RS codes. The present 
scheme requires substantially fewer multiplications and additions than the conventional 
method of computing the syndromes directly. 


I. Introduction 

It is shown in Ref, 1 that the finite field transform can be 
used to compute the syndromes for BCH and RS codes. The 
disadvantage of the transform method over GF(2'”) is that the 
transform length is an odd number, so that the most efficient 
FFT algorithm cannot be used. In this article a combination of 
the Chinese Remainder Theorem and Winograd’s algorithm is 
used to develop a fast finite field transform over GF(2^). 
Such a fast transform can be used to compute the syndromes 
needed when decoding BCH and RS codes (Ref. 2). It is also 
shown by example that the number of multiplications and 
additions of this new scheme for computing the syndromes is 
substantially fewer than used in the direct method. 

II. The Computation of the Syndromes for 
BCH and RS Codes 

Let n be the block length of a BCH or RS code of designed 
distance d in GF{2^). Also denote the received vector by r = 
(r^, Tj, . . . , If 0^^ 1 </< d - 1 are roots of the code’s 


generator polynominal, then the decoder needs to evaluate the 
syndromes Sj for 1 </ < d - 1 from the received vector by 

n-l 

S. = ^ r o'/ for 1 </<d- 1 (1) 

7 = 0 


Evidently, a direct computation of Sj for 1 < / < c? - 1 
involves (d - 1)(« - 1) multiplications and {d - l)(n - 1) 
additions. Using the Chinese Remainder Theorem, it is shown 
in this section that this number of multiplications and 
additions can be reduced substantially. To see this, if n=n^ 
«2 . . where (n.,ni)= 1 for i¥=j\ then by the Chinese 
Remainder Theorem it follows (Refs. 3,4) that the A7-point 
transform over GF{2^) defined in (1) can be decomposed 
into a multidimensional transform over GF(2"*) as follows: 

Given an integer /, define f(j) = 0 j » /2 » • * • » G)> where = / 
mod n^ (1 Then • • • » o,i,o, . . . , o)^ 

where 1 is in the position, is a primitive root of 
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unity. The computation of given by (1) can now be written 
as 


Assume the number of multiplications and additions needed to 
perform the -point transform for 1 < / < /• - 1 is and a,., 
respectively. Then by Eq. (3) it is evident that the number of 
multiplications and additions needed to compute Stage 1 up to 
Stage (r - 1) is given by 




i i 1 £ r 


'l=0 «2 = 0 


"l “ “2 " • • • 




Vr 


(2) 


for 1 < / < d - 1 


From Eq. (2) this algorithm consists of the following r stages: 


Stage 1: 


n-l 


> = V* ^ Vr 

Z 'r)“ 


'■ =0 


for 0 < / < n - 1 

J ^ K 


Stage 2: 


Si .. 


Vi-‘ 








V_.=o 




r-Vr- 


m. r-\ Q 

M = n 2^ A = n ^ , respectively. 

i=l ' /=! 

To compute the syndromes one needs only to compute the 
first d points of the transform over GF{2^) defined in Eq. (1). 
Thus, at the Stage r in Eq. (3), one needs only to compute = 
, /V) 1 Hence the number of both 

multiplications and additions needed to compute Stage r in (3) 
is [d- l)(rtj - 1). Combining this with the above results the 
total number of multiplications and additions needed to 
compute (1) is finally M {d - l)(«j - 1) and A {d - 1) 
(rtj - 1). respectively. 

Suppose that the transform length given in Eq. (3) of the 
first r - 1 stages is small. Then it was shown (Ref. 1) that such 
short length transforms can often be computed by a modifica- 
tion of Winograd’s algorithm. Algorithms for computing a 
transform over GF(2"*) of n points for n = 3,5,7,9.17 are given 
in detail in Ref. (1). The detailed algorithms for computing the 
/2-point transform over GF(2'”) for n - 3,5 are given in the 
appendix. These transforms are used in the following example 
to compute a transform of 2®- 1 points. 

Example 1: Let n = 255 be the block length of an RS code 
of designed distance d = 33. This code will correct t errors 
where 2/^ <33. The first step of the decoding process is to 
compute the 32 syndromes as follows: 


‘‘orO<Vj<Vi- 1 


255-1 

S. = r.o^' for 1 1 = 32 (4) 

1=0 


: where a is an element of order 255 in GF(2®). If Eq. (4) is 

computed directly, the number of multiplications and addi- 
tions needed is 32(255 - 1) = 8128 and 32(255 - 1) = 8128, 
respectively. Since n = 255 =22^/22/23 = 17 • 5 • 3, (4) for this 
Stage r: case reduces to 


^1-1 


r 

* Ij=0 * 

forl</<d-l ( 3 ) 


Stage 1: 




'■ 3 =® 
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In Eq. (5) Stage 1 is a 3-point transform and Stage 2 is a 
5-point transform. These two transforms are computed by a 


modification of Winograd’s algorithm. It is shown in the 
appendix that the number of multiplications and additions 
needed to compute the 3-point transform is 1 and 5, 
respectively. Thus, the number of multiplications and addi- 
tions needed to compute Stage 1 is (17)(5)(1) = 85, and 
(17)(5)(5) = 425, respectively. Again by the appendix one 
observes that the number of multiplications and additions 
needed to compute a 5-point transform is 5 and 17, 
respectively. Thus, the number of multiplications and addi- 
tions needed to compute Stage 2 is (17)(3)(5) = 255 and 
(17)(3)(17) = 867, respectively. Since one needs to compute 5y 
only for /= 1,2, .. . ,32, the number of multiplications and 
additions needed to compute Stage 3 is 32(17 - 1) = 512 and 
32(17 - 1) = 512, respectively. Thus the total number of 
multiplications and additions needed to compute the syn- 
dromes Sj for l</<32 is 85 +255 + 512= 852, and 
425 + 867 + 512 = 1 804, respectively. 
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Appendix 

A Summary of Transform Algorithm of n-Points for n = 3, 5 


For n-3, let 7 be a primitive cube root of unity in 
GF{2^), where 31(2"^ - 1). The 3-point transform over 
GF(2'”)is 


3-1 


^ 7*^ for0<^<2. 


1=0 


Algorithm for n = 3: 




=5^ • 7 , ^2 = 1 • 




2’ 


^ 3 =^ 0 +"^! ^^ 4 =^ 3 +"' 2 ’^ 5 =" 3 +"' 3 ^ 
v4q “*^q >^2 *'^2 ”'^5 ^ ^ 


Hence, from Eq. (A-1), the total number of multiplications 
and additions needed to perform a 3-point transform is 1 and 
5, respectively. 

Next consider the case n-S. Let 7 be a primitive 5^^ root 
of unity in GF(2'”), where 5 1(2'” - 1). The 5-point transform 
is 


5-1 

" 'll forO<A:<4 

1=0 


Algorithm for n = 5: 

Sj =(a2+fl3),S2=«i+V'3=''i +^^3 




Wq = 1 • ^6, = ( 7 ® + 7 ^) - s ^, m ^= ( y ^ + 7'*) • ; 


OT3 = (7* + 7^) • s^, = (7 + 7'*) • s^,m^=(y + 7“^) • s^; 


m^ = 1 ' a^,m^= I ’aym^=l • a^.w^ = 1 • a^; 


s,=mo + Wi,Sg=S7+m2,S9=S7+m3,Sjo = W5+m5; 


Sjj =^4 = mj +mg,5,3 =m4 +m7,s,4 =Sg +jj^; 


■^15 '*8 '*'^1 * 8 '*’* 12’*17 *9 “'" o ’ 


*14’"^2 *15’"^3 *I6’^4 *I7‘ 


(A-2) 


Hence, from Eq. (A-2), the total number of multiplications 
and additions needed to perform a 7-point transform is 5 and 
17, respectively. 
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Ne\^ X-band antenna feed assemblies with dual-polarization capability are being 
implemented in the DSN 64-meter stations. Together with dual X-band travelling wave 
masers, they permit the simultaneous reception of right- and left-hand circular 
polarization from the Voyager spacecraft. The new feed also includes a dual hybrid mode 
feedhom which increases the antenna gain by 0.36 dB over the present feedhom. 


I. Introduction 

The X-band feed cone assembly (Ref. 1) on the DSN 
64-meter antennas provides the capability for selectable 
right-hand (RCP) or left-hand (LCP) circular polarization 
(Ref. 2). Originally designed to work with a single travelling 
wave maser (TWM), the XRO cone assembly has recently been 
upgraded to include dual TWM’s (Ref. 3). As part of that 
reconfiguration, provision was made for a new feed assembly, 
the Mod III XRO, that would include two basic improvements 
over the feeds that have been in use since early 1977. The first 
improvement is the development of an orthogonal-mode 
transducer that permits simultaneous reception of two differ- 
ent polarizations. The second improvement is the incorpora- 
tion of a dual hybrid mode feedhom (Ref. 4) to increase the 
antenna gain. 

II. The Mod III XRO Feed Assembly 

The new feed assembly is shown in Figs. 1 and 2. It is 
designed to mount within the XRO cone assembly on the same 


circular mounting plate as the previous (Mod I and Mod II) 
feeds, so that replacement in the field is expedited. Each 
output arm of the orthogonal mode transducer has its own 
waveguide switch for TWM calibration. One of the precision 
ambient terminations may be seen at the bottom of the feed 
assembly. 

The new orthogonal-mode transducer is shown in Figs. 3 
and 4. The unit is based on a design by Microwave Develop- 
ment Laboratories, developed under contract to JPL, with the 
final units fabricated by the electroforming process. The unit 
has an isolation of 35 dB or better over the frequency band of 
8400 to 8500 MHz, with the return losses shown in Table 1. 


III. Noise Temperature Measurements 

Using a vertically oriented single hybrid mode (corrugated) 
horn on top of a TWM for reference, the noise temperature 
contribution of the Mod III XRO feed was determined in 
ground tests. By disconnecting the reference horn at the TWM 
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input and connecting the Mod III XRO feed at the same point, 
variables such as weather, smog, TWM noise temperature, 
receiver followup noise temperature, etc., are compared in a 
relative measurement and the loss contribution of the added 
feed components can be determined with high precision (est. 
±0.15 K). From the throat of the horn to the TWM interface, 
the following components contribute to the additional noise 
temperature compared to the simple horn: rotary joints (2 ea), 
quarter-wave plate polarizer, orthogonal mode transducer, 
waveguide twist (in the side path measured), and waveguide 
switch. The measurements were performed on two different 
feeds (for DSS 43 and DSS 63) and the resulting 2.43 K 
increase (over the simple horn) was the same for each feed 
within experimental tolerances. This indicates that the 64- 
meter antenna system temperature at zenith should not exceed 
25 K, including the dichroic plate losses and the normal 
antenna and receiver noise contributions. 


IV. DSN Implementation 

The first Mod III XRO feed assembly has been installed at 
DSS 14, and measurements (Ref. 5) have confirmed 0.36-dB 
antenna gain increase. For reasons yet to be determined, the 
system noise temperature is indicating several kelvins too high. 
Since the system noise temperature was also too high before 
the feed change, the problem is thought to be elsewhere, and 
efforts are underway to locate and correct the problem. 
Nevertheless, excellent data was received from Voyager 1 as it 
passed through the Jovian system, verifying the feed 
capability. 

The second and third Mod III feeds have been fabricated 
and tested and, together with some additional parts for the 
cone assembly, will be implemented at the overseas sites in late 
1979 and early 1980. 
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Table 1. Impedance match characteristics of X-band orthogonal 
mode transducer 


Frequency, MHz 

Return loss, dB 

Thru port 

Side port 

8400 

-24.0 

-28.0 

8425 

-25.5 

-33.5 

8450 

-28.5 

<-40. 

8475 

-29.0 

-32.0 

8500 

-26.0 

-27.0 



Fig. 1. XRO Mod III feed assembly 
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Fig. 2. Close-up of feed assembly showing orthogonal mode 
transducer and switch 


* # 

# • 

# 

# 



Fig. 3. Orthogonal mode transducer showing output flanges 



Fig. 4. Orthogonal mode transducer showing input flange 
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An introduction to the Multi-Mission Receiver (MMR) is presented. The MMR contains 
occulta tion receiver channels (OCC) and very long baseline interferometry (VLBI) chan- 
nels. The design considerations, performance and implementation of the OCC channels 
are discussed. The VLBI channels will be discussed in a subsequent report. 


I. Introduction 

Measurement of the signature of radio signals for obtaining 
radio science, radio astronomy and tracking data requires 
extremely stable receiver channels. The Multi-Mission Receiver 
(MMR) is being developed to provide this capability. 

There are two different types of signal signatures that are 
involved in radio science and in radio astronomy and tracking. 

(1) Occultation channels. The first type are narrowband 
spacecraft transmitted carrier signals. Occultation 
(OCC) channels used to receive these signals measure 
the effect of planetary and solar media, interplanetary 
media, and gravitational fields by detecting minor 
perturbations on these received signals. Real-time 
reduction techniques are used in these channels to 
minimize the data processing costs. This is accom- 
plished by programming the receiver local oscillator 
based on signal frequency prediction. 

(2) Very long baseline interferometry channels. The 
second type are wide-bandwidth signals received on 
the very long baseline interferometry (VLBI) channels. 


To support these measurements, the MMR contains both 
OCC and VLBI channels (Fig. 1). Functions such as control, 
instrumentation and power are common to both sets of 
channels. The functional requirements and performance, 
together with a detailed description of the OCC channels, 
are discussed in this report. The VLBI channels will be dis- 
cussed in a later report. 

II, MMR Occultation Channels 

A. General Description 

Four occultation channels are being implemented (Fig. 2), 
right and left hand circularly polarized channels at both S- and 
X-band. The right-hand channels are equipped with very nar- 
row, narrow and medium bands, while the left-hand channels 
are equipped with the medium band only. The very narrow- 
band channels are to be used for measurement of solar media 
and the narrow-band channels for measurement of planetary 
media. The medium bands of both the right- and left-hand 
channels are to be used for measurements during the Saturn 
ring experiment. In the future, a selection of narrow and very 
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narrow bandwidths might also be included in the left-hand 
channel. The selection of bandwidths available (Table 1) 
provides for all presently known requirements. 

B. Design Considerations 

The most significant features of the occultation channels 
are (1) a wide selection of operating bandwidths, (2) the use of 
one common oscillator to derive the first local oscillator of all 
channels, and (3) extreme phase stability. These features are 
discussed in more detail as follows: 


lator source is a Dana Synthesizer. With the control assembly, 
the Dana Synthesizer frequency can be programmed to com- 
pensate for doppler and maintain the signal within the receiver 
narrow bandwidth. At present, the controller interfaces 
directly with the computer in the Radio Science Subsystem. 
At a latter date, when the MMR controller is implemented, 
the synthesizer control will interface through the Star Switch 
Controller to the Radio Science Subsystem. 

The synthesizer frequency (Fig. 5) required to maintain the 
signal within the passband for a narrow or medium bandwidth 
S-band channel is 


1. Receiver bandwidth selection. The selection of the 
operating bandwidth is determined by such characteristics 
of the received signal as frequency uncertainty and spectral 
broadening. If a frequency uncertainty exists, the uncer- 
tainty at X-band is 1 1/3 that at S-band. In this situation, 
narrow and medium bandwidth channels are used and these 
filter bandwidths are designed with a S- to X-band ratio of 
approximately 3 to 1 1 . If the signal being tracked has essen- 
tially no frequency uncertainty, then the filter bandwidth 
for both S- and X-bands can be identical, the bandwidth 
required being determined by spectral broadening. This is 
the case for the very narrow bandwidth channels. 


BW„ 

fs-fLO^ = 300 MHz 


where 


= 5-band downlink frequency 
= 5-band local oscillator frequency 




BW. 

s 300 — 


48/, 


DANA 


a. Narrow and medium bandwidth channels. A typical 
filter bandpass response at the 10-MHz IF amplifier is shown 
in Fig. 3. This response translates directly down to video 
frequencies. The local oscillator synthesizer is programmed 
to maintain the signal as close to the filter bandwidth center 
as the accuracy of the spacecraft’s predicted frequency per- 
mits. The filters have been designed so that overflow into the 
operating bandwidth due to image noise in the last mixer 
or aliasing noise due to digitizing will be at least 20 dB 
down. This rejection will be increased to 23 dB on future 
implementations. 

b. Very narrow bandwidth channels. In contrast to the 
medium and narrow bands, the very narrow filter bandwidths 
are the same for S and X. However, since the X-band signal 
in the IF is offset from the IF center frequency 1 1 /3 as much 
as the S-band, the filter center frequencies for S and X differ 
(Fig. 4a). The passband at video frequencies for both the S- 
and X-bands do not go down to dc (Fig. 4b); however, with 
undersampling techniques aliasing causes the passband to 
appear as though it goes down to dc (Fig. 4c) (Ref. 1). With 
this technique a sampling rate of no greater than 2 BW is 
required. 

2. Common oscillator. The mechanization of the local 
oscillators is similar to the technique used in the Block III 
open loop receiver; however, the first IF is designed for 
300 MHz to accommodate the X-band channels. The oscil- 


BW^ 

2 

f 

■’DANA 


= y (5-band sampling bandwidth) 
= synthesizer frequency 
BW„ 

/y - 300 MHz — ^ 

48 


Likewise for the X-band channel: 


BW^ 


^LO 


= 300 MHz + — 

X 2 


fyr - 2T-band downlink frequency 
= X-band local oscillator frequency 
BW^ 


= 11(16/ 


DANA 11 


+ ^ MHz) 


BW 


X 1 

^ — = y (X-band sampling bandwidth) 


/, 


DANA 


BW. 

/y- 1 100 MHz 
11 X 16 
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or since = 11/3 = \\j3BW^ 


B. Operation Controls 

The following operation controls are required: (1) The 
s adjustment of the output video signal level (which is a func- 

- 300 2 tion of the received signal strength), and (2) the programming 

^DANA “ 48 of the local oscillator frequency to track the incoming signal. 

1 , Output level adjustment. The video output channels 

of the receiver are digitized on a converter with an input volt- 
age range of ±5 V. Under strong signal-to-noise ratios the 
signal output level is set to 5 (0.707) or 3.54 

(+24 dBm into 50 ohms). Likewise, with a low signal-to- 
noise ratio, the output is essentially noise and, allowing for 
5 X sigma noise peaks, the output level is set at 5/5 or 1 
(+13 dBm into 50 ohms). For any ratio of signal-to-noise 
between these two values, the output level is set in accor- 
dance with the graph of Fig. 6. The output level adjustment 
is made by locally setting manually controlled IF attenuators. 

2. Local oscillator frequency programming. The other 
control required to operate the OCC channels is the local oscil- 
lator frequency programming. The computer in the Radio 
Science Subsystem is used to program the frequency synthe- 
sizer of the local oscillator. No operator intervention at the 
receiver is required. 


A. Configuration Controls 

To set the receiver in the required configuration, it is neces- 
sary only to select the desired operating filter bandwidths. 
The filter select switch is a manual control and must be 
selected locally. This mode of operation will be required 
until the MMR controller is implemented, which will occur 
when the VLBI channels are added and the MMR implemen- 
tation is completed. 


Reference 

1 . Gregg, W. David, Analog and Digital Communication, John Wiley and Sons, New York, 
pp. 286-294. 


IV. Implementation 

The occultation channels are being implemented in phases 
as required to support mission tracking requirements. Table 2 
shows a summary of the implementation schedule. 


Equation (2) is identical to (1), demonstrating that the com- 
mon Dana Synthesizer oscillator can handle both S- and 
X-band channels. This can also be shown for the very narrow 
bandwidth channel. 

3. Improved phase stability. The most significant param- 
eter of the received signal is the information contained in its 
phase. To improve the accuracy of the phase data, the receiver 
phase instability is reduced to a minimum by enclosing the 
most phase-sensitive portions of the receiver (local oscillator 
multipliers) in a temperature-controlled assembly. Phase stabil- 
ity requirements are listed in Table 1 . 

III. Operations 

Two functional types of control are required, configuration 
and operation. 
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Table 1. Occultation channels 


Sampling bandwidth 
(BW), kHz 

S-band 


X-band 

RCP 

LCP 

RCP LCP 

Very narrow BW 

0.1 

_ 

0.1 


0.5 

- 

0.5 


1.0 

- 

1.0 

Narrow BW 

0.82 

_ 

3.0 


2.04 

- 

7.5 


4.09 

- 

15.0 


8.18 

- 

30.0 

Medium BW 

50 

50 

150 150 

Phase instability, 1000 s 




Single channel 

S-band 


< 36 deg rms 


X-band 


<132 deg rms 

Differential 

S-3/11 X 


<10 deg rms 


s-band RCP-LCP 

< 1 


X-band RCP-LCP 

< 1 


Table 2. Implementation schedule 


Operational date S - and X-bands 


DSS 

Project 


BW - RCP channel 



BW - LCP channel 




Medium 

Narrow 

Very narrow 

Medium 

Narrow 

Very narrow 

63 

Voyager 1 Jupiter encounter 
Pioneer 2 Solar conjunction 
Voyager 1 Saturn encounter 
Galileo Faraday rotation exp 

6-80 

10-78 

7-79 

6-80 

7-83 

7-83 

14 

Voyager 2 Jupiter encounter 
Pioneer 2 Solar conjunction 
Galileo Faraday rotation exp 

7-83 

5-79 

7-79 

7-83 

7-83 

7-83 

43 

Voyager 2 Saturn encounter 
Galileo Faraday rotation exp 

7-83 

5-81 

5-81 

7-83 

7-83 

7-83 
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A arrow indicates direction of output frequency 

CHANGE WITH INCREASING INPUT FREQUENCY 


Fig. 3. Medium and narrow band filter responses 





(c) APPARENT (ALIASED) S AND X-BAND VIDEO SPECTRUM 


A ARROW INDICATES DIRECTION OF OUTPUT FREQUENCY 
CHANGE WITH INCREASING INPUT FREQUENCY 


Fig. 4. Very narrow band filter channel spectrum: (a) lOCMcHz IF 
spectrum, (b) video spectrum, (c) apparent (aliased) S- and X-band 
video spectrums 
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Fig. 5. Occultatlon channels - first local oscillator 



Rg. 6. OCC channel operating level 
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Preliminary Report on DSN System Performance Under 

Local Weather Effects 

J. M. Urech 
Deep Space Station 62 


Local effects of precipitation are studied, and a simplified working model is developed. 
Experimental results obtained by simulation are in good agreement with the model, show- 
ing that it could be an important contribution to system degradation. If this is confirmed, 
some suggestions for improvement are presented. Nevertheless, definite results with actual 
rain are not yet available. 


I. Introduction 

The effects of the weather on X-band telecommunications 
system performance have generally been of major concern. At 
stations having only S-band capability, however, subject 
effects have proven to be very minor for normal tracking. 
Nevertheless, during DSS 62 continuous support of the plane- 
tary radio astronomy program (Jupiter Patrol), the weather 
problem became a major concern as the minor SNT fluctua- 
tions presented catastrophic results when measuring Jupiter 
temperature of approximately 0.5 K. In some cases the obser- 
vations had to be cancelled because of rain. 

As a result of this problem, it was decided to consult the 
DSN/Flight Project Interface Design Handbook (810-5) and 
Ref. 2. In both cases, statistical and probabilistic data are pre- 
sented for the effects (mainly at X-band) of nonprecipitating 
clouds. Although this information is especially important for 
long-term planning, the interest of the station was mainly on 
the total effects, including precipitation. Although the probab- 
ility of precipitation is small, some radio astronomy passes 


at S-band have undoubtedly been spoiled by rain and critical 
X-band reception has been badly hampered. 

In a further search for weather performance data, other 
available documents were consulted (Refs. 2-4). However, the 
information obtained presented a set of curves giving the path 
losses (dB/km) due to rain for different frequencies and rain- 
fall rates. By using these path losses (which do not include 
direct antenna effects), the equivalent system noise tempera- 
ture (SNT) increase may be estimated for certain rain assump- 
tions. The estimates for equivalent SNT increases are normally 
lower than the values encountered in actual rainfalls. 

The extra SNT during rain (not totally explained by cloud 
attenuation and path losses) as well as the fact that zenith SNT 
is exceptionally high once rain has subsided has led to the 
belief that there exists another important contribution related 
directly to rain falling on the antenna (Mylar and/or reflector). 
As this potential effect was not contemplated in the consulted 
documents (Refs, 2-4), a simple experiment was performed to 
clarify our assumptions. 
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With the DSS 62 antenna at zenith, a plastic tray was set in 
front of the horn in order to measure the S-band SNT as a 
function of the thickness of the various water layers. The 
initial result of the experiment was absolutely surprising - a 
system noise temperature of 180K was achieved with only 
about 0.5 mm of water. Because of this unexpected result, the 
entire system was completely checked out and the experiment 
repeated. However, the second attempt yielded the same 
results and produced alternating maximum and minimum SNT 
values for increasing water thickness. The results were not 
firmly accepted and doubts existed as to whether they were 
correct and the validity of their interpretation. Nevertheless, 
it was encouraging to think that the station was on the right 
track toward the understanding of the precipitation problem. 


II. Simplified Working Model 

As a first step to clarify the high SNT contribution, studies 
were made on the properties of water as a thermal radiation 
absorber and emitter. This, however, did not help as a water 
layer would only contribute less than 10 K through simple 
absorption and emission. A clear explanation would probably 
require an expert in electromagnetic theory, microwaves and 
radiation theory. In the absence of the latter, numerous docu- 
ments were consulted on the above fields in order to shed 
some light on the problem (Refs. 5-7). Although the effect 
is still not fully understood, it is believed that our simplified 
working model may help explain the increase in SNT. 

A. Effects of Water on Mylar Covering the Horn 

Water, as an imperfect dielectric will partially attenuate the 
electromagnetic radiation by absorption and consequently 
(Kirchoff law of radiation) partially emit thermal radiation. 
The most important fact in our case, however, is that the 
refraction index is fairly high for S- and X-band, thereby 
acting as a good reflector. If the horn is considered together 
with the waveguides, diplexer or polarizer (S- or X-band), 
etc., all the energy received within the design bandwidth will 
be absorbed with almost no reflection. Recriprocally, all 
internally generated energy will be radiated out of the horn 
with almost no back reflection. This may be equivalent to con- 
sidering the horn as a special black body cavity radiator emit- 
ting thermal radiation at ambient temperature. Transmitted 
to the outside, the radiation is normally unnoticeable. How- 
ever, if a water film or other reflective dielectric is present in 
front of the horn, part of the energy will be reflected back, 
some absorbed by the film, and the remainder will be trans- 
mitted outwards. If r is the fraction of incident power 
reflected, o: is the fraction of incident power absorbed, and 
t is the fraction of incident power transmitted, then r + a + 
1 . 


However, the thermal radiation power e emitted by the 
film (if in local thermodynamic equilibrium) is, according to 
the Kirchoff law of radiation, identical to the absorbed radia- 
tion OL, Therefore, if the horn and dielectric film are at ambient 
temperature Tq, the power (noise temperature) detected by 
the receiver for this effect will be: 

T = rjr + e) = T^ir^ot) = 

which is the same as the well-known formula generally used to 
determine the noise temperature contribution of any attenua- 
tion or mismatch loss. 

Obviously any signal received or transmitted will be atten- 
uated by 10 log / (dB’s). The power reflection and transmission 
coefficients for a dielectric film in air are (Ref. 8) 

[(1 - sin^ 

r = — 

(1 - ^ ^p^)^ + 4 v 4 ^p^sin^ (^ + x) 

_ [(l-p^)^+4p^ sin^ X] 

(1 - A^p^)^ + 4^^p^ sin^ (<I> + x) 

where 

A = exp (-2nd K) amplitude transmission coeffi- 

cient for single film transversing 

<p = IndnjX^ phase shift for single film transversing 

p = /? - l/« + 1 amplitude reflection coefficient for a 
single reflection 

X = tan”^ (InKln^ - 1) phase shift of single reflection 

d - film thickness 

n = refraction index 

= wavelength in air 
K = tgbjl 

tgb = a/cj€ loss tangent or dissipation factor 
a = dielectric conductivity 
e = dielectric constant 

The above expressions have been evaluated for the following 
cases of general interest: 

(Fig. 1) Temperature and attenuation of a water film in 
S-band 

(Fig. 2) Temperature and attenuation of a water film in 
X-band 
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(Fig. 3) Temperature and attenuation of snow and ice in 

S-band 

(Fig. 4) Temperature and attenuation of snow and ice in 

X-band 

These results should be considered indicative rather than 
totally accurate, because the values for the dielectric constant 
and dissipation factor at S- and X-band have been interpolated 
from the data presented for other frequencies in Table 4-19 of 
Ref. 2 (as the only source presently available). Nevertheless, 
it is considered that this questionable accuracy would not 
invalidate the following important conclusions: 

(1) A very thin water film on the Mylar will contribute to 
very high system noise temperatures, especially at 
X-band. This is mainly due to the high refraction index 
of water for these frequencies, thereby making it a 
good reflector, reducing also the effective wavelength 
(X = Xo/aj). 

(2) Figures 1 and 2 are evaluated up to very thick water 
films, not that they can be formed by rain, but in 
order to permit validation of the model by feasible 
experiments. In both cases, the puzzling result initially 
obtained of alternating maximum and minimum SNT 
is confirmed. This occurs at multiples of X/4 (wave- 
length in water) when the internal reflections of the 
film are in or out of phase with respect to the main 
reflection at the interface. 

(3) Figures 3 and 4 clearly show that SNT contributions 
for freshly fallen snow, hard packed snow, and ice are 
in any case much smaller than for liquid water. Subject 
contributions are a consequence of the lower refraction 
indexes: 1.1, 1.22 and 1.78, respectively, as compared 
to 8.83 and 6.5 for liquid water at S- and X-band. 

B, Effects of Water on Antenna Reflector Surface 

This could also be considered another potential contribu- 
tion to the general degradation during rain. However, it is 
difficult to confirm experimentally. Therefore, if the model 
used in the previous case is deemed valid, the reflection losses 
may be estimated using a similar expression. 

The power reflection coefficient for a dielectric film cover- 
ing a perfect reflector is (Ref. 8), using the same notation as 
before, 

^ AA'^ sin^ (0 - x/2) 

- (1 + pA - 4 p>l ^ sin^ (0 + x/2) 


Based on a temperature contribution of T = (1 - r'), the 

following values were estimated: 


Liquid water 

Film thickness, mm 

S-band temp, K 

X-band temp, K 

0.25 

2.6 

0.81 

0.50 

2.7 

7.4 

0.75 

4.3 

33.0 

Snow 

10 

0.005 

0.42 

20 

0.11 

0.58 


As can be seen, this potential contribution will be negligible 
when compared with similar films on the Mylar, as has been 
previously shown. 

It could be argued that the above models would not be 
applicable to real cases as a water film may not be formed. 
This is partially true because the water by surface tension may 
form dispersed drops instead of a film. Nevertheless, if the 
horn illumination is supposed to be uniform, and the tempera- 
ture contribution model practically linear for small thickness, 
the effect of the dispersed drops could be estimated assuming 
a film thickness equivalent to the total amount of water on 
the surface. 


III. Empirical Results 

Once the simplified working model was developed, a test, 
similar to the one mentioned in the introduction, was care- 
fully prepared for S-band at DSS 62. A flat plastic tray was 
fabricated having a slightly larger size than the horn window. 
With the antenna at zenith, the tray was perfectly leveled in 
front of the horn and barely touching the Mylar. Then each 
time a controlled amount of water was added, a measurement 
of the system noise temperature was taken with the noise 
adding radiometer. For the first two points obtained (equiv- 
alent to film thickness of 0.5 and 1 mm), the film was not 
actually formed, but the drops were purposely distributed 
more or less uniformly. 

All SNT data points have been directly plotted in Fig. 1, 
and the solid line curve is the theoretical water film con- 
tribution (not the total SNT including waveguide losses, 
receiver, etc.). The results, in general, are in good agreement 
with the simplified model, thereby confirming two aspects: 
First, that the model is valid for the study of these types of 
effects; second and more important, that the water film on the 
Mylar may have a catastrophic effect on telecommunications 
performance. 
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A similar test setup was prepared for X-band at DSS63. 
In this case, due to the geometry of antenna and X-band cone, 
the horn window could not be set horizontally and therefore 
not parallel to the levelled tray. This has some implication on 
the angle of reflection which may slightly jeopardize the test 
results. However, the SNT measurements, which are plotted 
in Fig. 2, are also in fairly good agreement with the model. As 
should be expected and can be seen from Figs. 1 and 2, the 
water temperature contribution of very thin films is much 
higher for X-band: 95 K vs 25 K for an equivalent film of 
one-tenth millimeter. 

No similar tests have been conducted to check Figs. 3 and 
4 (snow and ice) and therefore they should be taken only as 
indicative. 


IV. Suggestions for Degradation Reduction 

The possibility of an equivalent water film forming on the 
Mylar surface during rain will be directly related to two 
opposing factors: 

(1) The amount of water impinging on the surface which is 
a function of the rainfall rate direction and the antenna 
position. 

(2) The water drainage on the same surface which is also a 
function of the type of surface and its affinity to 
water, wind and rain intensity and direction, and 
antenna position. 

Although it is very difficult to quantify all these factors, one 
technique which will enhance the situation is to improve the 
water drainage as much as possible. In this regard, many dif- 
ferent ideas may be proposed. Some of our suggestions are 
presented below for further study, testing or development, 
based mainly on the X-band feed as the case of major concern. 

(1) The capability of Mylar to retain water drops seems to 
be low but there may exist another material with 
better performance, or a water repellent coating. 


(2) The ring connecting the Mylar to the X-band horn is 
quite thick. This could retain some water on the 
internal edge which may not be negligible when com- 
pared to the window size. 

(3) A special type of plastic shield wiper having low refrac- 
tion index and dissipation factor may be designed and 
tested. 

(4) A vibrating window (acoustically or ultrasonically 
activated) might be developed. 

(5) A centrifugal double window similar to the ones used 
in boats might be used. This method, if feasible, would 
be absolutely efficient in draining water (see sketch A 
of Fig. 5). 

(6) A well-designed system of nozzles (jets) blowing com- 
pressed air at high-speed over the Mylar window might 
be employed. 

As this system is relatively simple, it has been tentatively tested 
in a rudimentary form. A standard nozzle attached to a small 
plastic fixture (see sketch B of Fig. 5) was temporarily mounted 
near the X-band horn (outside its field of view) and connected 
by a hose to an air-compressor with a pressure of 90 psi. With 
a “rain emulator” (a plastic watering-can) a more or less uni- 
form shower was generated over the Mylar, obtaining a system 
noise temperature of about 1 80 K. Then with the nozzle blow- 
ing the SNT came down to about 100 K. The test was repeated 
with a lighter shower, changing the SNT from 85 K to 55 K 
when blowing. 

V. Conclusions 

From the theoretical and experimental results of this study, 
the logical conclusion would be to believe that the local rain 
effect may be of prime importance. However, to reach a firm 
conclusion, further tests have to be performed during actual 
rain in order to measure the relative contributions of the 
path and local effects and determine if they are separable. 
For a simple qualitative test, the rudimentary system suggested 
in IV-6 above has been temporarily left installed, but no test 
data is yet available. 
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Fig. 3. Temperature and attenuation of snow and Ice in S-band 
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Fig. 2. Temperature and attenuation of a water film in X-band 
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X-BAND HORN, SKETCH A: 
CENTRIFUGAL SUGGESTION 


Fig. 5. Proposed test 


Fig. 4. Temperature and attenuation of snow and ice in X-band 
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This article compares the capabilities and limitations, as well as the associated costs 
for two total energy systems for a diesel power generation plant Both systems utilize 
waste heat from engine cooling water and waste heat from exhaust gases. Pressurized 
water heat recovery system is simple in nature, requires no engine modifications, but 
operates at lower temperature ranges. On the other hand, a two-phase ebullient system 
operates the engine at constant temperature, provides higher temperature water or steam 
to the load, but is more expensive. 


I. Introduction 

Total energy systems have been widely studied as a more 
effective means of using the depletable energy sources. Waste 
heat from the diesel engine represents a potential source of 
energy that can be further utilized for heating and/or cooling 
purposes or, coupled with a Rankine bottoming cycle, to pro- 
duce extra power. Diesel engines for power generation have 
an efficiency between 30 to 35 percent, which means that 
65 to 70 percent is wasted to the atmosphere. 

First, about 40 percent of the total waste heat is exchanged 
during the combustion process in the engine, which tends to 
increase the cylinder wall temperature to such an extent that 
without cooling, the lubrication oil can evaporate resulting 
in engine damage. Water is usually employed to cool the 
engine. The temperature of the cooling water leaving the 
engines range from 87.7°C (190°F) to 104.4°C (220°F). 
However, to avoid excessive thermal stresses in the en^e, 
it is advisable to have a temperature difference of 5.55°C 
(10°F) to 8.33‘^C (15°F) between the inlet and the outlet 
water temperatures. The hot cooling water from the engine 
is then fed into an air-cooled or a water-cooled radiator 


before recirculating it back to the engine. Instead of reject- 
ing the heat to the atmosphere, as much as 90 to 100 percent 
can be recovered for heating or cooling purposes. Second, 
the other 60 percent of the waste heat is in the form of 
sensible heat in the exhaust gases at the end of the expan- 
sion stroke. The temperature of the exhaust gases varies 
from engine to engine, ranging between 371.1°C (700®F) 
and 698 .9°C (1200°F). The sensible heat in the exhaust 
gases can be extracted by means of a stack-type heat 
exchanger. To avoid water condensate in the exhaust ducts, 
the temperature of the exhaust gases leaving the stack cannot 
be lower than 176.7° (350°F), which means that only 60 to 
70 percent of the waste heat from exhaust gases is recoverable. 

DSN overseas stations in Spain and Australia operate their 
diesel engines continuously around the clock to provide power 
for the stations’ needs. Therefore, it is felt necessary to make 
use of the waste heat from the engines to conserve the valuable 
and depletable diesel fuel. DSS 62 is currently using the diesel 
waste heat for heating and cooling. DSS 61/63 is in the process 
of installing a waste heat recovery system for space heating, 
and an effort is underway to further use the diesel waste heat 
to provide cooling. 
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II. Relationship of Engine Load and the 
Available Waste Heat 

To utilize the waste heat from diesel engines for heating 
and/or cooling, it is important to investigate the engine operat- 
ing loads and the energy required for heating and/or cooling, 
since they determine the sufficiency or deficiency of the 
engine waste heat. By employing the waste heat for heating 
and cooling, the engine will be operating at a reduced load 
level since it no longer provides the electricity needed for the 
present heating and cooling devices such as electric boilers 
and vapor compression refrigeration units. In other words, 
the waste heat available will be less. On the other hand, the 
coefficient of performance (COP) of the absorption chiller 
and the efficiency of the heat exchangers should be taken 
into consideration to determine the heat needed to meet the 
heating and/or cooling load. For instance, a diesel engine 
operating at 500 kWe is providing heating and cooling loads 
of 150 kWt and 100 kWt, respectively. Utilizing the waste 
heat for heating and cooling, the engine load should be 
reduced by 200 kWe assuming an efficiency of 0.9 for the 
electric boiler and the COP of 3.0 for the vapor compression 
chiller. Further, taking the COP of the absorption chiller at 
0.6 and the heat exchanger efficiency at 0.9, the heat needed 
will be enhanced to approximately 350 kWt, and obviously, 
the waste heat from jacket water alone for an engine operating 
at 300 kWe load is not sufficient and the waste heat from 
exhaust gases will have to be used to compensate for the heat 
needed. Figure 1 shows the relationship between the engine 
load and the sufficiency of jacket waste heat for a one-day 
period. From hour 12 to hour 22, the waste heat from jacket 
water alone can no longer provide enough heat for both heat- 
ing and cooling; the heat available from exhaust gases or 
back-up systems should be further utilized to make up for the 
deficiency. 

In the event that the waste heat from exhaust gases is used 
to supplement the balance of heat required, several heat 
recovery systems can be used and their capabilities and limita- 
tions, as well as their associated costs, should be carefully 
investigated. 

III. Heat Recovery Systems 

Basically, there are two types of heat recovery systems 
that utilize both the engine jacket water and the exhaust 
gases: 1) a pressurized water system, and 2) an ebullient or 
2-phase steam/water, system. The pressurized hot-water sys- 
tem operates the engines at a temperature range of 82.2°C 
(180°F) to 92.9°C (210°F). One-third of the fuel energy 
input to the engine will dissipate into the jacket water, which 
would cause a temperature difference between the water 
coming into the engine and the water leaving the engines of - 


not more than 8.33®C (15°F). To achieve this, a temperature- 
controlled water flow rate should be maintained. The water 
leaving the engine is then fed to a gas-heat recovery unit to 
extract the waste heat from the exhaust gases. The gas-heat 
recovery units are typical gas-to-water heat exchangers with 
exhaust gases running in the hot side, and the engine jacket 
water in the cold side. Reference 1 hsts the common tempera- 
tures of exhaust gases for various models of caterpillar engines. 
Since the temperature of the exhaust gases leaving the heat 
recovery units should not be lower than 176.7°C (350°F), 
the amount of recoverable waste heat can easily be estimated. 
Moreover, the temperature of the water leaving the heat 
recovery units also can be determined. The temperature 
increase for the water across the heat recovery units will 
generally not exceed 5.55°C (10°F). As a result, a maximum 
water temperature of approximately 100°C (212°F) can be 
achieved only by using this scheme. For engines operating 
at temperatures close to 100°C (212°F), pressurized water 
systems should be employed throughout the cooling water 
loop to avoid flashing in the engines. It is also advisable to 
use in the external loop a heat exchanger to deliver the heat 
to the load to form a closed loop for the engine cooling water. 
A booster pump will also be installed since the existing engine 
pump will not be enough to circulate the water. At times 
when cooling is not needed, the heat can be removed using 
the existing radiators or converted to chilled water and stored 
in a chilled water storage tank. The pressurized hot -water sys- 
tem does not require any modifications to the existing engines. 
However, the maximum obtainable hot-water temperature 
(100°C) is lower than the absorption chiller manufacturer’s 
specification (116°C), which means an oversized chiller would 
be used at a lower capacity rather than the nominal size, even 
though it operates with the same coefficient of performance 
(COP). Table 1 lists the absorption chillers with the nominal 
sizes at 115.6°C (240°F) and the degraded capacities at 
90.5°C (195°F). Figure 2 depicts the schematic flow diagram 
for the pressurized hot -water system. 

The 2-phase steam/water ebullient system removes the heat 
from the engines by heat of vaporization. For an ebullient 
system, the engines are operating at a constant high tempera- 
ture ranging between 112.8°C (235°F) to 121. 1°C (250^F). 
The cooling water is fed to the engines at the operating tem- 
perature. As the water passes through the cylinders, the heat 
rejected will vaporize a portion of the water into wet steam. 
The wet steam is fed to heat recovery units which serve both 
as mufflers and as steam separators. Using the heat from 
the exhaust gases will further convert the wet steam to dry 
saturated or superheated steam, which feeds the absorption 
chiller. Theoretically, the coohng water into the engine should 
have the same temperature as the water leaving the engine, but 
a 1,1 °C (2°F) to 1.6®C (3°F) increase in temperature is 
normally experienced. Because of the smaller temperature 


90 



gradient across the engine, the engine is subject to smaller 
thermal stresses, which give the engines a longer life. Further, 
the cooling water circulation is by natural convection and a 
pump is not necessary. The water-steam flow rate is about 
1/70 of that of the pressurized system due to the large latent 
heat of water. However, this system requires a separate cooling 
loop for oil coolers, which uses water at 29.4°C (85°F). The 
cooling water loop should maintain a pressure of 69 kPa 
(10 psi) to 103.5 kPa (15 psi) to avoid flashing inside the 
engines; the temperature of the steam collected will be in the 
neighborhood of 115. 6°C (240°F) to 121. 1°C (250®F), which 
satisfies the chiller manufacturer’s specification. For engines 
running at such a high temperature, the seals, plumbing sys- 
tem, and even the cylinder head will have to be replaced for 
high-temperature applications. Thus, it is appropriate to 
convert the engines to the ebullient cooling system during a 
major overhaul period. When excess steam is not needed, 
the steam will be condensed in a bypass condensing unit and 
returned to the heat recovery units. Figure 3 shows the sche- 
matic diagram of a high-temperature ebullient system. Water 
treatment for the engine cooling water is necessary for an 
ebullient system, and it is critical to engine operation. It is 
desirable to have chemical testing of the water every other 
day to keep the water at a certain pH level. 

IV. Cost Comparison 

A pressurized hot-water system does not require any engine 
modification. However, the system needs a heat exchanger, 
a booster pump, and an oversized absorption chiller. On the 
other hand, an ebullient system requires.no heat exchanger, 
works with nominal size absorption chiller, but requires engine 
modification, an excess steam condensing unit, a holding tank 
for the condensate, and a more expensive heat recovery unit. 
It is necessary to have a thorough breakdown of the costs for 
each system to select the most cost-effective system. 

To demonstrate the cost comparison of the two systems, 
assume a diesel power generation plant that operates four 
diesel engines and utilizes the waste heat from the engines 
to feed an absorbing chiller of 700-kWt capacity. 


Table 2 lists the cost for each component for both systems. 
The installation cost and the maintenance cost for the two 
systems are assumed the same. The main differences in cost 
between the two systems fall into three categories: 1) engine 
modification, 2) heat recovery units, and 3) the absorption 
chillers. For ebullient systems, engines have to be modified 
for high-temperature application. The cost of the modifica- 
tion is taken as $7,000 per engine including the replacement 
of the cylinder head, seals and necessary changes in the plumb- 
ing system. The pressurized hot-water system uses a gas-water 
heat exchanger for the heat recovery unit, while the ebullient 
system employs a more expensive heat recovery unit acting 
as both a heat exchanger and a steam separator. Further, the 
pressurized hot-water system requires an oversized chiller to 
provide the same amount of cooling as the ebullient system. 
The material cost for hot-water system and the ebullient sys- 
tem are approximately $103K and $150K, respectively. 


V. Conclusion and Remarks 

A pressurized hot-water system and the high-temperature 
ebullient system are the two commonly used heat recovery 
systems associated with diesel total energy systems. Both 
systems recover waste heat from the engine jacket water as 
well as waste heat from the exhaust gases. The pressurized 
hot water system is simple, requires no engine modification, 
but needs an oversized absorption chiller. On the contrary, 
high-temperature ebullient systems require engine modifica- 
tions, more expensive heat recovery units and a water treat- 
ment system. The installation and maintenance costs for 
both systems are very much the same, but the material cost 
for the ebullient systems are 1 .5 times higher than that for the 
pressurized hot-water system. Furthermore, for future 
increase in station cooling consumption, the pressurized hot- 
water system has the capability of providing the additional 
cooling with a supplementary source such as solar energy, to 
enhance the hot water temperature to, in turn, increase the 
capacity of the absorption chiller, independent of the engine 
waste heat. 


Reference 
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Table 1. Degradation of chiller capacities and 
corresponding costs 


Nominal size 
at 115.6°C 
(240‘^F), kWt 

Operation capacity 
at 90.5°C 
(195°F),kWt 

Estimated 
cost, $K 

2180 

879 

70.0 

1976 

791 

64.5 

1835 ' 

738 

61.0 

1723 

703 

59.0 

1396 

562 

48.5 

1097 

440 

42.0 

819 

334 

36.0 


Table 2. Cost comparison for pressurized hot-water system and 
high-temperature ebullient system 


Component 

Pressurized 
hot-water system, 
$K 

High-temperature 
ebullient system, 
$K 

Engine modification 

- 

@ 7 X 4 = 28 

Heat recovery units and 
Insulation and controls 

@10 X 4 = 40 

@20 X 4 = 80 

Heat exchanger or 
condensing unit 

4 

10 

Absorption chiller and 
accessories 

59 

36 



Total Cost 103 154 
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Fig. 2. Pressurized hot-water system 



Fig. 3. High-temperature 2-phase water/steam ebullient system 
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Loadings of the 34-m HA-DEC Antenna for 
Deep Space Missions 
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The computer analysis of the 34-m HA-DEC antenna by the IDEAS program 
provided the rms distortions of the surface panels support points for full gravity 
loadings in the three directions of the basic coordinate system of the computer 
model. The rms distortions for the gravity vector not in line with any of the three 
basic directions were solved and contour plotted starting from three surface panels 
setting declination angle. By inspections of the plots, it was concluded that the 
setting or rigging angle of -15 degrees declination minimized the rms distortions 
for sky coverage of plus or minus 22 declination angles to 10 degrees of ground 
mask. 


I. Introduction 

A HA-DEC antenna has an axis of rotation parallel to the 
Earth’s axis of rotation combined with a second rotation 
axis that is normal to the first axis. Sky coverage is thus 
obtained by the despinning action of the first (polar) axis 
with the elevation or declination angle pointing provided by 
the rotation about the second (declination) axis. The polar 
axis rotation is measured in hour angle (HA) and the declina- 
tion axis rotation in plus or minus declination angle (DEC) 
from the equatorial plane of polar axis rotation. 

At an Earth’s location, the rotational motions of the 
antenna changes the direction of the gravity loading vector 
with respect to its structural symmetric plans and axis. This 
action results in variations of the reflective surface’s distor- 
tion values, after best fit of a paraboloid. 


II. Calculations 

By assembling a full or complete structural computer 
model of the 34-m reflector structure, the design and analysis 
were done in the IDEAS program (Ref. 1). The analysis 
results (Table 1) were the rms (root-mean-squared) distortions 
of the surface panels’ attach points with the applications of 
the full or 1.0-gravity loadings in the three directions of the 
cartesian coordinate system of the model. The cross-sectional 
areas of the reflector’s truss members were changed in the 
IDEAS program through its iterative steps of designing cycles 
to minimize the rms distortions. 

Using the method of analysis described in (Ref. 2), the rms 
distortion values of the 34-m HA-DEC antennas were plotted 
through the hour angle-declination angle ranges for several 
surface panels setting or rigging declination angle at zero hour 
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angle. Setting angles considered were at zenith look, at zero 
degrees declination, and at -15 degrees declination angles. 


Hi. Results 

The planetary missions normally confine the spacecraft 
to within plus or minus 22 degrees of declination angle values 
when the maximum gain capabilities of the ground antennas 
are required. By setting the surface panels at plus or minus 
15 degrees declination angle, with the sign dependent on the 
hemisphere location, the rms distortion was minimized for 
full mission sky coverage above 10 degrees from ground mask. 

This minimization is illustrated by Figs. 2,3, and 4 showing 
the rms distortions for DSS 12 with the surface panels setting 
angles of zenith look, zero degrees declination angle, and 
-15 degrees declination angle. Figures 5 and 6 illustrate the 
rms distortions with the setting at -15 degrees declination 
for DSS 61 and +15 degrees declination for DSS 42. 


Figure 1 describes the sign conventions used at the northern 
and southern hemisphere stations. 

The latitude values are shown in Table 1 with the full or 
1,0 gravity load values of the rms distortions for the three 
components loading directions as output by the IDEAS 
program. 


IV. Summary 

(1) The rms distortions quoted are only for the surface 
panels’ support points on the reflector structure. 

(2) For the selected observing ranges of plus or minus 
22 degrees declination angles and 10 degrees elevation 
ground mask, the maximum rms distortion is 0,4 mm 
for the surface panels set at plus or minus 15 degrees 
declination, as required for the hemispheric require- 
ment. 


References 
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Progress Report, Technical Report 32-1526, Vol. XII, pp. 100-111, December 15, 
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Table 1. RMS distortions for full-gravity loads 


Seq 

DSS 

Latitude, 

deg 

RMS distortion for 1.0 
gravity for component 
direction 

X, mm 

Y, mm 

Z, mm 

1 

12 (Goldstone) 

+35.20805 

0.307 

0.528 

0.648 

2 

42 (Australia) 

-35.21922 

0.307 

0.528 

0.648 

3 

61 (Spain) 

+40.23885 

0.307 

0.528 

0.648 



Fig. 1. Sign conventions - north and south hemisphere 



DECLINATION, deg ELEVATION, deg 


100 


100 



0 20 40 60 80 100 120 

HOUR ANGLE, deg 

Fig. 2. DSS 12 surface panels set at zenith look; panels 
supporting points; rms distortion contour map, mm 



0 20 40 60 80 100 120 

HOUR ANGLE, deg 


Fig. 3. DSS 12 surface panels set at 0 deg declination: panels 
supporting points; rms distortion contour map, mm 



HOUR ANGLE, deg 

Fig. 4. DSS 12 surface panels set at -15 deg declination; panels 
supporting points; rms distortion contour map, mm 



HOUR ANGLE, deg 

Fig. 5. DSS 42 surface panels set at -15 deg declination; panels 
supporting points; rms distortion contour map, mm 
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Radio frequency interference at NASA's deep space stations has become a serious 
problem. A radio spectrum surveillance capability at these stations is needed to determine 
the sources of interference so that preventative measures can be taken. The first phase of 
a program to develop this capability was the development of a low-cost surveillance 
station now in operation at the Goldstone Deep Space Communication Complex near 
Barstow, California. This interim surveillance station is described and findings from the 
use of this equipment are presented. 


I. Introduction 

In recent years there has been a significant increase in the 
number of occurrences of radio frequency interference (RFI) 
at the NASA deep space tracking stations. The interference in 
some cases degraded telemetry received from deep space 
spacecraft and in a few cases has caused the receivers to lose 
lock on the spacecraft carrier signal. Attempts to identify the 
sources of the interfering signals have been severely handi- 
capped because the stations have no means of determining the 
direction of the interfering source and very limited capability 
to view the interfering signal spectrum. Further, since the 
prime function of the station is to provide spacecraft data, the 
major activities of the station operators at times of signal 
degradation are to verify that a malfunction of station 
equipment has not occurred; to reacquire the spacecraft signal, 
if necessary; and to improve the quality of the spacecraft data. 

To facilitate the identification of RFI sources so that 
preventative measures can be taken, a radio spectrum surveil- 
lance station has been developed and installed at the Gold- 


stone Deep Space Communications Complex near Barstow, 
California. Because of the immediate need for RFI surveil- 
lance, a low-cost system using mostly off-the-shelf equipment 
was developed to provide surveillance capability. Later, a more 
sensitive and sophisticated system will be implemented. Its 
development will be influenced by the observations of this 
low-cost system. 

The purpose of this presentation is to describe the interim 
surveillance system and to present findings from field use of 
this equipment. 


II. Background 

Before describing the surveillance station, a brief discussion 
of the interference protection criteria for the deep space 
stations will point out the difficulty in providing them with 
interference protection. For brevity, the discussion is limited 
to the 64-m antenna stations. 
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Also some motivating factors for implementing the surveil- 
lance station and a general overview of the equipment cost and 
schedule are presented. 

A. Interference Criteria 

There are three major interference criteria for the deep 
space stations. These criteria were established to prevent 
receiver saturation and degradation of carrier tracking and 
telemetry performance. 

Power levels of - 120 dBW or stronger in the band of 2290 
to 2300 MHz will produce saturation of the DSN receiver. This 
condition causes an increase in the receiver-system noise 
temperature because of gain compression of the maser and can 
cause spurious signals to be introduced in the receiver because 
of operation in a nonlinear region of one or more stages of the 
receiver. Assuming that a signal source producing a - 1 20-dBW 
signal level at the receiver is 2 deg off the main beam axis of 
the station’s antenna, the power flux density (pfd) at the 
antenna corresponding to this power level is ~ 116 dBW/m^, 
using an off-axis antenna gain of: 

G = 32-25 log 0 1 < 0 < 48 

G = -10 0>48 

where 0 is the angle from the main beam axis in degrees and G 
is in dBi (Ref, 1). For the source having angles equal to or 
greater than 48 deg from the main beam axis, a pfd of -81 
dBW/m^ is required to saturate the receiver. 

The second criterion is a continuous wave (CW) signal no 
greater than -212 dBW in any 10-Hz band between 2290 and 
2300 MHz. This amplitude of interfering signal will induce 
phase modulation on the desired carrier signal of approxi- 
mately 10-deg amplitude, when the frequency separation 
between the two signals is comparable to, or less than, the 
carrier tracking loop bandwidth (assumed to be 10 Hz at 
receiver threshold). Using target locations of 2 deg and 48 deg 
off the antenna main beam axis, as above, the maximum 
allowable power flux densities are -208 dBW/m^ and -173 
dBW/m2. 

The third criterion is a power spectral density for noise-like 
interference no greater than -222 dBW/Hz. A -222-dBW/Hz 
amplitude interfering signal will cause approximately 1-dB 
degradation in telemetry performance. Both this and the 
second criterion are based upon C.C.I.R. Recommendation 
365-3 (Ref. 2) which assumes a 1 6-K system noise temperature 
and an antenna elevation angle of 30 deg. Corresponding 
power flux densities at 2 and 48 deg off pointing angles are 
-218dBW/m2-Hz and -183 dBW/m2-Hz. 


B. Motivating Factors 

At the time of considering the feasibility of a surveillance 
station at the Goldstone facility, there was evidence that 
aircraft in the vicinity of the stations produced the most 
detrimental radio frequency interference. Calculations show 
that an aircraft at a 160-km range and radiating 1-mW CW 
towards a deep space station can produce -145 dBW/m^ at the 
station’s antenna. Comparing this value with pfd values 
calculated from the interference protection criteria, it can be 
seen that emissions from aircraft can (and do) cause very 
serious interference. 

Of particular concern is the possibility that communica- 
tions with a spacecraft will be interrupted during critical 
operations such as a planetary encounter or trajectory change 
maneuver. At these times, it is very desirable to perform radio 
frequency surveillance of the spacecraft downlink frequency 
band. At the time of the decision to implement the surveil- 
lance station, the two upcoming encounters were the Pioneer- 
Venus spacecraft with Venus and the Voyager spacecraft with 
Jupiter. 

In July 1977, it was decided that, in the 18 months 
remaining before the Pioneer- Venus encounter, a surveillance 
station suitable for detecting aircraft emitters could be 
developed for reasonable cost. This station would provide 
protection against the major source of interference and yield 
radio environment information valuable for the development 
of future surveillance stations. 

C. Schedule and Cost 

The work on the development of the station began in July 

1977, and postinstallation testing was completed in November 

1978. Hardware costs were approximately $50,000, and 
slightly over 3 manyears of effort were expended. Much of the 
equipment such as power supplies, the equipment rack, and 
digital recorder were on hand and not purchased. 

The surveillance station did support both Pioneer-Venus 
spacecraft encounters and the Voyager 1 encounter. 


III. General Description 

The surveillance station is located on a mountain peak 
central to the three deep space stations and is normally 
unattended. In the unattended configuration, it monitors the 
frequency band of 2290 to 2300 MHz while scanning through 
360 deg of azimuth in approximately 1 min. The receiver 
sensitivity allows detection of power flux densities of -137 
dBW/m2 or greater at the antenna of the surveillance station. 
Upon detecting a signal, a photograph is made showing the 
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signal spectrum, spectrum analyzer control settings, time, and 
antenna azimuth. At the same time, the equipment provides a 
notification signal via landline to the Goldstone Radio 
Spectrum Coordinator. 

The surveillance station is normally operated continuously 
and serves a very useful function in aiding deep space station 
operators to determine if a radio frequency interferer is 
external or internal to their station. Upon occurrence of 
interference, the station operators may call the Radio Spec- 
trum Coordinator to determine if a signal was detected by the 
surveillance station. If so, the signal source is known to be 
external to the deep space station. If it was not detected by 
the surveillance station, however, it can not be concluded that 
the source is internal to the deep space station for it may be an 
external signal too weak to be detected by the surveillance 
station. 

The surveillance station can also be operated manually. In 
this mode the operator has full control of the spectrum 
analyzer settings, the antenna pointing direction, antenna 
rotation rate, and operation of the camera and digital recorder. 
The maximum frequency range of the equipment in this mode 
is 2250 to 2320 MHz and, by use of the narrowest resolution 
bandwidth (10 Hz) of the spectrum analyzer, sensitivites as 
high as - 1 77 dBW can be achieved. 

The manual mode is used primarily for real-time identifica- 
tion of interferers. Because the operator can determine the 
received signal characteristics and the direction of the source 
almost immediately, other agencies can be contacted while the 
source is still in the vicinity of Goldstone, thereby greatly 
facilitating the identification of the source. 

A third mode of operation, in which the spectrum analyzer 
is removed from the surveillance station and used with a 
hand-held antenna, is also provided. This mode is used for 
investigating radio frequency sources at locations other than 
the surveillance station site, for example, searching for radio 
emissions internal to one of the deep space stations. 

IV. Functional Description 

Figure 1 shows a functional block diagram of the surveil- 
lance station. The heart of the system is an analog spectrum 
analyzer (Hewlett Packard 8568A). Requisite features of this 
instrument are very high frequency stability, annotation of the 
display with important analyzer control settings and provision 
to display user-supplied information. The user-supplied infor- 
mation in this application is time, antenna-rotation rate, and 
antenna azimuth. The frequency range of the analyzer is 100 
Hz to 1.5 GHz, and down conversion of the S-band received 
signals is necessary. 


The down conversion is performed by the Microwave 
Assembly which contains the mixer and local oscillator 
circuits. To maintain high frequency stability, the local 
oscillator signal is generated by coherently multiplying a 
highly stable reference signal suppUed by the spectrum 
analyzer. 

The surveillance antenna is a standard gain horn having a 
maximum gain of 17.2 dBi. Its half power beamwidths are 28 
deg in azimuth (H plane) and 26 deg in elevation. The 
mounting bracket for the antenna allows the antenna elevation 
angle to be wrench-adjustable from -10 to +70 deg. Since the 
primary goal of this equipment is to detect the presence of 
aircraft emitters and sources at ground level, the antenna 
elevation is curretitly set to +5 deg. 

To achieve a useful receiver sensitivity at moderate cost, 
two uncooled, low-noise amplifiers [gallium arsenide field- 
effect transistor (GAs FET)] are used as a preamplifier to 
amplify the antenna output signal. The noise figure of the 
preamplifier is 1.2 dB. A bandpass filter is used before the 
'preamplifier to protect it from degradation by transmissions 
from a nearby deep space station and to restrict the number of 
mixer product terms out of the first mixer. A CW signal having 
a power level of -153 dBW (corresponding to a pfd of -138 
dBW/m^) will produce a barely discernible signal on the 
spectrum analyzer (12-dB signal-to-noise ratio) when the 
station is configured for unattended operation. The resolution 
noise bandwidth in this mode is 4400 Hz. 

Automatic recording of the received signal spectra is 
performed by a 35-mm camera triggered by a signal-to-noise 
ratio detector. Detecting the signal-to-noise ratio instead of the 
absolute signal level eliminates the problem of receiver gain 
variation with temperature (approximately 2 dB). 

The presence of a signal-to-noise ratio which exceeds a 
present threshold value, in addition to actuating the caniera, 
also: 

(1) Produces a printed record of the time and azimuth at 
the end of the sweep during which the signal was 
detected. 

(2) Notifies the Radio Spectrum Coordinator of the 
detection. 

(3) Increments counters at both the surveillance station 
and at the Radio Spectrum Coordinator’s position. 

Figure 2 is a photograph of the spectrum analyzer display 
taken by the recording camera. Time, antenna azimuth, and 
antenna rotation rate displayed at the top of the display are 
updated at the end of each spectrum analyzer sweep. With 
these parameters and the analyzer sweep time (displayed at the 
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end of the bottom line), the antenna azimuth at which the 
signal was received and the time of occurrence of any signal 
seen on the analyzer trace can be determined by interpolation. 
The Display Controller controls the timing of the equipment 
sequences, performs the antenna rotation rate calculation, and 
performs the data formatting necessary for communicating 
with the spectrum analyzer. 

The antenna motor is a DC motor with the field and 
armature separately excited. A switch is used in series with the 
armature to start, stop, or reverse the antenna rotation, and 
the motor speed is controlled by adjusting the armature power 
supply. 

An S‘Band signal source located 7.4 km from the surveil- 
lance station is provided for assessment of equipment perfor- 
mance. The source may be turned on from either the Radio 
Spectrum Coordinator’s position or the surveillance station 
thereby providing a test signal whose amplitude, direction, and 
frequency are accurately known. 

For transportable operation, the spectrum analyzer, the 
microwave assembly, and a hand-held antenna are used. The 
hand-held antenna assembly consists of a 1 2 dBi gain horn and 
a preamplifier having a noise figure of 3.1 dB. The preamplifier 
power supply is contained in the microwave assembly. 

The performance parameters of the equipment are listed in 
Table 1. 

Photographs of the major station elements are shown in 
Figs. 3 through 5. 

V. Operational Findings 

Since installation of the equipment in November 1978, 
there have been several hundred photographs taken by the 
equipment of radio spectra having frequency components in 
the deep-space-station downlink receive frequency band (2290 
to 2300 MHz), Although cataloging of these photographs is 
just beginning and many of the radio frequency sources have 
not yet been identified, some important facts have been 
learned. 

The results to date indicate that the number of detections 
correlate with military exercises in the vicinity of the 
Goldstone Deep Space Station Complex with the greatest 
number of detections associated with the use of FAA airspace 
R2508. During periods of critical deep space spacecraft 
operations, such as planetary encounters or spacecraft emer- 
gencies, military agencies are requested to impose radio silence 
on frequencies near the deep-space receive band and, during 
these periods, no sources have been detected. Further, it has 


been found that the weekends and evening hours, the hours 
not usually used for military exercises at the aircraft gunnery 
ranges, have the fewest number of detections. It should be 
pointed out, however, that, because of the large amount of 
effort required to identify potential interference sources, only 
sources which are considered as having a high RFI potential 
have been investigated. Of those investigated, approximately 
65 percent have been identified. The identified sources are all 
military aircraft. 

As an example of the benefit of the surveillance station at 
the Goldstone Deep Space Station Complex, it was learned by 
use of this equipment that a usual practice on board certain 
Air Force aircraft is to test some electronic equipment for a 
few minutes shortly after takeoff. The emissions from this 
equipment were first detected by the surveillance station at a 
range of approximately 170 km. From the surveillance-station 
photographs, the Air Force by very conscientious efforts 
identified the aircraft involved and modified operational 
procedure to eliminate this potential interference source. 


VI. Room for Improvement 

Although the unattended operating mode has satisfied the 
original requirements levied upon it to provide records 
showing the spectra and azimuth of potential interferers, there 
is a clear need to incorporate additional features which will 
facilitate the identification of the interferers. 

Identifying radio frequency sources is much like fighting a 
fire, the sooner you get to it the better are your chances. To 
solicit help from other agencies in identifying a radio 
frequency source, the Radio Spectrum Coordinator must have 
time, azimuth, and a description of the signal spectrum. With 
the present design, time is available to him, if he notes the 
time upon hearing the audible alarm. If he is not present at the 
time of detection, time is printed at the surveillance station 
(approximately 20 minutes away). Azimuth is not supplied to 
the Radio Spectrum Coordinator’s position, but this too is 
printed at the surveillance station. 

The most serious delay is that he must wait for develop- 
ment of the camera film to see the received signal spectrum. 
This can require from 6 hours to several days, depending on 
the quantity of film to be processed and to which facility it is 
sent for processing. 

A real-time display containing time, azimuth, and spectrum 
at the Radio Spectrum Coordinator’s position is the obvious 
solution. Since this position is not staffed on a 24-hour basis, 
the display device should be capable of making a permanent 
record. 
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Another deficiency is that the fixed spectruOi-analyzer 
control settings used in the unattended mode do not always 
provide an optimum display of the received signal spectrum. In 
several instances, the interfering signal carrier was outside the 
2290- to 2300-MHz range used and only sidebands falling 
within this range were photographed. Increasing the operating 
range results in lower frequency resolution and violates other 
requirements of the station. 

This problem can best be solved by providing adaptive 
control of the spectrum analyzer either by the Radio 
Spectrum Coordinator, by software routines in the Display 
Controller, or by both. 


VII. Conclusion 

The surveillance station described is providing highly useful 
experience in surveillance system requirements and techniques. 
It has identified several particular interferers and is yielding 
spectral signature data which, after cataloging, will serve as a 
library for rapid identification of frequently seen interference. 

It has helped deep-space station operators to determine 
whether or not RFI sources are external to their station. It has 
aided the designers of a more sophisticated surveillance 
station, employing fast Fourier transform techniques, in 
developing processing algorithms for their system. 
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Table 1. Performance parameters 




Operating Mode 



Parameter 

Unattended 

Manual 

Transportable 

Conditions 

Frequency 
coverage, MHz 

2290-2300 

2250-2320 

2250-2320 


Antenna 

rotation rate, RPM 


0-4 

N/A 


Receiver 

sensitivity, 

dBW 

<-153 

<-177 

<-156 

4.4-kHz resolution noise B.W. 
14-Hz resolution noise B.W. 
400-Hz resolution noise B.W. 

Detection 

sensitivity, 

dBW/m^ 

<-137 

<-164 

<-142 

90% confidence with 
above resolution 
band widths 

Source azimuth 
accuracy, deg 

±20 

at -135 dBW/m2 

±7 

at-140dBW/m2 

±12 

at 138 dBW/m^ 

95% confidence 
stationary target 

Frequency readout 
accuracy, MHz 

±0.2 

±0.06 

±0.06 


Amplitude 
accuracy, dB 

±5 

±4 

±4 


Time of detection 
accuracy, s 

±1 

±1 

N/A 


Permanent record 
capacity (records) 

>720 

>720 

N/A 

100-ft camera roll 

False-alarm 

rate 

< 1 per day 

N/A 

N/A 

SNR threshold setting 
12.5 dB 
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Fig. 1. Block diagram 
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Fig. 2. Recorded spectrum 



Fig. 3. Radio spectrum coordinator’s display 
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Fig. 4. Operator’s console 
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X-band noise temperature data have been taken at Goidstone DSS 13 continuously 
since August 1975. Presented here are sample data and cumulative distributions of 
atmospheric noise temperature increase above the quiescent baseline for the calendar 
years 1977 and 1978. Comparison is made with the existing Deep Space Network noise 
temperature statistics. 


I. Introduction 

Reference 1 presents statistics of zenith noise temperature 
increase above the quiescent (undisturbed) baseline at Gold- 
stone for the period August 1975 through January 1977. 
Included also are descriptions of elevation angle modelling and 
determination of atmospheric attenuation values from the 
given noise temperature increases. This article summarizes the 
X-band noise temperature statistics for the complete calendar 
years 1977 and 1978. 

II. System Description 

X-band atmospheric noise temperature data at DSS 1 3 (see 
Fig. la) are taken using an X-band noise-adding radiometer 
and horn antenna (no reflector) vyith a beamwidth of 15 
degrees. Numerous other instruments at the site provide 
measurements of rainfall, temperature, relative humidity, 
atmospheric pressure, and solar radiation. Data points are 


recorded every 2 minutes. The X-band horn is pointed at an 
elevation angle of 30 degrees and an azimuth of 325-degrees. A 
complete description of the X-band radiometer system is given 
in Ref. 2. 

The rain gauge at DSS 13 is an MRI (Meteorological 
Research, Inc.) tipping bucket unit. A funnel collects rain, 
which is channelled to the tipping bucket. The bucket tips 
only after it has collected the equivalent of 1 /100-inch (0.25 
mm) of rain. Each tip empties the bucket and registers as a 
step on the cumulative rainfall data (see Fig. lb). Avery light 
rain could continue for an extended period and only be 
registered as a single step (after 0.01 -inch had fallen) in the 
rain gauge data. Thus, the rain gauge gives inaccurate rain rate 
values for very light rain. The cumulative amount vyill be 
correct, however, to within 0.01 -inches. Heavy rains (>10 
mm/hr) will cause the bucket to empty more than once during 
one 2-minute data period. There is little loss of accuracy at 
these large rates. 
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III. Noise Temperature and Rainfall Data 

Figure 1 shows a five-day segment of raw antenna noise 
temperature and cumulative rainfall data taken during January 
1978. This is typical of a rainy period of the year. The 
characteristic “spiked” rain signature corresponds almost 
perfectly to the cumulative rainfall data shown directly below. 
The intensity of rain, or rain rate (mm/hr), can be determined 
from the slope of the cumulative rainfall curve. During the 
period 18-19 hours on day 16 the rain rate is nearly constant 
mm/hr), yet the noise temperature exhibits erratic 
changes. Possible causes of this could be alternate “beading” 
and “running-off” of raindrops on the antenna horn aperture, 
cover, or the radiometer “seeing” rain which is not actually 
falling into the rain gauge. Further investigation of rain and its 
effects on microwave systems will be carried out. Correlation 
of rain and noise temperature effects and the development of a 
weather/microwave model will be the subject of a future DSN 
Progress Report article. Figure 2 shows the correspondence of 
rain rate (slope of curve in Fig. lb) and increased antenna 
noise temperature for day 16. 

As described in Ref. 1, the noise temperature measured at 
an elevation angle other than zenith is referred to zenith for 
purposes of statistical consistency. Note in Fig. 1 that the 
baseline system antenna temperature is approximately 9 
kelvins for 30-deg elevation angle. This includes the minimum 
clear sky contributions from water vapor and oxygen and the 
noise temperature of the ground, cosmic background, and 
waveguide. Since the measurement of interest is the zenith 
increase above baseline due to atmospheric variables, the 
baseline is subtracted from the antenna temperature and the 
remainder divided by 1/sine (elevation angle) as a “weather- 
mass” correction. As an example, Fig. 3 shows the zenith 
equivalence (above baseline) of the data in Fig. 1 from Day 1 6, 
08 hours to Day 17, 08 hours. Times in both figures are 
referenced to local midnight (PST). The 67 kelvin peak at 30 
deg elevation in Fig. 1 becomes a 29 kelvin peak above 
baseline at zenith [(67-9) 2 = 29]. Nine-kelvin antenna 

temperatures in Fig. 1 become equivalent zenith zero-kelvin 
increases above quiescent baseline in Fig. 3. 

Figure 1 also shows typical clear weather (Day 15, 22 
hours) and cloudy weather (Day 14, 06 hours) noise tempera- 
ture measurements. The small noise temperature “humps” 
occurring at 10-14 hours were probably due to noise diode 
thermal effects and have since been eliminated and do not 
appear in the noise temperature statistics. 

IV. Noise Temperature Statistics 

Statistics have been generated for equivalent zenith X-band 
atmospheric noise temperature increase above quiescent base- 


line for various year-quarters and day-quarters at Goldstone. 
“Equivalent” means that the data are taken at elevation angles 
other than zenith and are converted to zenith values as 
described above. Statistics for system noise temperature can be 
obtained by adding the baseline system noise temperature for 
a particular antenna and elevation angle. Thus, for DSS 14, the 
64-meter antenna at Goldstone, approximately 25 K ,(at 
zenith) must be added to the atmospheric increase to obtain 
total system noise temperature. For example, if the atmos- 
pheric noise temperature increase at zenith is 5 K or less, 90% 
of the time, the total system noise temperature at zenith 
would be 30 K or less, 90% of the time. Different antenna 
elevation angles contribute different ground and clear sky 
atmospheric contributions and thus the baseline system noise 
temperature will be a function of elevation angle. A further 
discussion of elevation angle modelling and attenuation calcu- 
lation is given in Ref. 1 . 

Figure 4 shows the cumulative distribution of zenith noise 
temperature increase for a very wet period in 1978 (first 
year-quarter, January-March) and for a very dry period in 
1977 (second year-quarter, April-June). These curves have 1 K 
resolution. As an example for reading the curves, for the wet 
quarter, 95% of the time the zenith noise temperature increase 
was less than 5 K above the baseline. The wet/dry comparison 
shows that about 3 times as much data lie above 5 K in the 
wet quarter as in the dry quarter. Some amount of time the 
radiometer records values less than the accepted baseline value. 
This is due primarily to radiometer drift about the “zero-K” 
point. The baseline could alternatively be defined as that value 
below which only 5% or 1% of the data lie. 

Figure 5 shows a comparison of distribution functions for 
total-year Goldstone data in 1977 and 1978. In 1978, about 
70% more data lie above 5 K than in 1977. Indeed, 1978 was a 
much wetter year than 1977. Also shown in Fig. 5 is the 
total-year Goldstone noise temperature distribution function 
from the DSN Flight Project Interface Design Handbook. The 
data used to generate the Fig. 5 curve are pre-1970 and it is 
not known how these data compare descriptively (wet or dry) 
to the 1977 (very dry) and 1978 (relatively wet) data, 
although the curve indicates substantially worse weather than 
in either 1977 or 1978. As has been noted before (Ref. 3) the 
data appear to be “pessimistic,” i.e., they indicate that the 
weather effects are far worse than they appear to be based on 
actual measurements. 

Since a presentation of distributions for every year/day- 
quarter for 1977 and 1978 would be prohibitively lengthy, 
only tables of year-quarter and year-total cumulative distribu- 
tions are presented. Tables 1 and 2 present the 1977 and 1978 
Goldstone cumulative distributions, respectively, up to the 
1 0 K and/or 99% level. 
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Table 1. Cumulative distribution of zenith atmospheric noise temperature increases above 
Quiescent Baseline at Goldstone, 1977 


TZX» 

1st year quarter 

2nd year quarter 

3rd year quarter 

4th year quarter 

Year total 

0 

.487 

.464 

.061 

.017 

.300 

1 

.995 

.955 

.985 

.924 

.956 

2 

.999 

.980 

.992 

.951 

.976 

3 

.999 

.984 

.992 

.959 

.980 

4 

1.000 

.985 

.992 

.965 

.982 

5 

1.000 

.986 

.992 

.970 

.984 

6 

1.000 

.988 

.992 

.975 

.987 

7 

1.000 

.989 

.992 

.978 

.988 

8 

1.000 

.990 

.993 

.982 

.990 

9 

1.000 

.991 

.993 

.985 

.991 

10 

1.000 

.991 

.993 

.988 

.992 


^TZX = X-band zenith atmospheric noise temperature increase above quiescent baseline, K; e.g., 2nd year quarter, TZX = 3, 98.4% of data below 
3 K increase 


Table 2. Cumulative distribution of zenith atmospheric noise temperature Increases above 
quiescent baseline at Goldstone, 1978 


TZX® 

1st year quarter 

2nd year quarter 

3rd year quarter 

4th year quarter 

Year total 

0 

.207 

.053 

.694 

.512 

.216 

1 

.735 

.332 

.783 

.964 

.712 

2 

.879 

.947 

.945 

.980 

.931 

3 

.914 

.981 

.981 

.986 

.959 

4 

.933 

.983 

.984 

.988 

.968 

5 

.949 

.986 

.987 

.990 

.975 

6 

.962 

.988 

.989 

.992 

.980 

7 

.971 

.990 

.992 

.993 

.985 

8 

.979 

.992 

.994 

.993 

-988 

9 

.984 

.994 

.995 

,994 

.991 

10 

.987 

.995 

.996 

,995 

.992 


^TZX = X-band zenith atmospheric noise temperature increase above quiescent baseline, K; e.g., year total, TZX = 5, 97.5% of data below 5K 
increase. 
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ANTENNA TEMPERATURE, K RAIN RATE, mm/hr 



Fig. 2. X-band antenna temperature at Goldstone DSS 13 and 
corresponding rain rate, January 1978 



Fig. 3. Zenith equivalent X-band atmospheric noise temperature 
increase above baseline, K, corresponding to Figs. 1 and 2, 
January 1978 
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CUMULATIVE DISTRIBUTION 



ZENITH ATMOSPHERIC NOISE TEMPERATURE 
INCREASE ABOVE BASELINE, K 

Fig. 4. Wet-quarter/dry-quarter atmospheric noise temperature 
statistics at Goidstone DSS 13 
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CUMULATIVE DISTRIBUTION 



ZENITH ATMOSPHERIC NOISE TEMPERATURE 
INCREASE ABOVE BASELINE, K 


Fig. 5. 1977/1978 atmospheric noise temperature statistics at 
Goidstone DSS 13 and existing 810-5 distribution 
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Fig. 1. Raw noise temperature and rainfall data from Goldstone, January 1978, (upper row) X-band antenna noise temperature, K, 
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30-degree elevation angle, (lower row) cumulative rainfall, mm 
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